{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "480765be",
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "import HAWCtoHerschel as HtoH"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d5477674",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Source information\n",
    "SourceName = 'Snake'\n",
    "FullName = 'Filament 6'\n",
    "# Loading the fits files\n",
    "SourceT = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Tdust_10_10.fits')\n",
    "SourceN = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_N_10_10.fits')\n",
    "SourcePol = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Gal.fits')\n",
    "SourcePlanck = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Planck.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8fbf97a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot parameters\n",
    "ScaleLength = 0.0197572 # 0.0190986 # 1 pc in degrees for the target\n",
    "\n",
    "# Region location, width, and height\n",
    "l_center = 11.1309987 # Position in galactic longitude (degrees)\n",
    "b_center = -0.0979712 # Position in galactic latitude (degrees)\n",
    "l_width = 0.4241204 # Width in galactic longitude (degrees)\n",
    "b_height = 0.1803147 # Height in galactic latitude (degrees)\n",
    "Region = [l_center,b_center,l_width,b_height]\n",
    "\n",
    "# Intensity contours for the plot\n",
    "Icontours = [10, 20, 30, 40, 50, 60] # in mJy/arcsec^2\n",
    "# Ranges for plots\n",
    "Nscale = [0,2.0e+22] # Range for column density in cm^-2\n",
    "Tscale = [15.0,25.0] # Range for temperature map in K\n",
    "Iscale = [0.0,20.0] # Range for intensity map in mJy/arcsec^2\n",
    "# Vector scaling fraction in polarization plots\n",
    "Pvscale = 0.6\n",
    "\n",
    "# Figure size\n",
    "MapFigSize = [10.67,6] # in inches\n",
    "\n",
    "# Effective HAWC+ beam size to use\n",
    "HAWCBeam = 0.00519 # 18.7'' in Band E"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a2a163b6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Beam Size of polarization data\n",
      "18.2\n",
      "Beam size of target pojection\n",
      "36.395999999999994\n",
      "Pixel size of the polarization data\n",
      "4.550000000000001\n",
      "Pixel size of the target projection\n",
      "14.0000004\n",
      "Target standard deviation for smoothing in arcseconds\n",
      "13.384761981071023\n",
      "Target standard deviation for smoothing in pixels\n",
      "2.941705929905719\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: nan_treatment='interpolate', however, NaN values detected post convolution. A contiguous region of NaN values, larger than the kernel size, are present in the input array. Increase the kernel size to avoid this. [astropy.convolution.convolve]\n",
      "WARNING: FITSFixedWarning: 'datfix' made the change 'Set MJD-OBS to 58312.400043 from DATE-OBS'. [astropy.wcs.wcs]\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:96: RuntimeWarning: invalid value encountered in divide\n",
      "  dPI = ((StokesQ*dQ)**2.0 + (StokesU*dU)**2.0)**0.5/PI_biased\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:106: RuntimeWarning: invalid value encountered in divide\n",
      "  P_biased = 100.0*PI_biased/StokesI\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:109: RuntimeWarning: invalid value encountered in divide\n",
      "  dP = np.absolute(P_biased*((dPI/PI_biased)**2.0 + (dI/StokesI)**2.0)**0.5)\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:120: RuntimeWarning: invalid value encountered in divide\n",
      "  dO = (0.5*180.0/math.pi)*((StokesQ*dU)**2.0 +\n"
     ]
    }
   ],
   "source": [
    "# Reprojecting the HAWC+ Stokes data to the Herschel pixel scale\n",
    "NewPol = HtoH.ReprojectCube(SourcePol,SourceN, Smooth=True)\n",
    "#NewPol.writeto('test2.fits', overwrite=True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "31eb9087",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:228: RuntimeWarning: invalid value encountered in divide\n",
      "  imask_01 = np.where(Pdata[0].data/Pdata[1].data < IdI) # Total intensity SNR threshold\n"
     ]
    }
   ],
   "source": [
    "# Creating downsampled vector catalogs (independent vectors)\n",
    "# Creating the mask for the region\n",
    "DownsampleMask = HtoH.DownsampleVectors(NewPol,Step=2)\n",
    "# Creating the new catalogs for the region\n",
    "CatI, CatdI, CatQ, CatdQ, CatU, CatdU, CatP, CatdP, CatO, CatB,CatdO, CatPI, CatdPI, CatN, CatT = HtoH.MakePCats(SourceN, SourceT, SourcePol, NewPol, CatMask=DownsampleMask)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2c4fb0cb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "===\n",
      "Curve fit test\n",
      "Log scale fit: (-0.9408517995392947,47.67665403725213)\n",
      "Direct power law fit: (-0.9551343367378098,50.82902702393519)\n",
      "===\n",
      "\n",
      "Power law index: -0.9408517995392947 ± 0.03282108263542182\n",
      "Coefficient: 47.67665403725213 ± 3.722374462247168\n",
      "Chi-Squared: 295.6432534160674\n",
      "Number of elements: 146\n",
      "Reduced Chi-Squared: 2.0530781487226903\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfUAAAF2CAYAAAB3bAFnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB1tUlEQVR4nO3dd1xT1/sH8E8S9lSGQERQUUFFQala6/w5a62r1lKrVYZaLYJKna1KxVVt60CtVougdaF+1VZr1UpdtU4UNyiKgoIgslEZyfn9EZMSCJCEhAye9+t1X4Sbe899biB5cs499xwOY4yBEEIIITqPq+kACCGEEKIalNQJIYQQPUFJnRBCCNETlNQJIYQQPUFJnRBCCNETlNQJIYQQPUFJnRBCCNETlNQJIYQQPUFJnRBCCNETlNTLKS4uxrfffovi4mJNh0JInaP/f1Lf6cN7QCuS+saNG9G+fXtYWVnBysoKXbt2xZ9//il5/s2bNwgKCoKtrS0sLCwwcuRIZGRkqDyO4uJiLFq0SKf/oIQoi/7/SX2nD+8BrUjqzs7O+O677xAXF4erV6+iT58+GDZsGO7cuQMAmDFjBg4fPox9+/bhzJkzSEtLw0cffaTwcTZs2KDq0HWWtrwW6o5DVeXXthxl9ldkH3m31Za/u6Zpy+tQF3Fow3tA3f//8m6vLX93tWJaqmHDhuyXX35hubm5zNDQkO3bt0/y3L179xgAduHCBYXKbN26dbXP5+XlMQAsLy9PqZh1SU2vRV1RdxyqKr+25SizvyL7yLttddvR/3/dq4s4tOE9oO7/f3m3rw85wECzXykqEwgE2LdvH4qKitC1a1fExcWhtLQU/fr1k2zj4eEBFxcXXLhwAe+++67McoqLi6WaUIRCIYqLi5GXlwcOhyNzn/z8fKmf+kwgEGjFeao7DlWVX9tylNlfkX3k3ba67ej/Xz/j0Ib3gLr//+XdvqZttOU9wBjDy5cvYWNjAy73vwZ1Y2NjGBsb17izVrh58yYzNzdnPB6PWVtbsz/++IMxxtjOnTuZkZFRpe07derEZs+eXWV5YWFhDAAttNBCCy206MUSFhZWYy7Vmpq6u7s74uPjkZeXh/3792P8+PE4c+aM0uXNmzcPoaGhkt/z8vLg4uKC1NRUWFlZqSJkzZgyBdi1S/Tzu+9qVVRZWRmePXsGW1tb3Lt3D+3atYOJiYlkfePGjWFgoDX/IoQQUi/k5+ejSZMmSElJgbW1tWR9jbV0AFrziW1kZIQWLVoAAHx8fHDlyhWsXbsWvr6+KCkpQW5uLho0aCDZPiMjA46OjlWWV1UzhbiHvc76+GNRUj92DNiwAajiUoI8ysrKkJ+fDysrK5ibm8PKykqS1MXrKakTQohmWFtbK5yvtKL3uyzia+A+Pj4wNDREbGys5LnExESkpKSga9euGoxQQ/r3B4yNgeRk4O5dTUdDCCFEi2hFNWzevHkYNGgQXFxcUFBQgF27duH06dM4fvw4rK2tERgYiNDQUNjY2MDKygrBwcHo2rVrlZ3k9JqFBdC3L3D0KPD770DbtpqOiBBCiJbQiqSemZmJcePGIT09HdbW1mjfvj2OHz+O/v37AwBWr14NLpeLkSNHori4GAMHDsRPP/2k4ag1aOjQ/5L6vHmajoYQQoiW4DDGmKaDqAv5+fmwtrZGXl6ebl9TB4C0NKBxY9H19LQ0oJq+BdUpKyvD06dPYWdnh9u3b8Pb21tyTf3p06dwdnama+qEEFLHapOvtPaaOqkGnw906gQwBvzxh6ajIYQQoiUoqeuqoUNFP3/7TbNxEEII0RqU1HWVOKn/9Rfw6pVmYyGEEKIVKKnrqnbtgKZNgTdvgOPHNR0NIYQQLUBJXVdxOMCIEaLHBw5oNhZCCCFagZK6Lhs5UvTz8GGgpESzsRBCCNE4Suq6rGtX0e1seXnA339rOhpCCCEaRkldl3G5/zXB/+9/mo2FEEKIxlFS13XiJvhDhwCBQKOhEEII0SxK6rquZ0/AxgbIygLOndN0NIQQQjSIkrquMzQEhg0TPaZe8IQQUq9RUtcHH30k+nngACAU1rq4tLQ0rF27FmlpabUuixBCSN2hpK4P+vUDLC2BZ8+AK1dqXVxkZCQSEhIQGRmpguAIIYTUFUrq+sDEBBg8WPRYBb3gAwMD4eHhgcDAwFqXRQghpO5QUtcX4l7w//ufaPa2WuDz+Zg2bRr4fL4KAiOEEFJXKKnri0GDRDX2R4+Amzc1HQ0hhBANoKSuL8zNgfffFz2mgWgIIaReoqSuT8o3wddAKBQiNzcXwlr0lldFGYQQ7RYdHQ0OhwMOhyNZ17t3b3A4HPj5+WkuMCITJXV98uGHgJERcPeuaKmGUChEXl4eGGNIT09HZGSkwrewUVInRLeJk7Os5dChQwAAe3t7dOnSBV26dNFssCrk5+cHDoeD3r17y73PoUOH0LNnT1haWsLU1BQtW7bEd999p74glURJXZ80aAAMHCh6HBMj927r1q3DyZMnERERoZ64CCFazcjISJK4xYuNjQ0AYPDgwbh48SIuXryo4Sg158cff8SIESNw7tw5WFhYoHXr1nj9+jViY2M1HVollNT1ja+v6GdMTK17wRNC6gcnJydJ4hYvPXv2BCC7+V2W8ePHo2XLlrC0tISRkRFcXV0REhKC/Px8yTbiGnLTpk2xbds2NGnSBFZWVpg+fTpev36N6dOnw9raGi4uLti0aZNU+WlpaQgICACfz4eRkRGaN2+OxYsXo6ysTLKNuOVh3LhxCAsLg5OTExo2bIixY8eioKAAACTHBoAzZ85Izu306dMyzys1NRVz584FAERERCAtLQ3Xrl3D06dPcUALR/GkpK5vhg4V9YJPTARu3JBrl+DgYPTr1w8hISFqDo4Qoq9+++035OTkwM3NDU2aNEFKSgrWrVsnc7yL9PR0fPnllzA2NkZBQQHWrl0LHx8fREVFwdraGqmpqQgKCsK9e/cAAC9fvsS7776LqKgoFBYWonXr1khNTcXChQsxadKkSuXv2bMHq1evhqmpKXJzc7Fz505JU3mHDh1gZ2cHALC0tJS0TFhZWck8rwMHDqCsrAzm5ua4ePEi7Ozs4OTkhM8//xxFRUWqevlUhpK6vrG0BD74QPRYziZ4JycnBAYG0n3phNRTT548qXRNXVFnzpxBVlYW4uPj8fDhQ3zzzTcARNei37x5I7VtSUkJTpw4gfv376NJkyYAgKSkJNy4cQN3796FiYkJhEIhzpw5AwBYv349UlNT4eDggIcPH+LGjRvYv38/AFFLQlJSklT5JiYmuHfvHpKSkuDj4wMAkqbygwcPYvDbwbo6duwoaZno2LGjzPNKTEwEABQVFWHfvn1wcnLCy5cvsWPHDnzwwQcoLS1V+LVSJ0rq+oia4AkhCpB1TV1RJ0+ehKenJ0xNTcHhcLB06VIAQFlZGV68eCG1bcOGDdGtWzdwuVy4uLgAADw9PdG0aVNYWFigUaNGAICMjAwAwOXLlyW/N2rUCBwOB8OHDwcAMMZw6dIlqfL79OmDxo0bg8vlwsPDQ6osRZVv3t+6dStu374tGUL7+vXrOH/+vFLlqouBpgMgajB4MGBmBiQnA1evAp06aToiQogWE19TV9bOnTsxc+ZMSVlNmjRBVlYWHj16BAAQCARS25dv6jYwMKi0TtxSwCpUSiwtLdGmTZtKxzczM5P6vUGDBpXKr1iWvBo3bix53OntZ2nnzp0l6x4/fqxUuepCNXV9ZG4ODBkieqxAL3hCCFGG+AuBpaUlkpOTcenSJQwYMEBl5YuTqYGBAfbs2SNpMv/rr7/w5ZdfYsSIEQqVJ/4SIM818X79+kkeX716VeonALRs2VKhY6sbJXV9JW6C37tXJdOxEkJIVdq3bw8AKCgoQPPmzdG8eXPs3btXZeUHBQWhcePGyMnJgbu7O7y9veHm5gZbW1uMHz9e4fLETfJXr15Fu3bt8O677+L169cyt+3WrRuGDRsGAPD390e7du3g7+8PAOjbty+6deum5FmpByV1fTVokKjTXGoqUI/vLyWEqF9gYCBCQ0NhZ2eHgoIC9O7dG+Hh4Sor397eHhcvXoS/vz9sbW1x584dvH79Gj169MDq1asVLi8gIAAjR46EtbU1bt++jUuXLlW6RFDenj17MGfOHDg6OuLBgwdo1qwZFixYgMOHD9fmtNSCw5S90KBj8vPzYW1tjby8vCpvXdA7n38O7NgBhIQAa9dKPVVWVobLly8jJiYGU6dOxcuXL+Ht7Q0TExOUlZXh6dOncHZ2llyPkkXe7QghhMivNvmKaur6TNwEv28fIONbaExMDBITExEdHV23cRFCCFELSur6bMAA0dCx6enAP/9UetrX1xfu7u40KQMhhOgJSur6zMgIEPcK3bWr0tMODg4ICQmBk5NTHQdGCCFEHSip67vPPhP93LcPKCnRbCyEEELUipK6vvu//wOcnICcHODPPzUdDSGEEDWipK7veDxg9GjR4507NRsLIYQQtaKkXh+MHSv6+fvvQF6eZmMhhBCiNpTU6wNvb6B1a6C4GNDC+X8JIYSoBiX1+oDD+a+2vmOHZmMhhBCiNlqR1JcvX45OnTrB0tISjRo1wvDhwyVz2Ir17t270ny/kydP1lDEOkjcC/7UKeDZM83GQgjRWhs2bEDTpk1hYmKCLl26SKY9rUpBQQGmT58OV1dXmJqa4r333sOVK1eq3P67774Dh8PB9OnTVRy5iKLxK7qPuuOvLa1I6mfOnEFQUJBk1p3S0lIMGDCg0gw6EydORHp6umRZuXKlhiLWQU2bAt27i+ZX37NH09EQQrRQTEwMQkNDERYWhmvXrsHLywsDBw5EZmZmlftMmDABf/31F3799VfcunULAwYMQL9+/fBMRuXhypUr+PnnnyUTwMijd+/eco96qUz8iuyjTPx1jmmhzMxMBoCdOXNGsq5Xr15s2rRpSpeZl5fHALC8vDwVRKijNm5kDGDM25uVlpaypKQklpSUxAoKCtiFCxfY69evGWOMlZaWsuTkZFZaWlptcfJuR0h9d+nSJdarVy9mYmLC3N3d2ZUrV9jPP//MhgwZounQpHTu3JkFBQVJfhcIBIzP57Ply5fL3P7Vq1eMx+OxI0eOSK3v2LEj++abb6TWFRQUsJYtW7K//vpLoc/zXr16saioKLXEr8g+ysavjNrkK62oqVeU97aHto2NjdT6nTt3ws7ODp6enpg3bx5evXpVZRnFxcXIz8+XWuq9UaMAQ0MgPh64e1fT0RBSL1y8eBG9evXC4MGDcfPmTbRu3Rrh4eFYsWIFFi1aJLXtsmXLYGFhUe2SkpKiljhLSkoQFxcnNX84l8tFv379cOHCBZn7lJWVQSAQwMTERGq9qakp/qkwNHVQUBAGDx4sVb4qKRO/IvuoO35ZKuaw4uLiGvfRuqm1hEIhpk+fjm7dusHT01Oy/rPPPoOrqyv4fD5u3ryJOXPmIDExEQeq6M29fPnySm+Yes/WVjQl6++/g7t7NzBhgqYjIkTvhYaGYtSoUZg1axYAYPTo0Rg9ejSGDRuGDh06SG07efJkfPLJJ9WWx+fz1RJnVlYWBAIBHBwcpNY7ODggISFB5j6Wlpbo2rUrFi9ejNatW8PBwQG7d+/GhQsX0KJFC8l2e/bswbVr16q91i62bNkyLFu2TPL769evcfHiRUydOlWy7u7du3Bxcal1/PLuo0j8qtSkSROp38PCwvDtt99Wu4/WJfWgoCDcvn270re8SZMmSR63a9cOTk5O6Nu3Lx4+fAg3N7dK5cybNw+hoaGS3/Pz8yu9QPXS2LHA77+Ds3MnEBAAcLWysYYQvfD06VNcuHABP/zwg2SdgYEBGGMyKx02NjaVWigVNXfuXKxYsaLabe7duwcPD49aHUfs119/RUBAABo3bgwej4eOHTti9OjRiIuLAwCkpqZi2rRp+OuvvyrV6GWp+MVmzJgxGDlyJD766CPJOnV9sZFF0fhVfezyU68aGxvXuI9WJfWpU6fiyJEjOHv2LJydnavdtkuXLgCApKQkmUnd2NhYrheg3hkyBGjQAJzUVJhevIjX772n6YgI0Vv37t0DAHTs2FGyLjExEZ07d0a7du0qbV+xliqLrFpqeV999VWNMy82b9680jo7OzvweDxkZGRIrc/IyICjo2OVZbm5ueHMmTMoKipCfn4+nJyc4OvrKzlGXFwcMjMzpV4DgUCAs2fPYv369SguLgaPx5M8V/GLjampKRo1aiRV85dFmfjl2UfR+FXJyspK4fnUtSKpM8YQHByMgwcP4vTp02jWrFmN+8THxwMAzTCmKBMT4NNPgU2bYHngACV1QtQoLy8PPB4PHA4HAJCdnY0ffvgBXl5eMrdXRfO7vb097O3tFY7VyMgIPj4+iI2NxfDhwwGILofGxsZKNX1XxdzcHObm5sjJycHx48cldyf17dsXt27dktrW398fHh4emDNnjsoSojLxy7NPXcWvMirvtqeEKVOmMGtra3b69GmWnp4uWV69esUYYywpKYmFh4ezq1evsuTkZPbbb7+x5s2bs549e8p9DOr9Xs7Fi4wBTGBiwh5ev0693wlRkwcPHjAALDw8nN27d48NGDCAdezYkTk4OLDHjx9rOrxK9uzZw4yNjVl0dDS7e/cumzRpEmvQoAF7/vy5ZJt169axPn36SH4/duwY+/PPP9mjR4/YiRMnmJeXF+vSpQsrKSmp8jjV9R4vKCiQygOylrKyMqXirxi7vOesSPyqUJt8pRVJHYDMRXwbQ0pKCuvZsyezsbFhxsbGrEWLFmzWrFkKnTAl9XKEQib08GAMYBnLl1NSJ0SNwsPDma2tLTMxMWF+fn4sKyuLdezYkXl4eGg6NJnWrVvHXFxcmJGREevcuTO7ePGi1PNhYWHM1dVV8ntMTAxr3rw5MzIyYo6OjiwoKIjl5uZWe4zqkmJYWFiVOUG8JCcnKxV/xdjlPWdF4leF2uQrDmOM1XHjgEbk5+fD2toaeXl5Cl+j0EeCZcvA++YbvO7UCYK//8bt27fh7e0NExMTlJWV4enTp3B2doaBQdVXaOTdjhBCiPxqk6+o63M9xcaMAeNyYXrlCjgPH2o6HEIIISpQq6ReWlqK1NRUJCYmIjs7W1UxkbrQuDFed+sGADDcvVvDwRBCCFEFhZN6QUEBNm7ciF69esHKygpNmzZF69atYW9vD1dXV0ycOLHOb9AnyikYORIAYLBrFyAUajgaQgghtaVQUl+1ahWaNm2KqKgo9OvXD4cOHUJ8fDzu37+PCxcuICwsDGVlZRgwYADef/99PHjwQF1xExUo6tcPAktLcFNTYXXtmqbDIYQQUksK9W66cuUKzp49i7Zt28p8vnPnzggICMCmTZsQFRWFc+fOoWXLlioJlKgeMzFB4eDBsN6zB/Z//EHDxhJCiI5TqKa+e/fuKhN6ecbGxpg8eTICAgKUDozUjYK3Qy/anDoFvJ1IhxCi265fvw5DQ0P07t1b7ceqan7xpk2bgsPhVFqCgoJklqPo9srauHEj2rdvLxmtrWvXrvjzzz+r3F4gEGDBggVo1qwZTE1N4ebmhsWLF0NbbxxTSe/3rKws/PHHH/j999+Rnp6uiiJJHSnu0AFCd3fwiovB27tX0+EQQlQgJCQEM2fOxI0bN9R6nOrmF79y5QrS09Mly19//QUAGDVqVJVlKbJ9eYrMue7s7IzvvvsOcXFxuHr1Kvr06YNhw4bhzp07MrdfsWIFNm7ciPXr1+PevXtYsWIFVq5ciXXr1sl1vLpW66T+v//9Dy1atMCiRYsQFhYGNzc3REVFqSI2Uhc4HJS+HSeaR383QnTerl270LBhQwQFBSE3NxePHz9Wy3EKCwsxZswYbNmyBQ0bNqz0vL29PRwdHSXLkSNH4Obmhl69esksT9HtlTVkyBB88MEHaNmyJVq1aoWlS5fCwsICFy9elLn9v//+i2HDhmHw4MFo2rQpPv74YwwYMACXL19WaVyqonBSLywslPp90aJFuHz5Mi5fvozr169j3759+Oabb1QWIFG/0k8/hdDAANzr1wHqMEeIzioqKsLXX3+NFStWwNnZGdbW1pJ5Miqq7dztiswvXlJSgh07diAgIEAyDr4qt1eWQCDAnj17UFRUhK5du8rc5r333kNsbCzu378PALhx4wb++ecfDBo0SG1x1YbCw4D5+Phg5cqVGDZsmKgAAwNkZmaiVatWAESz2xgZGak2SqJednbI6dULtrGxQGQksHatpiMihChh2bJleP/999G6dWsAQJs2bRAfHy+ZrKS82kweo+j84ocOHUJubm6Ns8fJu72yc66L3bp1C127dsWbN29gYWGBgwcPok2bNjK3nTt3LvLz8+Hh4QEejweBQIClS5dizJgxcp1LXVM4qR8/fhxBQUGIjo7Ghg0bsHbtWvj6+kIgEKCsrAxcLlfuaxtEe2QMHSpK6jt3AsuXazocQoiCHj16hJ9//hm3b9+WrPP09Kyypq7s3O3KzC8eGRmJQYMGyT0Pek3b13bOdXd3d8THxyMvLw/79+/H+PHjcebMGZmJfe/evdi5cyd27dqFtm3bIj4+HtOnTwefz8f48ePlOp86peyA87t27WItWrRgERER7NWrV+zmzZvs+vXrkklBtA1N6CKttLSUJSUlsaSkJNGELufPM4GrK2MAK4uKUvuELgKBgOXk5DCBQKDkGRBCyhs6dCgDwHg8nmThcDgyJzBhjLGlS5cyc3PzapcnT55U2u/gwYOVjgOAcTgcxuPxKs2g9vjxY8blctmhQ4fkOg9Ft2dMNMGKeAIwZfTt25dNmjRJ5nPOzs5s/fr1UusWL17M3N3dlT5eTWqTr5SehWP06NEYNGgQZs6cid69e2Pz5s3w9vZWzTcNUve4XAj8/MBdtAjcyEhAzbfCCIVC5ObmwsLCAlwuTUFASG2cOHEC58+fx/Xr16UmV7py5QoCAgKQm5uLBg0aSO2jbPO7ovOLR0VFoVGjRhg8eLBc56Lo9qogFApRXFws87lXr15V+ozi8XgQaukonEol9aNHj+LevXvw8vLCL7/8gjNnzmDMmDEYNGgQwsPDYWpqquo4SR0QfP45DBcvBueff2Dw8CHg7KzpkAghNSgtLcX06dMxa9asShUr8Qxf8fHxle5ZV7b53dLSEp6enlLrzM3NYWtrW2m9UChEVFQUxo8fX2kmx/Xr1+PgwYOIjY2Va/vyCgsLpTpt79mzBwDw/PlzyTp7e/tKXzAAYN68eRg0aBBcXFxQUFCAXbt24fTp0zh+/LjMuIYMGYKlS5fCxcUFbdu2xfXr17Fq1SqtHYdF4SrSV199BX9/f1y5cgVffPEFFi9ejF69euHatWswMTFBhw4dqr2Rn2ixxo2BDz4AAFjGxGg4GEKIPNavX4+XL19KdRITa9KkCczMzKq8rq5uJ0+eREpKiswEmJWVhYcVZoisbvvyfvjhBzg5OVW7pKamytw3MzMT48aNg7u7O/r27YsrV67g+PHj6N+/v8y41q1bh48//hhffvklWrdujZkzZ0pynzZSeD51W1tbnDhxAj4+PsjOzsa7774r6eoPiHocfvHFFzh37pzKg60Nmk9dWllZGZ48eQIAcHBw+G8+9ePHgeHDIbCxAUtJgYG5eZVlpKSkYPXq1ZgxY0aVvUyrOz7NxU4IIZXV6Xzq5ubmSE5OBiDqBVmx92ObNm20LqETBQweDObsDF52NjgHDlS7aWRkJBISEhAZGVlHwRFCCKmOwkl9+fLlGDduHPh8Pnr16qW1TRBESQYGEAYGAgC4mzdXu2lgYCA8PDwQ+HZ7QgghmqVwUh8zZgxSU1Px22+/4fHjx5JBaIj+YAEBYDweOP/8A5S757UiPp+PadOmyX3vKSGEEPVSKKmLhwy0tbVFp06dKt0iUdGzZ8+UDoxoEJ+PV287jWDTJs3GQgghRG4KJfVOnTrhiy++qHZowLy8PGzZsgWenp743//+V+sAiWYUiIdA3L4dqDDef10S38+urfeEEkKINlGo2/Hdu3exdOlS9O/fHyYmJvDx8QGfz4eJiQlycnJw9+5d3LlzBx07dsTKlSvxwdvbo4juefPee2AtWoCTlATs3g1MnKiROGiQGkIIkZ9Cn5K2trZYtWoV0tPTsX79erRs2RJZWVl48OABANH19ri4OFy4cIESuq7jciGcNEn0eONGQLE7HwkhhGiAUjcIm5qa4uOPP8bHH3+s6niIBnC5XJk1YTZuHLBgAXD9OnD5MtCli4YiJIToutTUVHz++efIzMyEgYEBFixYgFGjRmk6LL1D7ZmkyqQOW1tAPDb0Tz/VfWCEEL1hYGCANWvW4O7duzhx4gSmT5+OoqIiTYeldyipk+qJh57cswfIzNRsLIQQneXk5CQZm97R0RF2dnbIzs7WbFB6iJI6qV7nzqKlpATYskXT0RBC9EBcXBwEAgGaNGmi0nKXL1+OTp06wdLSEo0aNcLw4cORmJhY434bNmxA06ZNYWJigi5duuDy5ctSz3/77bfgcDhSi4eHh0pjVxVK6qRmwcGinxs3AqWlmo2FEKLTsrOzMW7cOGyuYcRKZZw5cwZBQUG4ePEi/vrrL5SWlmLAgAHVNvPHxMQgNDQUYWFhuHbtGry8vDBw4EBkVmiZbNu2LdLT0yXLP//8o/L4VUGppP7kyRMcPny40rcZoqdGjQIaNQKePQMOHdJ0NIQQLTRq1CjY29tLJetLly7ByMgIJ06cAAAUFxdj+PDhmDt3Lt577z2Vx3Ds2DH4+fmhbdu28PLyQnR0NFJSUhAXF1flPqtWrcLEiRPh7++PNm3aYNOmTTAzM8PWrVultjMwMICjo6NksbOzU3n8qqBwUt+9ezdatWqFYcOGoWvXrnjnnXfw4sULdcRGtIWxMfDFF6LH69ZpNhZCiFaKiIjAyJEjER4eDkA05/nYsWMxZcoUDBgwAIwx+Pn5oU+fPvj888+rLGfZsmWwsLCodhGPblqTvLw8AKhy3viSkhLExcWhX79+knVcLhf9+vXDhQsXpLZ98OAB+Hw+mjdvjjFjxsgdQ11TOKkvWrQIn332GRISEiTfvubOnavywIjmcLlcNGjQQLo3/OTJgIEBcO4ccOOG5oIjhGglJycnTJ8+Hc+ePcPLly8REhICY2NjrFixAgBw/vx5xMTE4NChQ/D29oa3tzdu3bpVqZzJkycjPj6+2kWe+SaEQiGmT5+Obt26wdPTU+Y2WVlZEAgEcHBwkFrv4OCA58+fS37v0qULoqOjcezYMWzcuBHJycno0aMHCgoKFHmJ6oTC96k/evQIx44dQ9OmTdGqVSvs2LEDPj4+NP2mHhEndSl8PjByJBATI6qt//KLRmIjhGivVq1awczMDAsXLsTOnTtx+fJlyfTc3bt3l2u4Zxsbmypr1ooICgrC7du3VXLte9CgQZLH7du3R5cuXeDq6oq9e/dq3SyVCtfUy8rKYGZmJvndw8MDQqFQ6lsN0X5paWmIiIhARkaG/DuJO8zt3AnQJRdCSAVcLhft2rXDTz/9hCVLlsDLy0vhMlTR/D516lQcOXIEp06dgrOzc5Xb2dnZgcfjVfoczMjIgKOjY5X7NWjQAK1atUJSUpJiJ1cHlOoot23bNvz7778ofDvRh4GBAV69eqXSwIh6RUZGIjExETExMfLv9N57gI8P8OaNqCc8IYSUw94OJ92xY0d89dVXVW4nEAiqfK42ze+MMUydOhUHDx7E33//jWbNmlUbr5GREXx8fBAbGytZJxQKERsbi65du1a5X2FhIR4+fAgnJ6dqy9cIpqCePXsyKysrxuFwGI/HYy1atGAGBgbs66+/Zn///TfLz89XtMg6kZeXxwCwvLw8TYeiFZ48ecJCQkLY+fPn2evXr1lycjIrLS2tecfduxkDGLO3Z6UFBfLvV0Fpaalc+8q7HSFE81atWsUAMG9v70rPDRkyhE2ZMoW98847bMuWLWo5/pQpU5i1tTU7ffo0S09PlyyvXr1ijDG2bt061qdPH6l99uzZw4yNjVl0dDS7e/cumzRpEmvQoAF7/vy5ZJuvvvqKnT59miUnJ7Pz58+zfv36MTs7O5aZmamW86hNvlI4qYvdv3+f7d69m82aNYv17duXNWzYUJLoPTw8lC1WbSipSystLWVJSUksKSmJFRcXs5ycHCYQCOTZkTEXF8YAVrZxIyV1QghjjLGbN28yY2Nj9uWXXzIjI6NK79mmTZuyNWvWqDUGADKXqKgoxhhjYWFhzNXVtdJ+69atYy4uLszIyIh17tyZXbx4Uep5X19f5uTkxIyMjFjjxo2Zr68vS0pKUtt51CZfcRhT3fRbycnJuHr1Kq5fv45ly5apqliVyM/Ph7W1NfLy8mBlZaXpcDSurKwMT548AQC4urrCwECBPpOrVwOhoWAeHnjyxx9wdnFRbP+3x3/69CmcnZ2r3Vfe7QghmvPmzRt06tQJHTt2xJo1a2BjY4Nbt25Jep0XFBSgbdu2WnsbmLapTb5S6YhyzZo1w6hRo7QuoRMVCwwErKzASUiA6alTmo6GEKJhc+fORVFREdavX4+GDRvC1dUVa9asQVpaGgDgzp07ahlshlSmFcPEyjNe75s3bxAUFARbW1tYWFhg5MiRivXcJqpjZSUZjMZKi25tEwqFyM3Nleu2GUKIapw4cQIbNmzAjh07YGlpCQCYP38+Dh06hKCgIADArVu30K5dO02GWW9oRVKXZ7zeGTNm4PDhw9i3bx/OnDmDtLQ0fPTRRxqMup4LCQEzMIDpxYtANUMw1iVK6oTUvQEDBqC0tFSqJj5hwgRkZWXh4MGDACip1yWVXlNXlRcvXqBRo0Y4c+YMevbsiby8PNjb22PXrl34+OOPAQAJCQlo3bo1Lly4gHfffbfGMumaurRaXVN/SzhmDLi7dkH46afg7t6t8PFVfU2drr8TQvSB1lxTV5WK4/XGxcWhtLRUanxeDw8PuLi4VBqfV6y4uBj5+flSC1Et4YwZAADOvn0AdYAhhBCVqpjDiouLa9ynVtWZ2NhYxMbGIjMzs1KTZ8UZbuQla7ze58+fw8jIqNLQpRXH5y1v+fLlWLRokVIxEDl5e+N1t24wPX8eWLsW+PFHTUdECCF6o+J882FhYfj222+r3UfpmvqiRYswYMAAxMbGIisrCzk5OVKLssTj9e7Zs0fpMgBg3rx5yMvLkyypqam1Ko/Ilj9hgujBli3A2xYWQgghtZeamiqVx+bNm1fjPkrX1Ddt2oTo6Ohqp9BTlHi83rNnz0qN1+vo6IiSkhLk5uZK1darG5/X2NgYxsbGKouNyPa6Vy+wNm3AuXtXlNhnztR0SBJCoRD5+fmwsrKSnnGOEEJ0gJWVVd1dUy8pKVHZfYeshvF6fXx8YGhoKDU+b2JiIlJSUqodn5fUAQ4HwunTRY/XrgVKSzUaTnnUG54QUt8ondQnTJiAXbt2qSSIoKAg7NixA7t27YKlpSWeP3+O58+f4/Xr1wAAa2trBAYGIjQ0FKdOnUJcXBz8/f3RtWtXuXq+E/Vin30GODgAT58Ce/dqOhxCCKm3lG5+f/PmDTZv3oyTJ0+iffv2MDQ0lHp+1apVcpe18e2MX71795ZaHxUVBT8/PwDA6tWrweVyMXLkSBQXF2PgwIH46aeflA2fqJKxsWha1vnzgZUrgc8+AzgcTUdFCCH1jtJJ/ebNm/D29gYA3L59W+o5joIf6PLcKm9iYoINGzZgw4YNCpVNqsbhcFR3vfnLL4EVK4CbN4EjR4AhQ2pfJiGEEIUondRP0ZjfOo/L5aJBgwaqSeoNGwJBQcB33wFLlgAffki1dUIIqWPUJZiozowZgKkpcPkycPKkpqMhhJB6p1aDz+Tm5iIyMhL37t0DALRp0waBgYGwtrZWSXBExzRqJJroZc0aUW29f39NR0QIIfWK0jX1q1evws3NDatXr0Z2djays7OxevVquLm54dq1a6qMkeiSmTMBIyPg7FnRQgghpM4ondRnzJiBoUOH4vHjxzhw4AAOHDiA5ORkfPjhh5guvm+Z1D+NGwMBAaLHS5ZoNhZCCKlnalVTnzNnjtRsWAYGBpg9ezauXr2qkuCIjpozBzAwAP76C7h0SdPREEJIvaF0UreyskKKjJm5UlNTYWlpWaugiI5r2hQQDx+8dKlGQyGEkPpE6aTu6+uLwMBAxMTEIDU1FampqdizZw8mTJiA0aNHqzJGoovmzgW4XODwYSA+XmNh0FCxhJD6ROne7z/88AM4HA7GjRuHsrIyAIChoSGmTJmC7777TmUBEh3VqhXg6wvs3g0sW6ax4WMZY8jLy1Pd/fiEEKLFlP6UMzIywtq1a5GTk4P4+HjEx8dLesDT7GgEAPD116Kf+/cDb297JIQQoj61rrqYmZmhXbt2aNeuHczMzFQRE9EXnp7AiBEAY6LaOiGEELVSqPk9NDQUixcvhrm5OUJDQ6vdVpEJXYge++Yb4OBBYNcu0WMPD01HVCWaf50QousUSurXr19H6dv5sq9fv17ldopO6EL0mI8PMGwY8NtvQFgYEBOj6YiqJO5UZ2FhQUmdEKKTFErq5Sdx2bZtG5ydnSt9+DHGkJqaqproiH4IDwd+/13UWW7ePODt7H7qIE7MNjY2ajsGIYRoK6WrI82aNUNWVlal9dnZ2WjWrFmtgiJ6pn17UU94AFi4UK2HEvd2p1vYCCH1kdJJvao50AsLC2FiYqJ0QER3KDR166JFAI8num/94kX1B0cIIfWQwvepizvIcTgcLFy4UKrHu0AgwKVLl+CtxuZVoj3ESV0urVoB48cDW7cC8+cjbetWrF27FjNmzICLi4ta4ySEkPpC4aQu7iDHGMOtW7dgZGQkec7IyAheXl6YOXOm6iIk+mPBAuDXX4HYWEQuWICEzExERkZi0aJFmo6MEEL0gsJJXdxZzt/fH2vXroWVlZXKgyJ6qmlTYNIkYMMGTLp9G/k9eiAwMFBlxaelpSEiIgK+vr5wdXVVWbnyoNvhCCHaQOlPn6ioKEroRHHffAOYmqLxtWuY2749+Hy+yoqOjIxEYmIiYjRw2xyNMU8I0QZKJ/Xly5dj69atldZv3boVK1asqFVQRI85OQFTpwIAGvz4I6DCJBgYGAh3d3f4invaE0JIPaN0Uv/555/hIWN0sLZt22LTpk21CoroudmzwSwtYXz3LjgHD6qsWD6fj5CQEDg4OMi9D9WwCSH6ROmk/vz5czg5OVVab29vj/T09FoFRfScnR3YtGkAAO7ChcDbUQrVTVYCp6ROCNEnSif1Jk2a4Pz585XWnz9/XqXXSYluUDQ5CmfMgMDWFpz794HNm9Uc3dtjUgInhOg5pZP6xIkTMX36dERFReHJkyd48uQJtm7dihkzZmDixImqjJHoAIUTppUVct/W1rFoEZCfL3OztLQ0rF27FmlpaSqKlBBC9JfSSX3WrFkIDAzEl19+iebNm6NZs2YIDg5GSEgI5s6dq8oYiZ4q+PRTsFatgBcvgCo6V0ZGRiIhIQGRkZF1HB0hhOgepZM6h8PBihUr8OLFC1y8eBE3b95EdnY2Fi5cSLO0EfkYGkK4fLno8apVgIyJgAIDA+Hh4aHS+9kVVV0rBDXpE0K0icKDz1SUkpKCly9foqSkBI8ePZKsHzp0aG2LJvUAGzIE6NEDOHdONOJcdLTU83w+H9OmTdNoP43qpmQVP0fzHRBCtIHSSf3Ro0cYMWIEbt26BQ6HI5ngRVxLFwgEqomQ6DcOB/jhB6BLF2D7dmD6dLVOzUoIIfpM6eb3adOmoVmzZsjMzISZmRlu376Ns2fP4p133sHp06dVGCLRe507A59+CjAGzJol+kkIIURhSif1CxcuIDw8HHZ2duByueDxeOjevTuWL1+OkJAQVcZI6oNlywAjI+DkSeD4cU1HQwghOknppC4QCGBpaQkAsLOzk9xy5OrqisTERNVER+qPZs2A4GDR41mzgLIyzcZDCCE6SOmk7unpiRs3bgAAunTpgpUrV+L8+fMIDw9H8+bNVRYgqUe++QZo2BC4fRvYskXT0RBCiM5ROqnPnz9fchtPeHg4kpOT0aNHDxw9ehQREREqC5DUIw0bAosXix7Pnw+8fKnZeAghRMco3ft94MCBksctWrRAQkICsrOz0bBhQ7pPnSjviy+An38Gbt0CFi4E1q7VdESEEKIzlKqpl5aWom/fvnjw4IHUehsbG0ropHYMDCBcswYAwDZtAt5e4iGEEFIzpZK6oaEhbt68qepYCAEACHv2RNHgweAIheDNmEG3uBFCiJyUvqY+duxYGo+bqE32vHlgpqbgnDsHsz/+UKqMtLQ0REREICMjQ8XREUKIdlI6qZeVlWHjxo1455138MUXXyA0NFRqUcTZs2cxZMgQ8Pl8cDgcHDp0SOp5Pz8/cDgcqeX9999XNnSiAwSNG0M4Zw4AwGbZMqCoSOEyIiMjkZiYiJiYGLm2pxnhCCG6Tumkfvv2bXTs2BGWlpa4f/8+rl+/Llni4+MVKquoqAheXl7YsGFDldu8//77SE9Plyy7d+9WNnSiI1hoKFjTpjBITwd35UqF9w8MDIS7uzt8fX3l2p5mhCOE6Dqle7+fOnVKZUEMGjQIgwYNqnYbY2NjODo6quyYRAeYmkL4/ffgjRoFzo8/AgEBgJub3Lvz+XyFRjcMDAxEXl4eRo0aBaFQWGnyFkII0XYKf2o9evRIMnlLXTp9+jQaNWoEd3d3TJkyBS9ruIe5uLgY+fn5UgvRPWzYMLzu3h2c4mLgyy9V3mlOKBSisLAQQqEQfD4fwcHBMDU1palUCSEaVzGHFRcX17iPwkm9ZcuWePHiheR3X19ftXdEev/997F9+3bExsZixYoVOHPmDAYNGlTtTHDLly+HtbW1ZGnSpIlaYyRqwuHgZXg4mLExcOIEsHevyooWCoXIyclBfn4+JXFCiNZp0qSJVB5bvnx5jfsonNQr1tKPHj2KIiU6MSni008/xdChQ9GuXTsMHz4cR44cwZUrV6qdDW7evHnIy8uTLKmpqWqNkahPWbNmkk5zmD4dyMtTSblCoRB5eXkaaXkihJCapKamSuWxefPm1biPTl40bN68Oezs7JCUlFTlNsbGxrCyspJaiO5is2cDrVoBz5+LxognhBA9VzGHGRsb17iPwkldfEtZxXV16enTp3j58iWcnJzq9LhEg4yNgY0bRY9/+gm4ckWu3TIyMhAREUG3qRFC6gWFe78zxuDn5yf5xvDmzRtMnjwZ5ubmUtsdOHBA7jILCwulat3JycmIj4+HjY0NbGxssGjRIowcORKOjo54+PAhZs+ejRYtWkiNP0/qgT59gLFjgR07RGPEX74MGFT/LxwTE4PExERERkZiwYIFKgtFKBQiPz8fZmZmKiuTEEJqS+Ga+vjx49GoUSPJhfuxY8eCz+dLXcy3trZWqMyrV6+iQ4cO6NChAwAgNDQUHTp0wMKFC8Hj8XDz5k0MHToUrVq1QmBgIHx8fHDu3Dm5miKInvnhB6BBA+D6daCacQ3EfH194e7ujsDAQJWGIRQKkZubK+lgRwPXEEK0gcI19aioKJUH0bt372o7Kx0/flzlxyQ6ysEBWLFCVFOfPx8YPhxwda1mcweEhISAz+dLEu+MGTPg4uKi0rCioqIkA9csWrRIpWUTQoi8dLKjHKk9LpeLBg0a6OYAKxMmAN26AYWFouT+9gshl8uFtbV1lX081DlinL+/Pzw8PFTeIkAIIYrQwU90ogo6ndS5XCAyUtR57vhxYPv2t6urP6fAwEC1JV4+n49p06aBz+crtF/FZnxCCKkNHfxEJwSAuzvw7beix9OnA+npNe6ibOJVJ0rqhBBVoqROdNfMmYCPD5CbCwQF0bzrhJB6j5I60V0GBqJmeAMD4OBBYP9+hYvIyMjA1q1bq+21Tj3bCSG6olZJPTY2Fl9//TUmTJiAgIAAqYWQOuHlBYiHTpw6FWl37iiUgPft24eHDx9We1cHTclKCNEVSif1RYsWYcCAAYiNjUVWVhZycnKkFkLqzDffAG3aAJmZ2DpmjEIJeNSoUXBzc4O/v3+V26irgx1dTyeEqJrS86lv2rQJ0dHR+Pzzz1UZDyGKMzYGoqOBrl0x6c4dvHJxkTsBOzg4ICAgoNrOc+rqYCeeUIYQQlRF6Zp6SUkJ3nvvPVXGQojyOnUCvv4afADLL1wAv47nI1AHqskTQhSldFKfMGECdu3apcpYiA7Tis5k8+eDdegAXm4uuOUGpdFVlNQJIYpSuvn9zZs32Lx5M06ePIn27dvD0NBQ6vlVq1bVOjiiO8p3JtPYMKlGRhBERYHXuTO4f/4JFhkJ9O+v1kOmpaVh3bp1CA4O1s2BfAghekXpT6GbN2/C29sbXC4Xt2/fxvXr1yVLfHy8CkMkukCdo7UppG1b5MycCQDgzpoFg5QUtR5OPOa7OuZEIIQQRSldUz916pQq4yA6TptGa8sPCEDDs2fB+ecfNJo1C2k7d6rtWP7+/igsLKy29zwhhNQVai8k+ofHg2DrVjBLS5jGxaGhHFO0KkMoFMLMzAzBwcFa8WWGEEJqldRzc3Px448/YsKECZgwYQJWrVpFt+gQ7dCsGYRvk3nDDRuAc+fk2i0jIwMRERFydfhjjCE/P7/aaYMJIaQuKZ3Ur169Cjc3N6xevRrZ2dnIzs7G6tWr4ebmhmvXrqkyRlLPqKonPRs9GvkffQSOUAjeuHFAdnaN+8TExCAxMZFGjyOE6CSlk/qMGTMwdOhQPH78GAcOHMCBAweQnJyMDz/8ENOnT1dhiKS+UeWwrFlhYShp2hScp0/lus3N19cX7u7uddbhTygUorCwkG5bI4SoRK1q6nPmzIGBwX997QwMDDB79mxcvXpVJcGR+kmVPemZuTky1qwBMzQE99AhWFYYW6HihC4ODg4ICQlR2TXysrIypKSkoKysTHZ8jFWb1OledUKIIpRO6lZWVkiRcbtQamoqLC0taxUUqd9U3ZO+xNMTwmXLAAANFy8Gbt+WPFd+Qhculwtra2tw5BiNLi0tDREREcjIyKh2u7KyMqSlpVWZ1GvCGENeXh4ldUKIXJRO6r6+vggMDERMTAxSU1ORmpqKPXv2YMKECRg9erQqYySk1lhICITvvw9ucTF4Y8cCr18DkJ7QhcvlokGDBnINIhMZGYnExETs27dP3aETQojclL5P/YcffgCHw8G4ceMktRBDQ0NMmTIF3333ncoCJEQluFwIIyPBvLzAu3MHmDEDWL9ergldZAkMDERubi4+/vhjNQVMCCGKU7qmbmRkhLVr1yInJwfx8fGIj4+X9IA3NjZWZYyEqEajRnixahUYhwP8/DM4O3YoXRSfz0dISAgcHBxqFVLFa/qEEFIbtR58xszMDO3atUO7du1gZmamipgIUZs3PXqAffMNAID75ZcwSkjQaDzlr+kTQkhtKdT8HhoaisWLF8Pc3ByhoaHVbksTuhBtJZw/H9wrV8A5fhz8kBA8//lnjcUyatQoFBUV0TCzhBCVUCipX79+HaWlpZLHVZGn9zCp3xTplKZyPB6wcydYx44wevIELZcuBXr2rPs4gGqv6Yt72Pv6+sLV1VUD0RFCdI1CSb38JC7btm2Ds7NzpQ9lxhhSU1NVEx3RW+KkrjG2tkhdtw5bR4zApHPnYPXjj8D8+VVuLr5f3MbGptL/vFAolOu2M6FQiPz8fFhZWSnUwz4mJgadO3eW77wIIfWa0tWkZs2aISsrq9L67OxsNGvWrFZBEVIXIuPicNXDA1sAGISFAdXMPFjd/eLi52oaA17RgWQCAwPh7u4OX19fubYnhBClb2mr6gOssLAQJiYmSgdEiLpUbPIPDAxETnY2hjs6gvP338CnnwKXLqk1hvI1/pqIe9gTQoi8FE7q4g5yHA4HCxculOrxLhAIcOnSJXh7e6ssQEJUpWKTP5/Px9TgYLxMTYUwKwvcmzfB+/jjWt3qVhNxrb6qSw/lm+gJIURRCid1cQc5xhhu3boFIyMjyXNGRkbw8vLCzJkzVRchIWomNDFBaUwMjLt3BycuDrZz5gD792smlrc1eQsLC40cnxCi2xRO6uLOcv7+/li7di3VKIha1VUveda0KfC//4H16weL33+HYOVK4O397PLgcDiwtrZWKM60tDSsX78effv2VSJiQgipTOlr6uLBMu7evYuUlBSUlJRIPT906NDaRUYI5Oslr2yTdaVpT3v1gnDtWvCCgsCbPx/w9ASGDVMoTkWSurh3e1FREfr374/c3FwawIkQUitKJ/Xk5GQMHz4ct27dAofDkXScE9+jLhAIVBMhITVQtslaKBSiqKhIqtMnmzQJ+RcvwurXX4ExY4ALF4DWrVUdMgBRR73U1FS8evUKz549g6GhIXUyJYTUitJtmiEhIWjWrBkyMzNhZmaGO3fu4OzZs3jnnXdw+vRpFYZIiHIUmUq1vOwFCyDs0wcoKgKGDgX35UulyqkJn8+HhYUF0tPTER0drdKyCSH1k9JJ/cKFCwgPD4ednR24XC64XC66d++O5cuX0204RCsofT3e0BDC3buBFi2Ax4/B/fhjNDAxUeq6flpaGiIjI6ucsEU89aufn5/CZRNCSEVKJ3WBQABLS0sAgJ2dneRDy9XVFYmJiaqJjpBaUjqx29gAv/8OWFsD//wDbkAAIBQqPIBMVFQUkpKSqpywRTxMrJOTk2LxEUKIDEpfU/f09MSNGzfQrFkzdOnSBStXroSRkRE2b96M5s2bqzJGQpRWq+FoW7cGDh4EBg4Ed98+NGzYEMK1a5Gbmyv3tW9/f388f/5c7glbxJ33ysrKJJ34NDI+PiFEJyn9aTF//nxJbSU8PBzJycno0aMHjh49ioiICJUFSIhG/d//AVu3AgCsN28Gb9MmhXbn8/kIDAyUOWGLLOWTujxDzxJCSHlKJ/WBAwfio48+AgC0aNECCQkJyMrKQmZmJvr06aNQWWfPnsWQIUPA5/PB4XBw6NAhqecZY1i4cCGcnJxgamqKfv364cGDB8qGTogEh8OBhYVF9bXhsWMhCA8HABh89RVM//qrjqIDMjIyEBERUeU1eUIIKU+l7Xo2NjZK9RAuKiqCl5cXNmzYIPP5lStXIiIiAps2bcKlS5dgbm6OgQMH4s2bN7UNmdRzciV1AGzuXBT4+oIjFMI+JAScK1fqJL59+/YhMTERkZGRtS5L0f4AhBDdo3RS79OnDxYtWlRpfU5OjsI19UGDBmHJkiUYMWJEpecYY1izZg3mz5+PYcOGoX379ti+fTvS0tIq1egJqa0qO9ZxOHi5eDEE/fuD++YNjD76CNlXrmD9+vXIyMiQbKbqxDlq1Ci4u7sjMDCw1mVRUidE/ymd1E+fPo3169dj+PDhKCoqkqwvKSnBmTNnVBIcIBrk5vnz5+jXr59knbW1Nbp06YILFy5UuV9xcTHy8/OlFkJqUm1veUNDlO7cieJ27cDJysKRSZNw/8YNREVFSZrIVZU4xdfW7e3tERISIvc1eUKI/qiYw4qLi2vcp1bN7ydPnsTz58/x7rvv4vHjx7UpqkrPnz8HILr1pzwHBwfJc7IsX74c1tbWkqVJkyZqiY/UM5aWyIiKAnN3xxcFBeh86xY4b95U2URe033qVRGPdmdhYaHyQW8IIbqhSZMmUnls+fLlNe5Tq6Tu5OSEM2fOoF27dujUqZNWjSQ3b9485OXlSZbU1FRNh0T0hNDWFoKjR+Hk7IzwnBx8c/06PJo1k9lEXtN96tVd06+ryWwIIdopNTVVKo/Nmzevxn2Uvk9dXHswNjbGrl27sGTJErz//vuYM2eOskXK5OjoCEDUC7j8AB0ZGRnVzttubGwMY2NjlcZC9A+Hw4GVlZXiidPFBYKjR8Hp1QvNExKwzMYGRra2lTbz9/dHeno6Ro4cKbNJnsvlytVRjxBS/1hZWSk8UZXSnyQV75+dP38+du7ciR9//FHZImVq1qwZHB0dERsbK1mXn5+PS5cuoWvXrio9Fql/alUbbt0a6ZGREJqZwezff8EdPRqoMFshn8+Hv78/zMzMqIMaIUTtlE7qycnJsLOzk1o3cuRIXLx4EVvfDtYhr8LCQsTHxyM+Pl5Sdnx8PFJSUsDhcDB9+nQsWbIEv//+O27duoVx48aBz+dj+PDhyoZPtJguNTuXeHkhd9s2CI2NwT18GNyxY4HSUqXKqk0nO/G+ZWVl1MOdkHpM6eZ3V1dXmes9PT3h6empUFlXr17F//3f/0l+Dw0NBQCMHz8e0dHRmD17NoqKijBp0iTk5uaie/fuOHbsGE1TqUVUmYhrNbRrHeNyuTD78ENk/vwzHCZNAvfgQdiXlAD79wMG8r+90tLSsHbtWnz00UcwUGA/MXFSNzExkUxDqwtfigghqqXQp0doaCgWL14Mc3NzSeKtyqpVq+Qut3fv3tUOh8nhcBAeHo7wt6N6Ee2jLYlYU7X8N716Qbh/P7gjR8L8jz8gDAgAfv1Vapu0tDRERETA19e30pfiqKgoJCYmIjo6GlwuFyNGjMChQ4cwZ84cGBoa1uWpEEJ0mEJJ/fr16yh927R4/fr1KrejW3CIpmjyywUbNAjCPXvA9fUFd/duwNgYKDdKYlRUFG7cuIFLly5h7969Uvee+/v7Iy8vD9nZ2UhNTUV4eDjMzMwQHR2NiRMnqixGcY3exsaGavKE6CGF3tWnTp1CgwYNUFpaCi6Xi82bN+PUqVOVlr///ltd8RKiMsreQ14dNnQoXkREgPF4QHQ0DKZOBd5e3/b398ebN29gYmIiuaedw+HA2toazs7OmDp1Knr06IH8/HwEBASoZZ51xhjy8vLomjshekqpr+qGhoa4efOmqmMhpE7VdA95eeJr3vJ8AXg1aBCE27cDXC4MoqLQ4rvvAIEAfD4fq1atgpeXl+Se9oqXCy5cuAArKyvcunVL6XnW1T0cLA03S4j2Urr9bezYsSqZZIIQTfH390eLFi3kmus8MjISCQkJcn0BAAD2ySfAr7+C8XhwOHoUjjNnAqWlcHBwqHbY18GDB6N58+YYOXKkQucidewqauPia/rlx6pXBiV1QrSX0r3fy8rKsHXrVpw8eRI+Pj4wNzeXel6RjnKE1Ia4Fj1jxgy4uLjIvZ8ic50HBgYiPz9fri8AEp99hjIeD7yxY2F19CgEn30GzooVYNUMimRnZwc/Pz/Y2toiLy9P/mPJITIyEomJiYiJiUHnzp0V2lcoFCI/P1/hgTAIIXVL6aR++/ZtdOzYEQBw//59qeeooxypS+JadGRkpMyZA1WBz+dj2rRpcHR0rHbOgYo4I0cio7AQDkFB4B05AsecHDzftEllcYmTrZmZWY3bBgYGIjc3F76+vkodR3yrHCFEeymd1E+dOqXKOAhRmrgWrYrpSWtLVqvBq9698fyXX+A0eTLMzp+HU0AAcOwYoIJab/n702vC5/MREhJS62MSQrQX3dNCdJ64Fq0N05OWbzUo73XXrig5fBgCCwuYXrkC3oABQGamzDLE065WN3YDIYTIonRNXezu3btISUlBSYUxr4cOHVrbognROdW1GrCuXZG2Ywf4/v7gxcWB16sXbPbvB9fFpVKns1evXtGgM4QQhSmd1B89eoQRI0bg1q1b4HA4klqF+Hq6QCBQTYSE6JCaWg1KPD3xLCYGTSZOBCcpCVbvvw/8+SeEFYZWzsrKwrFjx+Do6ChXck9LS8O6desQHByskvMghOgmpZvfp02bhmbNmiEzMxNmZma4c+cOzp49i3feeUer5lUnpCrKTHsq6351RYemLW3eHIKzZ4H27YHnz4GePcGp0Efl6NGjePjwIaKjo+UqMyoqqtItd3TrGSH1j9JJ/cKFCwgPD4ednR24XC64XC66d++O5cuXU2ccohOUSeri5Fn+mrlS483z+cDZs0Dv3kBBAbgffgiLo0clT3/wwQcKjSjn7+8PDw8PqVvuZCX1jIwMREREqHQUPUKI9lA6qQsEAlhaWgIQ3Vsr/pBwdXVFYmKiaqIjRMuIk6dKetpbWwN//gl8/DE4JSXAjBk4NHMm7icm4ujRoxgyZAgaN24s1xcGeTsLxsTEIDExkQaOIkRPKZ3UPT09cePGDQBAly5dsHLlSpw/fx7h4eFo3ry5ygIkRJuovKe9iQmwZw+EwcHYAuDptWvY9vXXeJSUhOPHj8PAwECls875+vrC3d29Vl9KFBkylxBSt5T+pJg/f76kWS88PBzJycno0aMHjh49ioiICJUFSIimqf3aNI8H4apVGD5jBtpxOFifk4NuSUn4/MMPVT6TWk3D1Mqjqtv2CCGap3Tv94EDB0oet2jRAgkJCcjOzkbDhg1pRDmiV4RCoUqHbJU1/SmXy4XjnDkY4eiItuHh6JiRgZKpU4E//gDatlXZsVVBmwb7IYRIU2k1wMbGhhI6ITWQNeEKl8uFtbU18t59F3EREXjt6AijJ0/A694dUHAqY3V0hivf5K5Ng/0QQqQpVFMPDQ2Ve1ua0IUQkYyMDMTExGDatGnVbpeWlobo6GgMGjQIpRs2wGfJEphevw4MGAD88AMwbRogx5fmPXv24O7du9i6dSsmTJigknMo3+S+YMEClZRJCFE9hZL69evX5dqOauuE/Efc4zwqKqraW9SioqLw8OFDnDhxAmPGjMHTbdvg9t134O7aBcyYAVy5AmzZAlQzeUtaWhoKCwthb2+P8ePHSz0nbg3Iz89X+ByoyZ0Q3aBQUqdJXAj5j7z3p/v6+iImJgb+/v7Vjufu7++PtLQ0jB49GlwuF8zYGMLoaHA7dwa++grYtQu4cwc4eBBo1kxmGVFRUXj27BmcnZ3h5OSErKysSvEWFhYqfK7U5E6IbqjV2O+5ubmIjIzEvXv3AABt27ZFQEAArK2tVRIcIdpMnCRrUr7H+bNnz6rcjs/nS94/2dnZopUcjqjZ3dsb+OQT4MYNwMcH2LNH1CxfgZ+fH168eCHVkbUq5edIV3Uve0KIZij9Tr569Src3NywevVqZGdnIzs7G6tWrYKbmxuuXbumyhgJ0VniJm95L0lxOBxYWVlV3r5XLyAuDujcGcjJAd5/H1i8GGmpqVL3jDs4OGDMmDGwtbWt8VhV3apHw8sSoruUTuozZszA0KFD8fjxYxw4cAAHDhxAcnIyPvzwQ0yfPl2FIRKiuxQdQpbL5VZdc3Z2Bs6cASZMABgDFi7E1r59kXjjhtSY77VFSZ0Q3VWrmvqcOXNgYPBfC76BgQFmz56Nq1evqiQ4QtRJqTHbNc3ERNRZLjoaMDPDpORkdLp4ERPc3dV2LoomefpSQIjmKP0JYGVlhZSUlErrU1NTJWPCE6LN5E3qaWlpiIiIQEZGRq3KlLcpPj09HVu3bq3+eOPHA1evwsnTE4tev4bLuHGw/v57cGtIpMp8kREPvlNdJ7/y24ovx1FSJ6TuKZ3UfX19ERgYiJiYGKSmpiI1NRV79uzBhAkTMHr0aFXGSIhGRUZGIjExETExMTVuW1NSt7GxqXHUxe3bt+Phw4fYt29f9Qdr3RqCf/9F3qefgsMYOMuWwWjgQBhlZgIQfTlQdJrY2tSyxQk9JydHri8AhBDVU7r3+w8//AAOh4Nx48ahrKwMAGBoaIgpU6bgu+++U1mAhGhaYGAgcnNz4evrW+uy5LmtzM/PD7m5ufjkk09qLtDUFFlLluDNu++i0YIF4P77LzreuoWk0FBEl5ZKBoxZtGiRXPHVZkhcRWr0hBD1ULqmbmRkhLVr1yInJwfx8fGIj49HdnY2Vq9eDWNjY1XGSIhG8fl8hISEwMHBoU6O17hxY8ycORNOTk5ybS8UCvG8d2+UXrwIYceOMCwoQOtFizDp1i20btpU5QPG6GRfBELqCYXflUKhECtWrEC3bt3QqVMnhIeHo0WLFmjXrh3MqhnpihB1qe9JhjGGwsJCCJs3R8np07j8ySf4lsOB4eHDWHnoEBq/HUdCVer7602INlP4Xbl06VJ8/fXXsLCwQOPGjbF27VoEBQWpIzZC5KIvSYbL5cLCwgJcLhccDgfW1taKn5OhITaZm+Oclxc2N2wIg+fPwXv/fdEANq9fqydwGagHPCGaofCn4Pbt2/HTTz/h+PHjOHToEA4fPoydO3fSm5cQOZWf8ay88km9Nl9UhgwZAtcOHTDm1Cnkjx0rWhkRAXh5AWfPquIUqpSRkYHvvvsOy5cvx/379+lzgZA6pvAnRkpKCj744APJ7/369QOHw1HpNI+EaBula84ylJ/xTB3s7e0xdepUNG7ZEtmLF0Nw5AgYnw88eCAame7LLwEZk7pUVbtOS0vDunXr8Pr16xrPPyYmBmfPnsXZs2cldwtQrZ2QuqPwJ1RZWRlMTEyk1hkaGqK0tFRlQRGibVTZxB8YGAgPDw+1zXjG4XAkNX4AYAMHQnDjBgo+/VS0wcaNQNu2SP/1V6kWg6p6r4tv6du3b1+15y8UCvHhhx+iR48e6Nmzp+RuAXUldfqyQEhlCt/SxhiDn5+fVA/3N2/eYPLkyTA3N5esO3DggGoiJETPqHPGs7S0NGzduhV+fn5wcXH574kGDfBy+XKYBQaCN3ky8PAhtgYE4JGTEyJ5PCz64Ycqy5R32lWhUAgzMzPMnj0bPB4PPB5PrpgrTiwj70Qz4qRe/gsMIfWdwu+E8ePHo1GjRrC2tpYsY8eOBZ/Pl1pHCKl74jnZqxq4hvXuDdy8CcyciQkcDjqmp2PS1q3AL78Ab2u8GRkZiIiIkNTg1T3tasUaN9XACVGewjV1VU4cQYguKN/0ru2JRjwn+6hRo6reyMwM+P57OHz0EeYFBMAoIQGYOBG8LVtgOG8edh47huTkZKlBa8SJ1sbGhmrFhGgxencSUoPySV2dt8+pomzxnOxyDZTTqRPSfv8dgu+/BywswLl8GU1GjsTEJ0/gUWHQGsYY8vLytP5LDSH1HSV1QhSg7UldYYaGYNOnA4mJEPr6giMUwufYMfz4v/+h8V9/AUIhhEKhaHAbSuiEaD2dSOrffvstOByO1OLh4aHpsAjRGHFzOAC8fv0aUVFRCt9WKpWs+XwId+xAanQ0Xrm4wODlS/AmTAA6dQI7exaFhYU1jumelpaGqKgoFBUVSU1YU9V9+RWff/r0KXJzc5GSkoLvv/8eT58+Veh8CCE6ktQBoG3btkhPT5cs//zzj6ZDIkRjyif1/fv34+HDhwr3d5FVA3/cvDmm9u6NW1OmgFlaAteuwbh/f7h/8w0MnzyptjxxJ71Dhw5JzUQn67788p3hyj+fm5uLyMhI3L17V2338ROiz3QmqRsYGMDR0VGy2NnZaTokQrTCBx98gPz8fHzwwQfgcrmwtLRUqBlfnGCfPn2KmTNn4m5iIrYYGECQkABMngzG5cLuzBk0HTwY3NmzgbdfJiry8/ODq6sr/P39pY4fGBgId3d3jBo1SmYPd/F9+/7+/pJyWrRoIfm9KjW1ABBSH+lMUn/w4AH4fD6aN2+OMWPGICUlpdrti4uLkZ+fL7UQognqvlZ+9OhRWFlZ4ejRo+ByubCyslKol7o4wUZGRsLExATFxcWi3vONGgEbN6Ls6lUUdO0KTmkpuKtXAy1aAD/+WGkseQcHB3zyySewt7dHQkIC1qxZg7S0NPD5fAQHB8PU1LTSdXnxve3BwcGSW+acnJwQGBhY4y10kZGRuHfvHjZs2EDX+4leqpjDiouLa9xHJ5J6ly5dEB0djWPHjmHjxo1ITk5Gjx49UFBQUOU+y5cvl7pvvkmTJnUYMSH/UXdS9/f3l6tmK085Xl5eWLBggVTveU67dsj89VekR0aCtW4NvHwJzJwJtGwJbN4MvB1NMj09Hbt27cKzZ8+wZcsWuYbCrc0c7OIWgGHDhlFSJ3qpSZMmUnls+fLlNe6jE0l90KBBGDVqFNq3b4+BAwfi6NGjyM3Nxd69e6vcZ968ecjLy5MsqampdRgxIcpT9EsAn8+XWbMtX448TdV8Ph9Tp06VeWkrIyMD3924gdQjR4DISKBJE+DZM+CLL4A2bYDduxEdFYXk5GRs374dvr6+cHd3V2oo3PT0dERGRsqMtXyzfV3Pc09IXUtNTZXKY/PmzatxH51I6hU1aNAArVq1QlJSUpXbGBsbw8rKSmohRBfUlNTlHXGtfDm1nUQmJiYGiYmJiNy2DQgIEE0Os2YNYG8PJCUBn32GyQcOwNPAAOM+/xwODg4ICQlRaBQ68ReP9evX4/Hjx4iOjq60Tflzp1vtiL6rmMPKD89eFZ1M6oWFhXj48CGcnJw0HQohdU6ZJuvqJpGRJzlWqnkbG4vmaH/0CFi8GLCyQpPERPx89ixaffIJzI8elQw7WxPxDHjR0dFISEgAh8OBp6cnJkyYIBVjxS8y4riVabonRF/pRFKfOXMmzpw5g8ePH+Pff//FiBEjwOPxMHr0aE2HRojCNDHIjKzx28XztwOoMalXWfO2sADmzweSk1Hy1VcoNTUF784dOIaEoMmgQeDs3AmUlUntUjFBi1+PCRMmwMPDAyEhIZVipfHgCZGPTiT1p0+fYvTo0XB3d8cnn3wCW1tbXLx4Efb29poOjRCFacuodOKkXn7bjIwMREdHIyMjQ7ED29ig5NtvcToqCm/mzIHAygpGDx+C5+cHXrt2aHTkCFBSgtzcXJSVlclsaVBm4pisrCysX79e6vo7fQEg9ZlOJPU9e/YgLS0NxcXFePr0Kfbs2QM3NzdNh0WI1qntF4Z9+/ZVmuUtLS0NERERNSb69PR0bDt8GJcGDcJXn32G2198AWZnB05SEiynT4fr//0fOCtXQpidrVRssvz++++ia/1VDGxDSH2jE0mdECKt4vSoyhB3TEtPT5esGzVqFNzc3BAQECD5YhAZGYnExETExMRUW150dDSSk5Mxf/58JCQnY4uxMQRJScAPP+BZo0ZYmpGB4u+/h3HLlij5+musX7Kkyh7u+fn5yM7OrjExDx06VOle9oToI0rqhOggSW/0WgylKu4Rv27dOkRGRiI9PR0ODg6STnUVR4Tz9fWtVEb5WrGfnx+aNWuGJUuW/Le9uTnw1VfYPHEiLrdti59sbMApLMTevXuREh2NbUOGwDA+vlLNWiAQ4NmzZyircD2+/HHz8/Ph4uIiNXAN1dJJfUdJnRAdVJv7wMXEyRsAkpKSsG3bNlhbW0tNxgKg2vvByydRJycnfPbZZ/Dy8qq0feCkSXDr0wedIiPx+tAhjH3nHbQH8MXdu2g8YgQs+vUDfv0VnOJicLlcmJub49WrV5Jb18THKD/5S0FBAczNzaUuNdRmMBtC9AEldUJ0kDL3gVck7pgWHByMFi1aICAgANbW1gpdj5f3XnHJwDb29hD07Qvs2YOAw4chGD4c33K5yL5xA0YTJqBJt24wXLAABuVmaCufqKOiopCQkICoqChkZGRg/fr1CnXqo5o80XeU1Amp5+Qda12sqgFguFwuTE1Nq/xSkJGRga1btyI+Ph4RERFItbHBZk9PXOrZE2u7dgVr0gTcly/B++EHNOvXD63nzAH37ZzuYv7+/pLJX8SXICpe66+uvwEldaLvKKkTomMq9kY3MDAAn8+HgYGB2o4pHiCGy+VWmRjLJ/Xy24uJe9Z//fXXkmTs7++PVu3bY9gPP0Bw/z5w8CCEffsinTGsO38eWUOHgufujjfLl2P90qUAILntTXwJouK1/or9DSiRk/qEkjohOqZib3QDAwO4uLgondTLj7Uu6951QL5b5crvK2t7cc/6ZcuWSZKx1PV6AwNg+HA83boVH7Vvj7OOjthsbAzO48fYHxmJlKgobOvXD2ZHjgDFxXBwcEBwcDBatWoldZyK/Q0oqZP6hJI6ITqmut7oyoiOjkZSUhKioqKqTOpVSUtLq/ELgZiDgwMCAgLg7e0t1ZGuYnN5ZGQkjBs0QI6TE8bFxUGwbRvGdugg6liXnIxGwcF42bQpfhk5Ern//IMGFVoEqutvUFdzsNMXCaIplNQJ0TGqnp3Mz89P6albIyIicPLkSUREREjWVVWrF19TL39fPFC5uTwwMBBeXl4ICwsD380NqT174vvOndHpp59gN3s2ypyc8EtuLu7fuIGjAQHgtW8vGn/+4cMa4y0/sY06Eq+4zLKyMkrqRCMoqRNSD1TXfK5oR7majlPVscTX1CvOvlaxubziFLCRkZG4d+8e9l+5gpIFC/D03DmM37oVHi1aINDICJyEBGDhQqBFC2R07oyfAwORlZBQKTahUIhRo0ZJjqXOpE7JnGiK+nrWEEK0hjjRqlpISAiKiooQEhICANUms1GjRqGoqAh+fn5S62U1l4vHoHdycoK/vz+ePXuGDz/8UPQkjwfH0aMR0LMnSgoKIIiLAy8mBoiNxS/Xr+MhgD/PnkWn/v2BTz5Bmo8P1m7bhuDgYJiamkoNViOLeGAbCwsLFBYWwsrKqtrLEeLtaXpnog2opk4IUUj569J8Pl/uWr74mro8UyaLa/VRUVHg8/nw9/eHnZ1dpVYAZmkJNn48cOIE8PQpAhYuRBsHB0xiDNwTJ4AJE7C1Y0ck//orooOCwHvxosZjK9qELqt2XlfX7gmpiJI6IfVcTb3WKyp/XVqRYzRo0ACWlpYAgKKiIlhaWlZ5HHFPefF1/qysLERFReH+/ftV15ydnOC0YAECL15Eo7t3gfBwwMsLkxhDh5cvMfnYMbi+9x74n34Kztq1QEqK3PGL1dS8Lk7mERERCr9GtT02IQAldUJ0kqz7wJVVU1KvuE48vKwiQ9RyuVzY2NhImqhfvXpVbbO2o6MjQkJC4OzsDAA4fPgwHj58iK1bt0qSWrWvQcuWwIIFQHw8GiUkIHjOHAjbt0c4Y8i5ehW8mTNh4OYG/gcfgLtwIXDxIiAQSBUhq7ZdU61c/IUHgMKvUU0oqRN50DV1QnSQKuZkl7eMitfjy897LmvCFXHyqSpp13TrG5fLRcOGDcHj8SSD3QwZMgS///47Ro0apXD8cHND/uTJiDQ0xNWGDVHC42FJWRnSzp7FL/fuYeK9e+AvXw7Y24MNGoSH7dtjY1ISYGws6ZW/aNGiKosv33IRGBiI/Px8BAcHg8vlqqTzISGKoKROSD0lTopVzYQmbxniBC2usYoTWlWJu+JzHA5H6guAOK7CwkIAoprwkSNHMGrUKDg6Okpto8iXGn9/fxQWFuKLGTPAcXHB5lmzEHf8OJhQiG+fPgVevABn+3bsAZDK4cC1cWN4uboisEsX0VC1VRxLnMjFfQuCg4NhZmaGV69eKfJS1kj8+s6YMQPOzs6SznmqaK0h+oP+GwghSjMwMICzszMMDAwkNdaoqCiZ21bVtG9jYwMbG5sqk1NUVBQePnyIffv2Yf369ZKBbuRN6uJtnZ2dJS0MABAYHIyWffsi4OhR4MUL4NQpCGfMwHgXF7RnDDOfPsV358/DZfBgwNERGD0a6atXI6LCPPCOjo7w9/eHiYkJhEIhGGMoKChQecJVxz321KSvfyipE1LP1aYpv/y+4mvtVQ1iI8/1elkCAwPRunVrAKg0h7w8San85YPyM8qVv4wAQ0Ogd28IV66E4O+/4XfyJOx/+AGv+vQBs7AQJf09e7B19mw8iYzE1o4dgeBg4LffIMzKQkFBAXJzcwEA1tbW4PF4tb48Iut1EF+nr03v+vKv2dOnT7Fo0SI8LTcrHtFtlNQJ0TGquJ6ujvKkkmQ1FE1Izs7OmDlzJgICAirNIa9ITbOsrAwvXryo9nKDeLKcZ6amEAQFITMyEoKMDODsWWDBAkzo0AHtOBxMevkSWL8eGD4cJo0bw3v8eNgvXgzewYN49eQJ1q1bJzk/VdWGy7++ytyBIFY+kdemHKKdKKkTomNUndTrmrKJRJ455Kt7bYRCIV6/fi2ZMlZWoi0/WY6kLBMToEcPIDwcDpcvIyA+Hg779gFTpgAeHkgDsPDhQyzdsQMZn36KqHfeweMtW7B1xAhg1y4InzxRSVKX1SqiTO/6ih37VN1Ln2gWdZQjhKhUWloa1q1bhxkzZsDFxaXS8+U7llWlfAITCoWSznQ1dT6rbuS8jIwM7N27F+7u7nB1dUVeXp7M2PLy8jBhwgQYGBjILEtoZQU2fDjw8ccAgJ+/+gqHd+4E580bWBgbY2JmJn4pKsLEa9eAMWNgAMDZ0RHc994D3nsPePddoGNHwNS02nOp7tzkbRWRpfzr7+zsjLCwMBoNT49wGGNM00HUhfz8fFhbWyMvL4/+gQmphfLDopavEZeVlUmadK9evYrOnTvLvBVMvJ24g11NxNs7Ojri+fPnUvspUtaCBQvwzz//oHv37ggLC8OTJ08AAK6urnKXJ+vcU1JS8O2336KsrAzh4eFwNjPDy99+g/29e+CeO4dn167hF6EQEwFIUrChIeDtLUrw4qVZM4DDkeu1rmq9PGqzL6kbtclXVFMnhCikpnHkxbeO1UWTbvnbvGS1ClSM6/nz55KOfBVvpZOH+NzL34vP5/Mxd+5cvHjxQnTLnYEBXvfvD6G/P7gGBtj89de4GhsLoaUlFllaAhcuABkZwJUromXdOlHhNjaiGryPj+Sn0MUFubm5lW4PrM1Y/uqaB4BoB/qaRghRiapuHVN1+eWTm7LX58WXCF69elXjZC2yrodX1/mt0gh8kyej1XvvIXDrVuDgQSA9HXj8GNi9G5g2DejSRVRzz84GTp4EVqwAfH1Fs87Z2uKXHj3w/MsvRdsnJlYa+Y6Q8qimTghRCVUMZiNP+eXJc31eLCoqCklJSYiKigKXy5V8GSh/iaBiQhYn74o15fItBHw+H5mZmdi5cycaN26Mpk2bVjkCHwBRE7urq2j59FPRuuJi4PZt4No1IC5O9PPmTfxSUIB7BQWI3LIFYVu2iLY1NQXatgXatZNeHBwUf1GJ3qGkTgipU4r23q9ue0U6jAUGBiI7OxuBgYHgcrkyvwzI2zRdvoUgLCwMhw8fRnJyMqKioir1I5DrfI2NRc3uPj7AxImidaWl8I+NxYYff0Qgnw/cvw/cuAG8fg1cvSpayrO3Bzw9/0vynp6AhwdATe31CnWUI4SolKId4erqWOW3BSDXfikpKVi9enWla/YV11e1XW1VOj+BAHj4ELh1S3pJSgKq+ih3dATc3UUJvvzi4lLl0LdEs2qTryipE0JUSp+SelhYGC5fvlypJ3/F46rrnOXuqf7qFXD3rnSiv3NHdP2+KiYmQKtWogTv7g60aAG4uYkWB4dKPfFJ3aHe74QQraHrg+OUp8g1e3WQu6e6mRnwzjuipbz8fFHnusREICHhv+XBA+DNG+DmTdFSkbk50Lz5f0m+/OLiIurYR7QS1dQJITpLkXuuy9emxdfUa9qvqhp4xePq3L3fAoGoB375RP/woWhJTa26KR8AeLz/Ovq5uPz3U/y4SROFB9Yh0qj5XQ6U1Amp35RJvHV5KUFrFBcDT578l+TLL48eiWr4NbG3l0724qVxY4DPF13nNzJS/7noKErqcqCkTghRlM7VwNVNKBRdp3/0CEhJESX/lJT/Hj95AhQVyVeWvb0owfP5/yX78kvjxqJteDz1npMWomvqhBCiBjT6WgVcrijZNm4s+3nGgJwc6UQv/pmaCqSlib4UlJaKprN98UJ0m15VeDxRp71GjURL+ceyfjc2Vs956xBK6oQQQlSDwxENd2tjIxrbXhahEHj5UpTgKy7Pnv33OCNDdO1f/Ls8rK2lE76dHWBrK1psbCo/trEB9Oyyin6dDSGEEO3G5Yqa1e3tAS+vqrcrKwMyM0U1+8zM/5aMDNm/l5UBeXmi5cED+eOxtpad8MWPra1FA/hYW1d+rIWXBiipE0II0T4GBv9dX68JY0BubuWE/+KFqFUgO1v0s/zj3FzRvuIvAsnJisdoYSGd6Cv+FD+2sgIsLWUv5uYqHROAkjohhBDdxuEADRuKFg8P+fYpKxMl9orJvvzj7GxRws/N/S/55+aKhuoFgMJC0fLsWe1it7CQTvRmZkoXR0mdEEJI/WNgILrmbmen+L4lJdJJvuLjij/z8oCCAumlsFDUwsDYf+tUcVoqKaWObNiwAd9//z2eP38OLy8vrFu3Dp07d9Z0WIQQQuoTI6P/+gUoizHR7X8Vk31BgegygnhiHwXpTFKPiYlBaGgoNm3ahC5dumDNmjUYOHAgEhMT0ahRI02HRwghhMhP3OxuYQE4OUk/l5+vdFLXmdEUVq1ahYkTJ8Lf3x9t2rTBpk2bYGZmhq1bt2o6NEIIIUQr6ERNvaSkBHFxcZg3b55kHZfLRb9+/XDhwgWZ+xQXF6O4uFjye15eHgDRSD2EEEKIthLnKXHeEjM2NoZxDQPs6ERSz8rKgkAggIODg9R6BwcHJCQkyNxn+fLlUlMlijVp0kQtMRJCCCGq5OLiIvV7WFgYvv3222r30Ymkrox58+YhNDRU8rtQKISPjw+uXbsGThX3BObn56NJkyZITU3V+/HhO3XqhCtXrmg6DLXHoarya1uOMvsrso+821a3Hf3/62cc2vAeUPf/v7zb17SNtrwHGGN4+fIlbGxspOYcqKmWDuhIUrezswOPx0NGRobU+oyMDDg6OsrcR1YzhbGxMaytrWs8npWVld5/qPF4PK04R3XHoarya1uOMvsrso+828qzHf3/61cc2vAeUPf/v7zby1umNrwH5MlVsuhERzkjIyP4+PggNjZWsk4oFCI2NhZdu3aVu5ygoCB1hKeTtOW1UHccqiq/tuUos78i+8i7rbb83TVNW16HuohDG94D6v7/l3d7bfm7q5POTL0aExOD8ePH4+eff0bnzp2xZs0a7N27FwkJCZWutSuLpmcl9Rn9/5P6Th/eAzrR/A4Avr6+ePHiBRYuXIjnz5/D29sbx44dU1lCB0TN82FhYXJdtyBE39D/P6nv9OE9oDM1dUIIIYRUTyeuqRNCCCGkZpTUCSGEED1BSZ0QQgjRE5TUCSGEED1BSV1OR44cgbu7O1q2bIlffvlF0+EQUudGjBiBhg0b4uOPP9Z0KITUudTUVPTu3Rtt2rRB+/btsW/fPk2HJBP1fpdDWVkZ2rRpg1OnTsHa2ho+Pj74999/YWtrq+nQCKkzp0+fRkFBAbZt24b9+/drOhxC6lR6ejoyMjLg7e2N58+fw8fHB/fv34e5ubmmQ5NCNXU5XL58GW3btkXjxo1hYWGBQYMG4cSJE5oOi5A61bt3b1haWmo6DEI0wsnJCd7e3gAAR0dH2NnZITs7W7NByVAvkvrZs2cxZMgQ8Pl8cDgcHDp0qNI2GzZsQNOmTWFiYoIuXbrg8uXLkufS0tLQuHFjye+NGzfGs2fP6iJ0QlSitu8BQnSdKt8DcXFxEAgEWjnrZ71I6kVFRfDy8sKGDRtkPh8TE4PQ0FCEhYXh2rVr8PLywsCBA5GZmVnHkRKiHvQeIPWdqt4D2dnZGDduHDZv3lwXYSuO1TMA2MGDB6XWde7cmQUFBUl+FwgEjM/ns+XLlzPGGDt//jwbPny45Plp06axnTt31km8hKiaMu8BsVOnTrGRI0fWRZiEqI2y74E3b96wHj16sO3bt9dVqAqrFzX16pSUlCAuLg79+vWTrONyuejXrx8uXLgAAOjcuTNu376NZ8+eobCwEH/++ScGDhyoqZAJUSl53gOE6DN53gOMMfj5+aFPnz74/PPPNRVqjep9Us/KyoJAIKg0MYyDgwOeP38OADAwMMCPP/6I//u//4O3tze++uor6vlO9IY87wEA6NevH0aNGoWjR4/C2dmZEj7RG/K8B86fP4+YmBgcOnQI3t7e8Pb2xq1btzQRbrV0ZpY2TRs6dCiGDh2q6TAI0ZiTJ09qOgRCNKZ79+4QCoWaDqNG9b6mbmdnBx6Ph4yMDKn1GRkZcHR01FBUhNQdeg+Q+k6f3gP1PqkbGRnBx8cHsbGxknVCoRCxsbHo2rWrBiMjpG7Qe4DUd/r0HqgXze+FhYVISkqS/J6cnIz4+HjY2NjAxcUFoaGhGD9+PN555x107twZa9asQVFREfz9/TUYNSGqQ+8BUt/Vm/eAprvf14VTp04xAJWW8ePHS7ZZt24dc3FxYUZGRqxz587s4sWLmguYEBWj9wCp7+rLe4DGfieEEEL0RL2/pk4IIYToC0rqhBBCiJ6gpE4IIYToCUrqhBBCiJ6gpE4IIYToCUrqhBBCiJ6gpE4IIYToCUrqhBBCiJ6gpE4IIYToCUrqhBBCiJ6gpE4IIYToCUrqhBBCqpSamorevXujTZs2aN++Pfbt26fpkEg1aEIXQgghVUpPT0dGRga8vb3x/Plz+Pj44P79+zA3N9d0aEQGqqkTrdK7d29Mnz5d02HUCX07V02ez8uXL9GoUSM8fvxYI8fXZ05OTvD29gYAODo6ws7ODtnZ2ZLnP/30U/z4448aio5UREmdKOXFixeYMmUKXFxcYGxsDEdHRwwcOBDnz58HoBsJi8PhVLt8++23cpWjznP18/PD8OHDFdpHU6/9gQMHsHjxYpXH0KtXLwQEBFS7zdKlSzFs2DA0bdq0VsdS5vWuT+Li4iAQCNCkSRPJuvnz52Pp0qXIy8vTYGREzEDTARDdNHLkSJSUlGDbtm1o3rw5MjIyEBsbi5cvX2o6NLmlp6dLHsfExGDhwoVITEyUrLOwsNBEWDrLxsZG5WUyxnD9+nWMGjWqym1evXqFyMhIHD9+XOXHrwslJSUwMjLSaAze3t4oKyurtP7EiRPg8/kAgOzsbIwbNw5btmyR2sbT0xNubm7YsWMHgoKC6iReUg1GiIJycnIYAHb69GmZz48fP54BkFqSk5MZY4y9efOGBQcHM3t7e2ZsbMy6devGLl++LNm3V69ebNq0aZLfjxw5wqysrNiOHTsYY4wJBAK2bNky1rRpU2ZiYsLat2/P9u3bJ3X8ffv2MU9PT2ZiYsJsbGxY3759WWFhYbXnFBUVxaytrWU+V13M1Z3rn3/+ybp168asra2ZjY0NGzx4MEtKSqryXKt6LYcNGya1T3BwMJs1axZr2LAhc3BwYGFhYVLby4qnptetpnLleV3F51NVDNu2bWM2NjbszZs3UuUOGzaMjR07Vub5JyYmMgDs33//rfI12rdvH7O3t5da16tXLzZ16lQ2bdo01qBBA9aoUSO2efNmVlhYyPz8/JiFhQVzc3NjR48erfL1ViZexuT7uwcFBbFp06YxW1tb1rt3b8aY6H97xYoVzM3NjRkZGbEmTZqwJUuWSM6xute+pr9vdWXL482bN6xHjx5s+/btMp9ftGgR6969u9zlEfWhpE4UVlpayiwsLNj06dMrfeAxxlhubi7r2rUrmzhxIktPT2fp6emsrKyMMcZYSEgI4/P57OjRo+zOnTts/PjxrGHDhuzly5eMMelEt3PnTmZpackOHz4sKXvJkiXMw8ODHTt2jD18+JBFRUUxY2NjyReMtLQ0ZmBgwFatWsWSk5PZzZs32YYNG1hBQUG151RdUq8u5urOdf/+/ex///sfe/DgAbt+/TobMmQIa9euHRMIBJXOtSqykrqVlRX79ttv2f3799m2bdsYh8NhJ06cqPa1r+l1q6lceV5X8flUFcOrV6+YtbU127t3r2SfjIwMZmBgwP7++2+Z579r1y7G4/FYUVFRla9RSEgIe//996XW9erVi1laWrLFixez+/fvs8WLFzMej8cGDRrENm/ezO7fv8+mTJnCbG1tpcou/3orEy9j8v3dLSws2KxZs1hCQgJLSEhgjDE2e/Zs1rBhQxYdHc2SkpLYuXPn2JYtW+R67Wv6+1ZVtjyEQiH79NNPK33JK+/PP/9kRkZGMj8PSN2ipE6Usn//ftawYUNmYmLC3nvvPTZv3jx248YNyfOyElZhYSEzNDRkO3fulKwrKSlhfD6frVy5Umq/9evXM2tra6nWgDdv3jAzM7NKtbbAwEA2evRoxhhjcXFxDAB7/PixQudTVVJXJOaavHjxggFgt27dkns/WUm9Yo2oU6dObM6cOVLblC9XntetpnLleV3LH7eqc5syZQobNGiQ5Pcff/yRNW/enAmFQpllzpw5k7Vp06bKYzImqjkHBARUiqX8+ZSVlTFzc3P2+eefS9alp6czAOzChQuSdRVfb0XjlUXW371Dhw5S2+Tn5zNjY2OZibam176mv291Zcvj3LlzjMPhMC8vL8ly8+ZNqW1u3Lih1PuOqB5dUydKGTlyJAYPHoxz587h4sWL+PPPP7Fy5Ur88ssv8PPzk7nPw4cPUVpaim7duknWGRoaonPnzrh3755k3f79+5GZmYnz58+jU6dOkvVJSUl49eoV+vfvL1VuSUkJOnToAADw8vJC37590a5dOwwcOBADBgzAxx9/jIYNGyp1nvLGLMuDBw+wcOFCXLp0CVlZWRAKhQCAlJQUeHp6KhUPALRv317qdycnJ2RmZla5vTyvW03lqup1nThxIjp16oRnz56hcePGiI6Ohp+fHzgcjsztr127Bh8fn2rLfP36NUxMTCqtL38+PB4Ptra2aNeunWSdg4MDAFT72ikaLyDf373iOd27dw/FxcXo27dvpfJqeu1r+vtWV7Y8unfvLjmHqpiamgIQ9W8gmkVJnSjNxMQE/fv3R//+/bFgwQJMmDABYWFhVSZ1eXXo0AHXrl3D1q1b8c4770g+QAsLCwEAf/zxBxo3biy1j7GxMQDRh/dff/2Ff//9FydOnMC6devwzTff4NKlS2jWrFmt4lLUkCFD4Orqii1btoDP50MoFMLT0xMlJSW1KtfQ0FDqdw6HU+2HrjyvW03lqup17dChA7y8vLB9+3YMGDAAd+7cwR9//FHl9teuXcOQIUOqLdPOzg45OTmV1ss6n/LrxP9X1b12isYLyPd3r3iPtzgpylLTa1/T3zc3N7faeFVBfIubvb292o9Fqke3tBGVadOmDYqKigAARkZGEAgEUs+7ubnByMhIctsbAJSWluLKlSto06aN1HanTp3Cb7/9huDgYKnyjY2NkZKSghYtWkgt5W+x4XA46NatGxYtWoTr16/DyMgIBw8eVOqc5IlZ1rm+fPkSiYmJmD9/Pvr27YvWrVvLTDzqUDEeeV+3mijyusp6TcQmTJiA6OhoREVFoV+/flXG8OjRI+Tm5tZYU+/QoQPu3r0r93koSt54AeX/7i1btoSpqSliY2NlPl/da1/T37emslXh9u3bcHZ2hp2dndqOQeRDNXWisJcvX2LUqFEICAhA+/btYWlpiatXr2LlypUYNmwYAKBp06a4dOkSHj9+DAsLC9jY2MDc3BxTpkzBrFmzYGNjAxcXF6xcuRKvXr1CYGCg1DFatWqFU6dOoXfv3jAwMMCaNWtgaWmJmTNnYsaMGRAKhejevTvy8vJw/vx5WFlZYfz48bh06RJiY2MxYMAANGrUCJcuXcKLFy/QunVrpc5VnphlnWvDhg1ha2uLzZs3w8nJCSkpKZg7d27tXng5yYqnptetJoq+rrJi4HJFdYjPPvsMM2fOxJYtW7B9+/YqjxkXFwcOhyMZ+KQqAwcOxLx585CTk6P0ZZbqyBsvAKX/7iYmJpgzZw5mz54NIyMjdOvWDS9evMCdO3fg6elZ7Wsvz/uiqrIrvu+Ude7cOQwYMEAlZZFa0vRFfaJ73rx5w+bOncs6duzIrK2tmZmZGXN3d2fz589nr169YoyJbkV69913mampqdRtXq9fv2bBwcHMzs5Orlva7t69yxo1asRCQ0MZY6KeuGvWrGHu7u7M0NCQ2dvbs4EDB7IzZ85Ith84cKDk9rNWrVqxdevW1XhO1fV+rynmqs71r7/+Yq1bt2bGxsasffv27PTp0wwAO3jwoMxzlUVWR7mK+wwbNoyNHz++2nhqet1qKlee17V8GVW9JmKff/65zNvFyps7dy5r1apVta+PWOfOndmmTZtkxiLm6urKVq9eLbWu/N9DHNfIkSMrlS9PvGLK/t0FAgFbsmQJc3V1ZYaGhszFxYUtW7ZMrte+pr9vVWWrwuvXr5m1tbVUh0OiOTT2OyGkzvXt2xdt27ZFRESESsr7448/MGvWLNy+fVvSIqCM999/Hy1atMD69eul1qs6Xn2yceNGHDx4ECdOnNB0KATU/E4IqUM5OTk4ffo0Tp8+jZ9++kll5Q4ePBgPHjzAs2fPFOonUD6u8+fP4/Tp05g8ebLa49UnhoaGWLdunabDIG9RTZ0QUmeaNm2KnJwcLFiwADNnztR0OBIjRozAlStXMH78eCxZskTSM15b4yWkKpTUCSGEED1Bt7QRQggheoKSOiGEEKInKKkTQggheoKSOiGEEKInKKkTQggheoKSOiGEEKInKKkTQggheoKSOiGEEKInKKkTQggheoKSOiGEEKInKKkTQggheoKSOiGEEKIn/h8TO82jskHP9wAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "PvI_plot_Step2 = HtoH.PvI(CatP, CatdP, CatI, Iscale=[0.9,150.0], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "59006181",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Power law index: -0.6029396034448222 ± 0.07562184472476109\n",
      "Coefficient: 126812861965214.86 ± 489213661705306.2\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgAAAAF2CAYAAAAY6yC7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABpnUlEQVR4nO3deVxU9frA8c8MCLgAiiCLIO5b7rhkapmZZqtp/sxKBbFuZW5kli2ampHeytRMb4ZYN02z1CzNFkotc19z31BQFjOFERcU5vz+OHcmdmaGM8wMPO/X67x0Dud855mBmfOc76pTFEVBCCGEEJWK3tEBCCGEEKL8SQIghBBCVEKSAAghhBCVkCQAQgghRCUkCYAQQghRCUkCIIQQQlRCkgAIIYQQlZAkAEIIIUQlJAmAEEIIUQlJAiCEC8jOzubNN98kOzvb0aEIIeysvD7vTpEALFiwgDZt2uDj44OPjw9du3bl+++/N//8xo0bjBo1itq1a1OjRg0GDhxIenq6AyMWonxlZ2czdepUSQCEqATK6/PuFAlAaGgo77zzDrt372bXrl306tWLRx55hEOHDgEwfvx4vv32W1auXMmmTZtISUlhwIABVj/P/PnztQ7drhwdr72fX+vytSivLGXYcq6jf8cVmau9t46OVz7v9j/X0b/jQhQnVatWLeWTTz5RMjIylCpVqigrV640/+zIkSMKoGzdutWqMlu0aKF1mHbl6Hjt/fxal69FeWUpw5ZzLT0nMzNTAZTMzEyrn6OycvTnx1qOjlc+7/Y/19k+7+4Ozj8Kyc3NZeXKlVy9epWuXbuye/dubt26Re/evc3HNG/enHr16rF161Zuv/32IsvJzs7OV31iNBrJzs4mMzMTnU5n99ehhdzcXAwGQ4V9fq3L16K8spRhy7mWnmM6xpF/D67G0Z8fazk6Xvm82/9cS8/JzMwEICMjI99+T09PPD09rXrOEtk1vbDCgQMHlOrVqytubm6Kr6+vsm7dOkVRFGXp0qWKh4dHoeM7deqkTJw4sdjypkyZogCyySabbLLJViG2KVOmaHrddZoagGbNmrFv3z4yMzP56quvGD58OJs2bbK5vEmTJhETE2N+nJmZSb169UhOTsbHx0eLkMX/5OTkcP78eerWrYu7e+l/UtYebzonKSkJgHr16ll8nhBCuBqDwUBYWBhJSUn4+vqa92t69w84zbeoh4cHjRs3BiAiIoKdO3cyZ84cBg8ezM2bN8nIyKBmzZrm49PT0wkKCiq2vOKqSkwjDYR2cnJyMBgM+Pj4WJwAWHO86Rxvb28Aq84TQghX5evra9frlVOMAiiKqc0+IiKCKlWqkJCQYP7ZsWPHSEpKomvXrg6MUAghhHBdTnEbNWnSJPr160e9evW4cuUKy5YtY+PGjfzwww/4+voSHR1NTEwMfn5++Pj4MHr0aLp27VpsB0AhhBBClMwpEoALFy4wbNgwUlNT8fX1pU2bNvzwww/ce++9AMyePRu9Xs/AgQPJzs6mb9++fPTRRw6OWgghhHBdOkVRFEcHUR4MBgO+vr5kZmZKHwCN5eTkcO7cOUJDQy3uA2DN8aZzzp49C0B4eLj0ARBCVFjldb1y2j4AQgghhLAfSQCEEEKISkgSACGEEKISkgRACCGEqIQkARBCCCEqIUkAhBBCiEpIEgAhhBCiEpIEQAghhKiEJAEQQgghKiFJAIQQQohKSBIAIYQQohKSBECUu5SUFObMmUNKSoqjQxFCiEpLEgBR7uLi4jh69ChxcXGODkUIISotWVJNlLvo6GgMBgPR0dGODkUIISotqQEQ5S4kJISxY8cSEhLi6FCEEKLSkgRACCGEqIQkARBCCCEqIUkARIVhNBrJyMjAaDQ6RTlCVARLlixBp9Oh0+nM+3r27IlOpyMyMtJxgYkykwRAuISUlBTmzp1Lenp6scdIAiCEdUwX8qK2NWvWABAQEECXLl3o0qWLY4PVUGRkJDqdjp49e1p8zpo1a7jzzjvx9vamatWqNGnShHfeecd+QZYDGQUgXEJcXBz79+9n+/btfPnll9SrV8/RIQlRYXh4eNC+fft8+/z8/AB44IEHeOCBBxwRltN47733mDBhAgBBQUEEBwdz4cIFEhISeOWVVxwcne2kBkC4hOjoaG7cuIGXl5fMHyCExoKDg9m2bVu+7c477wSKbgIoyvDhw2nSpAne3t54eHgQHh7OmDFjMBgM5mNMd97169fn008/JSwsDB8fH8aNG8f169cZN24cvr6+1KtXj4ULF+YrPyUlhREjRhASEoKHhwcNGzZk+vTp5OTkmI8x1WgMGzaMKVOmEBwcTK1atXjqqae4cuUKgPm5ATZt2mR+bRs3bizydSUnJ5sv8nPnziUlJYU9e/Zw7tw5Vq1aZd0b7WQkARAuISQkhPfff5+2bdvK/AFCOKFvvvmGy5cv06hRI8LCwkhKSmLevHlFfl5TU1N5/vnn8fT05MqVK8yZM4eIiAji4+Px9fUlOTmZUaNGceTIEQD+/vtvbr/9duLj48nKyqJFixYkJyczefJknnnmmULlL1++nNmzZ1O1alUyMjJYunSpubq+ffv2+Pv7A+Dt7W1u3vDx8Snyda1atYqcnByqV6/Otm3b8Pf3Jzg4mKFDh3L16lWt3j6HkARAuIzAwEDGjBkj8wcIobGzZ88W6gNgrU2bNnHx4kX27dvHqVOneO211wC17fzGjRv5jr158yY//vgjx48fJywsDICTJ0+yf/9+Dh8+jJeXF0ajkU2bNgHw4YcfkpycTGBgIKdOnWL//v189dVXgFpDcfLkyXzle3l5ceTIEU6ePElERAQACQkJAKxevdrcpNGhQwdzjUeHDh2KfF3Hjh0D4OrVq6xcuZLg4GD+/vtvPv/8c+6//35u3bpl9XvlLCQBEEKISs7Dw8N8J2xrh7+ff/6ZVq1aUbVqVXQ6HTNmzAAgJyeHv/76K9+xtWrVolu3buj1enN/nlatWlG/fn1q1KhBnTp1AMydfnfs2GF+XKdOHXQ6Hf379wdAURS2b9+er/xevXpRt25d9Ho9zZs3z1eWtfI2MSxevJiDBw+amyH37t3Lli1bbCrXGUgnQCGEqORMfQBstXTpUnMnueDgYMLCwrh48SKnT58GIDc3N9/xeavb3d3dC+0z1UAoipLvPG9vb1q2bFno+atVq5bvcc2aNQuVX7AsS9WtW9f8/06dOgHQuXNn874zZ87YVK4zkBoAIYQQZWJKHry9vUlMTGT79u306dNHs/JNF153d3eWL19urrb/6aefeP7553n00UetKs+UMFjSht+7d2/z/3ft2pXvX4AmTZpY9dzORBIAIYQQZdKmTRsArly5QsOGDWnYsCFffvmlZuWPGjWKunXrcvnyZZo1a0a7du1o1KgRtWvXZvjw4VaXZ2oW2LVrF61bt+b222/n+vXrRR7brVs3HnnkEQCioqJo3bo1UVFRANxzzz1069bNxlfleJIACCGEKJPo6GhiYmLw9/fnypUr9OzZk2nTpmlWfkBAANu2bSMqKoratWtz6NAhrl+/To8ePZg9e7bV5Y0YMYKBAwfi6+vLwYMH2b59e6FmiryWL1/Oyy+/TFBQECdOnKBBgwa88cYbfPvtt2V5WQ6nU2xtGHExBoMBX19fMjMzix3uIWyTlJTE7NmzGT9+vEUT9OTk5HDu3DlCQ0PN7XOWnHP27FkAwsPDizzPlnK1ik8IIbRSXtcrqQEQZRYXF8fRo0dlgh4hhHAhcnsjyiw6OhqDwSAT9AghhAuRGgBRZiEhIYwdO1Ym6BFCCBciCYAQQghRCUkCIIQQQlRCkgAIIYQQlZAkAEIIIUQlJAmAEEIIUQlJAiCEEEJUQk6RAMTGxtKpUye8vb2pU6cO/fv3N6/BbNKzZ89C61U/++yzDopYCCEqpvnz51O/fn28vLzo0qWLeSnekpw/f56nnnqK2rVrU7VqVVq3bp1vwRxbyy2v+Es7Jzc3lzfeeIMGDRpQtWpVGjVqxPTp021eYdBZOEUCsGnTJkaNGmVe3enWrVv06dOn0EpNTz/9NKmpqeZt1qxZDopYCCEqnhUrVhATE8OUKVPYs2cPbdu2pW/fvly4cKHYcy5fvky3bt2oUqUK33//PYcPH+a9996jVq1aZSrXpGfPnixZssRu8VtyzsyZM1mwYAEffvghR44cYebMmcyaNYt58+ZZFJfTUpzQhQsXFEDZtGmTed9dd92ljB071uYyMzMzFUDJzMzUIEKR161bt5TExETl1q1bdjnedM7JkyeVkydPFnueLeXasxwhTJo0aaLcfvvtyrVr18z7jEaj0qVLF+WVV15xYGT5de7cWRk1apT5cW5urhISEqLExsYWe87LL7+sdO/eXfNyTe666y4lPj6+9OBtfB5LznnggQeUESNG5DtvwIABypNPPmlRXNYqr+uVU9QAFJSZmQmAn59fvv1Lly7F39+fVq1aMWnSJK5du1ZsGdnZ2RgMhnybEEI4wooVK9izZw9btmwx71u6dClnz57l1VdfNe97++23qVGjRolbUlKSXWK8efMmu3fvpnfv3uZ9er2e3r17s3Xr1mLPW7t2LR07dmTQoEHUqVOH9u3bs2jRojKXWx7xW3rOHXfcQUJCAsePHwdg//79/P777/Tr10+z+ItS8BqWnZ2taflOtxaA0Whk3LhxdOvWjVatWpn3P/HEE4SHhxMSEsKBAwd4+eWXOXbsGKtWrSqynNjYWKZOnVpeYQshRLHat29Pu3btOHr0KL179+batWtMmjSJt956C29vb/Nxzz77LP/3f/9XYln2mnL74sWL5ObmEhgYmG9/YGAgR48eLfa806dPs2DBAmJiYnj11VfZuXMnY8aMwcPDg+HDh1td7ttvv83bb79tfnz9+nW2bdvGCy+8YN53+PDhQiuP2hK/pee88sorGAwGmjdvjpubG7m5ucyYMYMnn3yy2PdFC2FhYfkeT5kyhTfffFOz8p0uARg1ahQHDx7k999/z7f/mWeeMf+/devWBAcHc88993Dq1CkaNWpUqJxJkyYRExNjfmwwGAq9mUIIUV6aNm1q7tw8a9Ys/P39iYqKyneMn59foZpPa73yyivMnDmzxGOOHDlC8+bNy/Q8JkajkY4dO5ov2u3bt+fgwYMsXLiQ4cOHW11ewSToySefZODAgQwYMMC8r7zXHfnyyy9ZunQpy5Yt47bbbmPfvn2MGzeOkJAQm16jpZKTk/MtB+zp6alp+U6VALzwwgt89913bN68mdDQ0BKP7dKlCwAnT54sMgHw9PTU/M0SQghbNWvWjM2bN3Pu3Dn+/e9/s27dOvT6/K2wBe9+i1LU3W9eL774IpGRkSWW0bBhw0L7/P39cXNzIz09Pd/+9PR0goKCii0rODiYli1b5tvXokULvv76a5vKLZgEVa1alTp16tC4ceMSX5Mt8Vt6zksvvcQrr7zC448/Dqg3oWfPniU2NtauCYCPj0++BEBrTpEAKIrC6NGjWb16NRs3bqRBgwalnrNv3z5A/eMTQghn17RpUxYtWsQrr7xCnz596NmzZ6FjtGgCCAgIICAgwOr4PDw8iIiIICEhgf79+wPq3X1CQkK+6veCunXrVmjY9vHjxwkPDy9TueURv6XnXLt2rVCy5ubmhtFo1Cx+h7BrF0MLPffcc4qvr6+yceNGJTU11byZesyePHlSmTZtmrJr1y4lMTFR+eabb5SGDRsqd955p8XPIaMA7EdGAQhRur179yo6nU7x9PRUTpw44ehwirR8+XLF09NTWbJkiXL48GHlmWeeUWrWrKmkpaWZj5k3b57Sq1cv8+MdO3Yo7u7uyowZM5QTJ04oS5cuVapVq6Z8/vnnVpVrcuXKlXzXgaK2nJwczeK35Jzhw4crdevWVb777jslMTFRWbVqleLv769MnDjRtje6FOV1vXKKBAAocjMN/UhKSlLuvPNOxc/PT/H09FQaN26svPTSS1a9OZIA2I8kAEKU7urVq4pOp1NefPFFR4dSonnz5in16tVTPDw8lM6dOyvbtm3L9/MpU6Yo4eHh+fZ9++23SqtWrRRPT0+lefPmyscff2x1uXnLL+6aYNoSExM1jb+0cwwGgzJ27FilXr16ipeXl9KwYUPltddeU7Kzs4uNoyzK63qlUxQXn8rIQgaDAV9fXzIzM+3aplIZ5eTkcO7cOUJDQ3F3L71VydrjTeecPXsWgPDw8CLPs6VcreITojSXLl2idu3a7N+/nzZt2jg6HOHEyut65ZTzAAghREWzf/9+PDw8aNGihaNDEQIoYyfAW7dukZaWxrVr1wgICCjz8BUhhKio9u/fT8uWLalSpYqjQxECsKEG4MqVKyxYsIC77roLHx8f6tevT4sWLQgICCA8PJynn36anTt32iNWIYRwWePGjWPv3r2ODkMIM6sSgPfff5/69esTHx9P7969WbNmDfv27eP48eNs3bqVKVOmkJOTQ58+fbjvvvs4ceKEveIWQgghRBlY1QSwc+dONm/ezG233Vbkzzt37syIESNYuHAh8fHx/PbbbzRp0kSTQIUQQgihHasSgC+++MKi4zw9PXn22WdtCkgIIYQQ9qfJGKeLFy+yfft2cnNz6dSpk8zOJ4QQQji5MicAX3/9NdHR0TRt2pRbt25x7Ngx5s+fX2iRCyGEEEI4D6tHAWRlZeV7PHXqVHbs2MGOHTvYu3cvK1eu5LXXXtMsQCGEEEJoz+oEICIigm+++cb82N3dnQsXLpgfp6en4+HhoU10QgghhLALq5sAfvjhB0aNGsWSJUuYP38+c+bMYfDgweTm5pKTk4Ner2fJkiV2CFUIIYQQWrE6Aahfvz7r1q3jiy++4K677mLMmDGcPHmSkydPkpubS/PmzfHy8rJHrEJYxGg0kpGRgY+PT6ElPIUQQqhs/nYcMmQIO3fuZP/+/fTs2ROj0Ui7du3k4i8czpQAuPxa3UIIYUc2jQJYv349R44coW3btnzyySds2rSJJ598kn79+jFt2jSqVq2qdZxCCCGE0JDVNQAvvvgiUVFR7Ny5k3/9619Mnz6du+66iz179uDl5UX79u35/vvv7RGrEEIIITRidQKwZMkS1q9fz/Lly9m5cyf//e9/AfDw8GD69OmsWrWKt99+W/NAhSgvKSkpzJkzh5SUFEeHIoQQdmN1AlC9enUSExMBSE5OLtTm37JlS3777TdtohPCAeLi4jh69ChxcXGODkUIIezG6j4AsbGxDBs2jDFjxnDt2jU+/fRTe8QlhMNER0djMBiIjo52dChCCGE3VicATz75JPfddx+nT5+mSZMm1KxZ0w5hCeE4ISEhjB07lpCQEEeHIoQQdmNVE0BSUhIAtWvXplOnTqVe/M+fP29zYEIIIYSwH6sSgE6dOvGvf/2LnTt3FntMZmYmixYtolWrVnz99ddlDlAI4TxkjgUhKg6rmgAOHz7MjBkzuPfee/Hy8iIiIoKQkBC8vLy4fPkyhw8f5tChQ3To0IFZs2Zx//332ytuIYQDmBKAGjVqyCyLQrg4qz7BtWvX5v333yc1NZUPP/yQJk2acPHiRU6cOAGo/QN2797N1q1b5eIvhBBCODGbZgKsWrUqjz32GI899pjW8QhRJL1ej6+vLwaDwdGhCCFEhSB1eMIl6PV6atasKdXOQgihEfk2FUIIISohSQCEEEKISkgSACGEEKISsqkT4NmzZzlw4ACBgYF07txZ65iEEEIIYWdWJwBffPEFkZGR3Lp1C51OZ17+NyAgwB7xCSGEEMIOrG4CmDp1Kk888QRHjx7lxx9/BOCVV17RPDBRcUmPfiGEcDyrawBOnz7Nhg0bqF+/Pk2bNuXzzz8nIiJClk4VFjMlAEIIIRzH6luwnJwcqlWrZn7cvHlzjEYjaWlpmgYmhLVSUlKYM2cOKSkpjg5FCCGcnk11sJ9++il//PEHWVlZALi7u3Pt2jVNAxPCWnFxcRw9epT4+HhHhyKEEE7P6iaAHj168NZbb3HlyhX0ej0NGjTgxo0bxMXF0bt3bzp27Ii3t7c9YhWVXGl9B6KjozEYDERFRZVzZEII4Xp0iqIotpx44sQJdu/ezZ49e8xbRkYGer2eJk2acOTIEa1jLRODwYCvry+ZmZn4+Pg4OpwKJScnh3PnzhEaGoq7u00jSzWNIygoiLS0NJvjcZbX44zkvRHC/srremXzJ7hJkyY0adKExx9/3LwvMTGRXbt2sXfvXk2CE0IIIYR9aJrCN2jQgAYNGjBo0CAtixXCYqb16o1Go6NDEUIIp+YUA7FjY2Pp1KkT3t7e1KlTh/79+3Ps2LF8x9y4cYNRo0ZRu3ZtatSowcCBA0lPT3dQxMJZKYqCwWDAxpYtIYSoNJwiAdi0aROjRo1i27Zt/PTTT9y6dYs+ffpw9epV8zHjx4/n22+/ZeXKlWzatImUlBQGDBjgwKiFEEII1+UUvXg2bNiQ7/GSJUuoU6cOu3fv5s477yQzM5O4uDiWLVtGr169AIiPj6dFixZs27aN22+/vVCZ2dnZZGdnmx8bDAb7vgghhBBCQwWvW56ennh6empWvlPUABSUmZkJgJ+fHwC7d+/m1q1b9O7d23xM8+bNqVevHlu3bi2yjNjYWHx9fc1bWFiY/QMXQgghNBIWFpbvOhYbG6tp+WWqAUhISCAhIYELFy4U6nS1ePFim8o0Go2MGzeObt260apVKwDS0tLw8PAoNH1sYGBgsTMQTpo0iZiYGPNjg8EgSYAQQgiXkZycnG8YoJZ3/1CGBGDq1KlMmzaNjh07EhwcjE6n0ySgUaNGcfDgQX7//fcylaN1VYkQQghRnnx8fJxzHoCFCxeyZMkShg4dqlkwL7zwAt999x2bN28mNDTUvD8oKIibN2+SkZGRrxYgPT2doKAgzZ5fCK0ZjUYMBgM+Pj6y+qEQwqnY/I108+ZN7rjjDk2CUBSFF154gdWrV/PLL7/QoEGDfD+PiIigSpUqJCQkmPcdO3aMpKQkunbtqkkMQtiDzEsghHBWNicAI0eOZNmyZZoEMWrUKD7//HOWLVuGt7c3aWlppKWlcf36dQB8fX2Jjo4mJiaGX3/9ld27dxMVFUXXrl2LHAEghBBCiJLZ3ARw48YNPv74Y37++WfatGlDlSpV8v38/ffft7isBQsWANCzZ898++Pj44mMjARg9uzZ6PV6Bg4cSHZ2Nn379uWjjz6yNXwhhBCiUrM5AThw4ADt2rUD4ODBg/l+Zm2HQEtmbfPy8mL+/PnMnz/fqrKF/ZW2Sp8QQgjnY3MC8Ouvv2oZh3BhpgRACCGE65BbNiGEEKISKtNEQBkZGcTFxXHkyBEAWrZsSXR0NL6+vpoEJ4QQQgj7sLkGYNeuXTRq1IjZs2dz6dIlLl26xOzZs2nUqBF79uzRMkYhhBBCaMzmGoDx48fz8MMPs2jRItzd1WJycnIYOXIk48aNY/PmzZoFKYQQQght2ZwA7Nq1K9/FH8Dd3Z2JEyfSsWNHTYITQgghhH3Y3ATg4+NDUlJSof3Jycl4e3uXKSghhBBC2JfNCcDgwYOJjo5mxYoVJCcnk5yczPLlyxk5ciRDhgzRMkYhbCLT8AohRPFsbgJ499130el0DBs2jJycHACqVKnCc889xzvvvKNZgELYypQA1KhRQyYpEkKIAmxOADw8PJgzZw6xsbGcOnUKgEaNGlGtWjXNgrOLixfBjssrCiGEEK6gTPMAAFSrVo3WrVtrEUv5uO8++PVXqFvX0ZEIJ2Nautfpk1ghhNCAVQlATEwM06dPp3r16sTExJR4rDWLAZWrEyegRw9ISIACyw6Lys3UZODl5eXoUIQQwu6sSgD27t3LrVu3zP8vjrWLAZWrBg0gMVFNAn7+GZo3d3REFZ7pztrHx0fa4oUQwklYlQDkXQDo008/JTQ0tNAXuqIoJCcnaxOdPXz/PQwYAIcPw513wo8/wv9WNRT2IZ3xhBDC+dj8bdygQQMuXrxYaP+lS5do4MxV68HBsGkTdOgAf/0Fd98N27Y5OirhRGR5YyFEZWDzN5yiKEXuz8rKcv42VH9/+OUX6NYNMjKgd2/1sRBIAiCEqBysHgVg6vyn0+mYPHlyvh7Tubm5bN++nXauUKXu6ws//AD9+6t9Ae6/H77+Gh54wNGRiTJKT0/nq6++YuzYsXIRF0KIYlidAJg6/ymKwp9//omHh4f5Zx4eHrRt25YJEyZoF6E9Va8O334LgwfD2rVqMrBsGQwa5OjIhI3S09OZMGECVatWJT4+nujoaEeHJIQQTsnqBMDUETAqKoo5c+bg4+qT6nh5wVdfwfDh8MUX8PjjkJUFUVGOjkzYYMWKFXh5eXHjxg2i5HconIiMhhHOxuaJgOLj47WMw7GqVIH//letEfjkExgxAi5fhlLmOhDOZ/DgwSiKwuOPP05ISAhpaWmODkkIQEbDCOdj819hbGwsixcvLrR/8eLFzJw5s0xBOYSbG3z88T8X/RdfhAkTQBaScSmBgYG88MILBAYGOjoUIYRwajYnAP/5z39oXsQkOrfddhsLFy4sU1AOo9PBu+/CrFnq4/feg2HD4OZNx8ZVycmqfs5Bfg9CVCw2JwBpaWkEBwcX2h8QEEBqamqZgnIonQ5eegk++wzc3WHpUnjwQbhyxdGRVVr2vPDIRc1yRqORzMzMYocACyFci80JQFhYGFu2bCm0f8uWLYSEhJQpKKcwdKg6QqB6dfjpJ+jZE9LTHR2V0JgkAEKIysrmBODpp59m3LhxxMfHc/bsWc6ePcvixYsZP348Tz/9tJYxOo5p5cCAANizB+64A06edHRUQjhMeno6c+fOJSUlxdGhCCHKyOZRAC+99BJ///03zz//PDdv3kRRFKpWrcrLL7/MK6+8omWMjtWpE2zZAn37wunT6uyB69dDRISjI6v0ZFiV5bR6r1asWMGxY8eIi4tj6tSpGkYohChvNn8T6HQ6Zs6cyV9//cW2bds4cOAAly5dYvLkyc69GqAtmjSBP/6A9u3hwgW46y51ESHhUFJ9bzmt3qvBgwfTrFkzmWBJiArA5hoAk6SkJP7++29u3rzJ6dOnzfsffvjhshbtXIKC1EWEBgxQpw5+4AGIj4ennnJ0ZEKUm8DAQMaMGVMx+vkIUcnZnACcPn2aRx99lD///BOdTmfuGWy6+8/NzdUmQmfi7Q3r1kFkpDpr4NChcOoUTJ6sjh4QQgghXITNTQBjx46lQYMGXLhwgWrVqnHw4EE2b95Mx44d2bhxo4YhOhkPD/j8c5g4UX385pvqXAHZ2Q4NSwghhLCGzQnA1q1bmTZtGv7+/uj1etzc3OjevTuxsbGMGTNGyxidj14PM2eqMwe6uakJwb33wt9/OzoyIYQQwiI2JwC5ubl4e3sD4O/vbx4WFB4ezrFjx7SJztk9/TRs2KAuLfzbb3D77XD8uKOjEkIIIUplcwLQqlUr9u/fD0CXLl2YNWsWW7ZsYdq0aTRs2FCzAJ1e797qCIH69dU5Am6/Xe0sKIQQQjgxmxOA119/3TykaNq0aSQmJtKjRw/Wr1/P3LlzNQvQJbRsCdu3qxf/y5fV5oDPPnN0VEIIIUSxbB4F0LdvX/P/GzduzNGjR7l06RK1atWqePMAWKJOHfjlF3WEwJdfwvDhcOIETJsmIwSEsJFM9iSE/dj0ibp16xb33HMPJ06cyLffz8+vcl78TapWVYcHvvqq+vitt2DwYLh61bFxCeGiZLInIezHpgSgSpUqHDhwQOtYKga9HmbMUCcJqlIFVq6E7t3h7FlHRyaEEEKY2Vyn9tRTTxEXF6dJEJs3b+ahhx4iJCQEnU7HmjVr8v08MjISnU6Xb7vvvvs0eW67iYxUmwTq1IF9+9Q1BTZvdnRUogQpKSnMmTNH04Vu7FGmEEJoweYEICcnhwULFtCxY0f+9a9/ERMTk2+zxtWrV2nbti3z588v9pj77ruP1NRU8/bFF1/YGnr56d4ddu5U1xD46y+45x74z38cHZUoRlxcHEePHtUssbVXmUIIoQWbOwEePHiQDh06AHC8wNh3a/sB9OvXj379+pV4jKenJ0FBQdYF6Qzq1YPff4cRI2DFCnj2Wdi/H+bMUZsIhFMwGo0MGjSIzMxMTRe6iY6OxmAwyOI5QginY3MC8Ouvv2oZR6k2btxInTp1qFWrFr169eKtt96idu3axR6fnZ1Ndp7peQ0GQ3mEWbRq1dTOgW3bwmuvwYIFcOgQfPUVBAQ4Li5hZjQaqVq1KqNGjaJatWoYjUZNep2HhIQwduxYWTxHCGG1gtctT09PPD09NSvf6m+406dPmxf+KS/33Xcfn332GQkJCcycOZNNmzbRr1+/Ehccio2NxdfX17yFhYWVY8RF0Olg0iT45ht1UaHNm9V+Afv2OTYukY+iKGRmZhbqdS690YUQ5S0sLCzfdSw2NlbT8q1OAJo0acJff/1lfjx48GDS09M1Daqgxx9/nIcffpjWrVvTv39/vvvuO3bu3FniokOTJk0iMzPTvCUnJ9s1Ros99BBs2waNG6sjA+64A/77X0dHJUohCYAQorwlJyfnu45NmjRJ0/KtTgAK3v2vX7+eq+U8zr1hw4b4+/tz8uTJYo/x9PTEx8cn3+Y0WraEHTvgvvvg+nV1NcFRo+DmTUdHJoQQwkkUvIZpWf0PZRgF4Ejnzp3j77//Jjg42NGh2K5WLfjuO5g8WX380Udw111w7pxj46pAZAieEEIUz+oEwDQOv+C+ssjKymLfvn3s+197eGJiIvv27SMpKYmsrCxeeukltm3bxpkzZ0hISOCRRx6hcePG+aYjdklubjB1Knz7LdSsqTYNRERAOXewrKji4+OdbgieNCUIIZyF1aMAFEUhMjLSXBVx48YNnn32WapXr57vuFWrVllc5q5du7j77rvNj03zCAwfPpwFCxZw4MABPv30UzIyMggJCaFPnz5Mnz5d8+oQh3nwQdi1CwYOVIcI9u4N77wDEybIOgJlEBUVRVZWllMNwTMlADVq1JC57YUQDmV1AjB8+PB8j5966qkyB9GzZ88SRxb88MMPZX4Op9eokbqs8LPPqp0CJ05UVxhcvBicqf+CC5EheNpKSUlh7ty5DB48mPDwcEeH43JMTVLjx4+nXr16jg5HCOsTgPj4eHvEIUCdL+DTT6FrVxg7Fr7+Wh0muHKlOpugKJFer6dWrVrk5OSQlZXl6HAqnLi4OI4dO8aKFSvo3Lmzo8NxOXlnhZw6daqjwxHCNTsBVmg6HTz3nDpPQL16cOqUmhAsXAjlPP+Cq9Hr9fj5+eHn52dz9brRaCQrK0va6IsQHR1Ns2bNGDx4sKNDcUnR0dE0b97cqZqkROUmCYCzuv122LtX7R+Qna0mBUOGgCNnNKwEFEWRBKAYISEhjBkzhsDAQEeH4pKkSUo4G0kAnJmfH6xdC+++C+7u6loCEREye6AQQogykwTA2el08OKLapNAWBicPKnWDlSiJgEZz1928h4KIQqSBMBVdO2q3vnnbRIYNAguXXJ0ZKUq68VHltQtO3kPhRAF2bwaIEBCQgIJCQlcuHChUJvp4sWLyxSYKIKpSeD99+GVV9RRAtu3q8MGe/Z0dHTFKmvvZ0ctqWs0GjEYDM41jbSNZFliIURBNtcATJ06lT59+pCQkMDFixe5fPlyvk3YialJYNs2aNJEnTq4Vy91meFbtxwdXZHK2vvZUZ2nKtKsfdIBTQhRkM01AAsXLmTJkiUMHTpUy3iEpSIiYM8eGDcO4uLg7bfh559h6VJ1pUEnIhef8pW35kJmGxRCFMfmb4ebN29yxx13aBmLsFaNGvDJJ/DVV+riQjt2qBMGffpppekgKAqrSDUXQgj7sTkBGDlyJMuWLdMyFmEr0xoCd90FWVkQGanOGZCR4ejInJ70jhdCVFY2NwHcuHGDjz/+mJ9//pk2bdpQpUqVfD9///33yxycsEJYGCQkwMyZ6hLDK1bA1q1qk0D37o6Ozi60mFvd1EExPj6eyMhIbQMUQggnZnMNwIEDB2jXrh16vZ6DBw+yd+9e87ZPJqpxDDc3ePVV2LIFGjaEpCS1VuCVV9ShgxWMFkPbTB0Uo6KiNIxMCCGcn801AL/KmvXOq0sXdc6AMWNgyRK1VmDdOrVvQIcOjo5OM1oMbTN1UAwKCuL8+fMaRieEEM5NughXVN7eEB8Pq1dDnTpw8KCaGEyb5rTDBa1lzegC6RgnhBD5lSkByMjI4L333mPkyJGMHDmS999/n8zMTK1iE1ro31+9+A8cCDk5MGWKOqvg4cOOjqxclZYAmNa6v3DhQjlHJkoinTSFsB+bE4Bdu3bRqFEjZs+ezaVLl7h06RKzZ8+mUaNG7NmzR8sYRVkFBMDKlWqHwJo1YfdutSng3XchN9fR0TmF+Ph4jh07xsqVKx0dishDpjAWwn5sTgDGjx/Pww8/zJkzZ1i1ahWrVq0iMTGRBx98kHHjxmkYotCETgdPPAGHDkG/fmqnwJdeUjsJnjzp6OgcLioqimbNmnH33XezePHiEu84XbU5wVnjLimuss4iKYQoXplqAF5++WXc3f/pR+ju7s7EiRPZtWuXJsEJOwgJUTsELlqkTiS0ZQu0bQsffFCpawNMa93/+uuvnDp1ivj4+GKPLcuFNCUlhQ8++ICjR4+W+4VYURQyMzNdKgGQWSSFsB+bEwAfHx+SkpIK7U9OTsbb27tMQQk70+lg5Eg4cECtAbh2DcaPV+cLqGR9A/LS6XRERkbSuHFjuw0LNFVpL1q0yOkuxEKIysXmBGDw4MFER0ezYsUKkpOTSU5OZvny5YwcOZIhQ4ZoGaOwlwYN4JdfYMECddTAtm3qVMLTp8PNm46Ortzp9XqaN2/OyJEj7XbHGR0dTbNmzRg8eLBdyhdCCEvZnAC8++67DBgwgGHDhlG/fn3q169PZGQkjz32GDNnztQyRmFPej08+6zaN+CBB9QL/+TJ0KkTVMKmnNTUVOLi4uzW69zU1BAYGGiX8oUQwlI2JwAeHh7MmTOHy5cvs2/fPvbt22ceCeDp6alljKI8hIXBt9+qIwVq11abB7p0gYkT1SaCSmLJkiWcPHmyxD4AQghREZR5IqBq1arRunVrWrduTbVq1bSISTiKaaTAkSPqv0Yj/PvfaifBX35xdHTlwt59AIQQwllYNRVwTEwM06dPp3r16sTExJR4rCwG5MICAtSagCFD1OaBkyfhnnvgqafUuQOcvPpar9dTs2ZN9Hrr89vg4GCio6Ndute5FoskCSEqPqu+Iffu3cut/00jm3fxn4KbLAZUQTz4oNo3YPRotXbg88+heXNYuFCtHXBSZUkAKgKtJ88pOEzPWecTsFVFez1CWMqqGoC8CwB9+umnhIaGFvqSVRSF5ORkbaITjufrC3PnwrBham3A7t3w3HPqIkMLF0K7do6OsFzkvat29toBLRZJyst0gfTy8gL+mU+goiRZptdXo0aNCvF6hLCUzX/tDRo04OLFi4X2X7p0iQYNGpQpKOGEOnaE7dth3jx1yOD27RARATExcOWKo6OzOy3vqo1GI1lZWXa747TX5Dl6vR5fX190Op2m5QohHMPmBEBRlCL3Z2Vlme8URAXj5gYvvABHj8LgwWozwOzZ0KIFfP01FPM3URFoOSWtoihcuXKFS5cuuVS1c2VvWhGiorGqCQAwd/7T6XRMnjw5X8//3Nxctm/fTrtKUi1caYWEwPLlEBUFo0bBqVPw2GNw//1qDUHDho6OUHNa31UrioLBYMDf318uqEIIh7A6Adi7dy+gfoH9+eefeHh4mH/m4eFB27ZtmTBhgnYRCufVty/8+SfExsI778D69ZCQABMmwKRJUL26oyMUGjM1AxgMBkeH4nKkBkU4G6sTAFNHwKioKObMmYOPj4/mQQkXUrUqTJsGTz6pNg/8/DPMmKF2Evz3v+Hxxx0doc1cqeNfXrb2MTAajRgMBnx8fIq9SJkuYllZWVqE6lCm11te85eY3jshnIXNqWh8fDw+Pj4cPnyYDRs2sHbt2nybqGSaNYMff4TVq9U1Bs6fVycTuvNO+F+tkatx1bXoi0sA0tPTmTt3brHTHFe24XCV7fUKUZDVNQAmiYmJ9O/fnz///BOdTmfuFGjqIZxbiZeWrbR0OujfH+67D957D95+G37/HbcuXaj9+OPqvuBgR0dpMS2G0zlLta9er+ebb77h2LFjxMXFMXXqVIfGI4RwPJu/lcaMGUODBg24cOEC1apV49ChQ2zevJmOHTuyceNGDUMULsfLC157DY4dgyFD0CkK3l98gVvLluqcAv+bTMrZadHxz5kSgFGjRtGiRQvN5gfQiqmpxV4LMAkhimbzt9LWrVuZNm2auRezXq+ne/fuxMbGMmbMGC1jFK4qNBSWLSPn11/JbtkSXUYGjB2L0q6d2mGwAg8bdEb2mh+grFy1qUUIV2dzApCbm4u3tzcA/v7+5uw9PDycY8eOaROdqBi6dyd17VpuffghubVqoTt8WF16uE8fsOO00dbcfdt7ch5RPC3nWLCEtP0LobI5AWjVqhX79+8HoEuXLsyaNYstW7Ywbdo0GlbAceCijNzcyI2O5vzGjRhffBE8PNQRAx06qPMJnD+v+VNKAuAayrtmQhIAIVQ2JwCvv/66+QM0bdo0EhMT6dGjB+vXr2fu3LlWlbV582YeeughQkJC0Ol0rFmzJt/PFUVh8uTJBAcHU7VqVXr37s2JEydsDV04kNHHB+M776izCT7+uNoMsGQJNGkCkydDEcPLnKUdXSvp6enEx8eTnp7u6FCEEJWYzd+offv2ZcCAAQA0btyYo0ePcvHiRS5cuECvXr2sKuvq1au0bduW+fPnF/nzWbNmMXfuXBYuXMj27dupXr06ffv25caNG7aGLxytQQP44gvYtg26dYPr12H6dGjcGBYtgpwc86EVLQFYuXIlp06dYuXKlY4OpczkbloI16XpN6qfn59NC4X069ePt956i0cffbTQzxRF4YMPPuD111/nkUceoU2bNnz22WekpKQUqikQLqhLF/jtN3UtgcaNIT0dnnkGWrfWdH2B0nqa6/X6clsNbtCgQTRq1IhBgwbZ/bnsTRIAbcn7KcqTzd92vXr1KnIs8eXLl62uAShJYmIiaWlp9O7d27zP19eXLl26sHXr1mLPy87OxmAw5NuEk9LpYMAAOHQIPvgA/PzUJoLHHlMThJ9/LvNTlNbTvLwSAL1eT5MmTRgxYgSBgYF2fS6tlJQ8Sd8JbUkCIPIqeA3Lzs7WtHybv+02btzIhx9+SP/+/bl69ap5/82bN9m0aZMmwQGkpaUBFPqyDAwMNP+sKLGxsfj6+pq3sLAwzWISduLhAWPHwunT8MYb6loCO3fCvffCPffAjh02F13ePc2LY2rOcKUldUtKniQBEMJ+wsLC8l3HYmNjNS2/TLc7P//8M2lpadx+++2cOXNGo5C0MWnSJDIzM81bcnKyo0MSlvL1VdcXOH0axoxRE4NfflFrAwYMgMOHrS5Sy57mzj5xjda1Gc6SPAlR2SQnJ+e7jk2aNEnT8sv0DREcHMymTZto3bo1nTp1sssMgEFBQQCFekynp6ebf1YUT09PfHx88m3CxdSpA3PmwPHjEBkJer261kDr1urQwbNnNX9Kd3d3QkJCcHcvfpZsZ5+4RusEwFknEBKioit4DfP09NS0fJu/IUxVmJ6enixbtoyxY8dy33338dFHH2kWHECDBg0ICgoiISHBvM9gMLB9+3a6du2q6XMJJxUeDvHx6tLDjz4KRiMsWYJby5b4TZ4M585p9lTu7u7Uq1evxATAGe6Ipa1YCFFWNicASoHe2a+//jpLly7lvffes7qsrKws9u3bx77/zQqXmJjIvn37SEpKQqfTMW7cON566y3Wrl3Ln3/+ybBhwwgJCaF///62hi9cUcuWsGoVbN8OvXqhu3kTn//+F7dmzWD0aLtMJlQUZ5i4RhKAwirCeyJ9KkR5sjkBSExMxN/fP9++gQMHsm3bNhYvXmxVWbt27aJ9+/a0b98egJiYGNq3b8/kyZMBmDhxIqNHj+aZZ56hU6dOZGVlsWHDBry8vGwNX7iyzp0hIQF+/RXuvBPdzZvw4YfQqJHaZ8CKtnlXmGOgvC5szti3wdrZHCUBEMJyNn/rhYeHF/mhbNWqFcOHD7eqrJ49e6IoSqFtyZIlgNrcMG3aNNLS0rhx4wY///wzTZs2tTV0Uc7sdpHt2RM2blQ7CPboAdnZMG8eNGyojiaw4ELmqATg4sWLfPjhh051sXV034aiLuD2+P0UTHScMfERojwU39BZhJiYGKZPn0716tWJiYkp8dj333+/TIGJisP0JZ6TZ3Y/zeh0cPfdajLw668wZQr8/ru67PB//gP/+hdMnAhONuZ+7dq1nDt3jri4uCLn03CEqKgoLl68SFRUlF3KNxqNGAwGfHx8irygmxIAe8/HYEp04uPjiY6OJj4+3pz4OMvvQojyYFUCsHfvXm79by33vXv3FnucK41xFhWETge9eqnJwC+/qInAli1qIrBwIfqhQ3F/6il1ieIimO4Cx48fT7169ewaakpKCteuXaNu3brl2pGwtLvpoKAghg4dWuLomrIorwt8aaKjozEYDOZEJyoqiqysLBnmKCodqz6Fv/76KzVr1uTWrVvo9Xo+/vhjfv3110LbL7/8Yq94hSiZTqdOGvTbb/DTT2rTwM2b6OPiqHvPPeiHDoUDBwqdVp7V3/Hx8aSmplKjRo1yHVrnCv0dtFRcnwBTJ06AOXPmAMgwR1Ep2fRNUKVKFQ4U8SUqhNPQ6aB3b9i8GX77DWO/fuiMRvTLl0PbtvDQQ5BnKunyHNoXFRXlkmsBmC6o586dM7eZp6SkEBcX55Tt56V1CjRV/cfHx5dzZEI4B5tvBZ566imnnQhFiHy6d8e4di0p69ZhHDRITQ6++w7uuENtMvjpJ0KCg8vtLjAkJMSl1gIwMV1Q89aWxMfHc/LkSasvonkvzo6qmYiKiqJ58+Z26/MghLOzqg9AXjk5OSxevJiff/6ZiIgIqlevnu/n0glQOJubLVtiXLYM/VtvwaxZ8Nln6iiCjRtxi4igWmQkjBwJJUwC5Ars3Z8hb5t5Tk4OaWlpVl9E8/YHcHd3t18n0RKYmgKCgoJKXFdEiIrK5pT74MGDdOjQAW9vb44fP87evXvNm2lCHyGcUtOm8MkncOqUOlywalV0u3dTZ/Ro3Fq0UFckvHLFbk9vNBq5fv263cZ627s/Q96JkEJCQoiOjpb2cyFckM23Or/++quWcQhR/sLC1Iv9a69hnDsX5aOPcDtzBsaPhzffRP/007gNHJhv5IBpKFtZLt6mBKDgbJpaMfVyl17trsfUp2LSpEk0bNjQ0eGICq5ydAcWoiQBARinTOHcli3kzp+v1hBkZqJ/911Ce/TALSoK9u8H1It3Zmam3S7eWpDFe1yXrX0qhLBFmROAw4cPs2HDBtauXZtvE8LVKF5eKM88A0eOwNq1cNdd6HJy0H3+Oe4dOxL41FPof/gBnPjiX5lVhBn9oqKiaNy4sXRMFOXC5gTg9OnTtG3bllatWvHAAw/Qv39/+vfvz6OPPsqjjz6qZYyignCZL2i9Xh0muHEj7NwJjz+O4uZG1S1b8Ojfn7D77sNn6VLIyiq1qIowP31JvfSdaRigo6cy1oL0qRDlyeYEYOzYsTRo0IALFy5QrVo1Dh06xObNm+nYsSMbN27UMERRUbjkF3THjvDFFygnTnDj+edRvL3xOHWKgClTcKtfX+0vcPJksadX9ATAmipreyeAzrBMsxCuxOYEYOvWrUybNg1/f3/0ej16vZ7u3bsTGxvLmDFjtIxRVBCu/AWtb9AAr/nzyU1M5OIbb3AzPBxdZqbaibBpU3jgAdiwASy40Ov1eqpWrWrxlNnOXHNiTZV1aQlgWV+n9H0Qwjo2JwC5ubl4e3sD4O/vb/7QhoeHc+zYMW2iExVKhfiC9vUlc/hwkn/6idxvv4V+/dQ+AevXq/9v3lxdf8BgKLYIUwJg6cQ3zlxzYk2VdWkJYHm9zso2JbIQxbH5E9CqVSv2/69ndJcuXZg1axZbtmxh2rRpMnxFVHx6Pcp996kX/uPH1fkEfHzgxAn1/3XrwnPPgQZzYrhyzUlepSWAWr/O4moUJAEQQmXzJ+D11183t2tOmzaNxMREevTowfr165k7d65mAQrh9Jo0UZsCzp+H+fOhRQu1g+DChbh36kTwo4+iW7IErl2zqXhLak4qQl8DrWuInLnmRAhnYHMC0LdvXwYMGABA48aNOXr0KBcvXuTChQv06tVLswCFcBk1asDzz8OhQ5CQAP/3fyhVquC5bx9uTz8NISEwejS6Q4c0f+q8CYDRaCQrK8vqZECv15e4VK+WfRHKo19DRak5EcJeNK0D8/Pzs7hjkxAVlk4HvXrBihXkJiZyeeJElAYNIDMTPvwQz44daf3cc3ivWQPXr9v0FCVdQG1NANLS0oiLiyt2XnzTHXVZJ6kxGo3Mnz+fI0eO2PXuvEL0ORHCjqyaCjgmJsbiY2UxIOFK7NYuHBhI5nPP4T19Ou4bN8J//oPyzTf4/PknPi+/jPLOOzB0KIwYoS5TbKG81dtTp07VJNTSyjRNMZy3x39ptQZFMRqNPPLII9y4cUPuzguQqYBFebIqAdi7d69Fx0ktgHA1pgTAjk8AffpAnz5kJyaSHhtLyPr1VDl/Xh01MHcudOgAUVHwxBPg51dicfaY77+0MotaPc+WBAAgMDCQMWPGyN15AXnnVZg+fbqjwxEVnFUJgCwAJFyVU/X8Dg7m3PDhZI8fT6NTp3D77DNYswb27FG3F1+E/v3VWoHevYssoqjqbVM/gGrVqtkUlqtUmZsWZPLx8XGO36eGoqKibFpeWQhblOnTk5GRwXvvvcfIkSMZOXIks2fPJjMzU6vYhNCMlgmATqfD19fX5rJSUlKIj4/n+s2b6O6/H778ElJT1VqAdu3g5k113333Qf366CdPxv3MmVLLVRSFzMxMlx4JYImSRjy4+mgImQpYlCebvw137dpFo0aNmD17NpcuXeLSpUu8//77NGrUiD179mgZoxBOpazJRHx8PKdOneKrr776p4zatWH0aNi7V60FGD1abQY4dw59bCyhd9+N2913w6JFcPmyhq+mYqksSZAQWrA5ARg/fjwPP/wwZ86cYdWqVaxatYrExEQefPBBxo0bp2GIQlQspU6f2769WhuQkgJffomxb18UvR7d77/DM89AUBD6QYOotmEDZGeXb/BOTK/X4+vra3EfJKdqFqLwUE5XrskQrqFMNQAvv/wy7u7/dCNwd3dn4sSJ7Nq1S5PghHA2Woxft7ia19MTBg3C+N13nPv9d3LfeQfatIGbN9GvWUOd557DLTQUnnmGlNWrmTtnDunp6TbHVRQtL0T2vuBaW77peHd390LnOeICbDQauXLlCtWrVweQBEDYnc2fRB8fH5KSkgrtT05ONq8RIIQzs+VLXovZ5UrrOV9UXLnBwSgvvgj798P+/RgnTCAnKAhdRgYsWkT8//0fSYsWsXbsWHRHj9ocmyWx2MrZ7rhNiorLUXfgto6qEMIWNv+VDR48mOjoaFasWEFycjLJycksX76ckSNHMmTIEC1jFMIubPmSt3Z2uaJqDGxJAPJp0wZjbKxaK/DTTzBiBCNr1KDdjRuM2bWL6p060XrYMNzeew8SEy1+bUKIysXmBODdd99lwIABDBs2jPr161O/fn0iIyN57LHHmDlzppYxCuE0rB0qZ9f56N3cUHr2hLg4AlNS+Nfcufjecw+KuzvVT5ygyuuvQ8OG0KULvPceFFFjZ2/WJFla1xA4a42DEM7C5k+Gh4cHc+bM4fLly+zbt499+/Zx6dIlZs+ejaenp5YxCuGyym0++qpVuXr//aT95z9cPXmS0xMnktuzpzoB0Y4dMGEChIfDHXeoCxedO2ffeP6nqAQgJSWFuXPnFuqvIAmAEOXL6k+G0Whk5syZdOvWjU6dOjFt2jQaN25M69atbZ6ARIiKSovJdazu3BYQwLWhQ8n94Qd1JMFHH8Fdd6lrFGzdCuPHQ1gY9OgB8+apcxCUU2yg1oocO3aMFStW2Py8BUmveSGsZ3UCMGPGDF599VVq1KhB3bp1mTNnDqNGjbJHbEIIrL/I5lvUJzAQnnsONm5UlyueNw+6d1cP/P13GDMG6taFHj3QzZ6Ne3KyXWMDtVakWbNmDB482KrnKonWCUB5rFbozM8vKgerE4DPPvuMjz76iB9++IE1a9bw7bffsnTpUsm8hXAS8fHxRfc7CA6GF16A335TmwA++EBtElAU+P133CZOJPTOO3GLiICpU+HAAfVnBZS1at1UK9K0aVOby7D3Hb9d+25QetNIsb9DITRk9acvKSmJ+++/3/y4d+/e6HQ6yVRFheZKVcxRUVGl9zuoWxfGjoUtW9TOgfPmYbz7bhQ3N3QHDsCbb0Lbtrg1b06tGTNgyxb0YB437+Pjg8FgsPn9KEsSYbo7Pn78uN1+H/buu1HU39Mnn3zCoUOHWL58uWW/QyHKyOpPX05ODl5eXvn2ValShVu3bmkWlBDOxhUSgPT0dObOnQtgXb+DsDB44QWMP/5I8s6d5H7yCTz8MHh5oTt9Gt9PPsG9Z0/0oaHUfPll9D/+iPH6dYe9H/boQ1CQIxZGioyMJDw8nEGDBjlsYSZX+DsviavHX96sWg0Q1Lm2IyMj8/X0v3HjBs8++6x5BiuAVatWaROhEMIiK1as4NixY8THx9t852isVQtl+HCIjoarV8ldt47ry5ZRfeNGdOnp8PHH8PHHuFWvTkD37ugGDoSHHoKgII1fTfGio6PJyMgosg9BSkoK8+bNY/z48dSrV6/cYtJCSEgII0aMICAgwGExmC6grjoZkavHX96sTgCGDx9eaN9TTz2lSTBCVHT2ukPR6/U8/fTTfPLJJ9otJVu9OsqAAVzs3BmvOnVw//13WL0avvkGXWoq1X/4AX74QT22Uyd48EF44AF1LQM7fvmGhIQwZsyYIn+Wt+186tSpdouhLOQuVTgLqxOA+Ph4e8QhKgEZl61++dtjyWy9Xk/z5s0ZN24cQUFB6ggAG8oo9vfj4QF9+qjb/Pnk7NzJlRUrqPnbb+h27YKdO9VtyhS1s+EDD6gJQc+eZX9xVoiKiiIrK8up287t9TcghLUq7zexKHeSAGjPIXP16/UQEUHmmDHkbt2qziMQFwePPgrVq6uPP/kE+vfHLTCQwOHD0c2dC0ePFjmqQEuOajsXwhW5xDfxm2++iU6ny7c1b97c0WEJF6flWGtrLsQ6nQ5fX19NEiF7Vydb9B4FBcGIEbBqFfz9t9osMHo0NGiALjubqps34/bii9CiBdSvj/7ZZ6n2/feQkWFzXEajkaysLLtWo9szUTV12HS20VOuPv9ASkoKH3zwAUePHpUmFgu4RAIAcNttt5Gammrefv/9d0eHJFyclmO9tZzz3pm+hK1+jzw91WaCuXPh1ClyDhzg75df5uadd6J4eEBSEvq4OOo8/zxuQUHQrRtMmwbbt0NursVxKYpi9wTAnmWbOmwW9746qrbM3vMf2Jsp/kWLFkkCYAGXSQDc3d0JCgoyb/7+/o4OSbiI4i7O5TZPv5Wc6Uu4uPfIoiRFp+OctzevZ2ayY8YMci9cgHXrUEaPJrdJE3S5ufDHH2q/gdtvhzp1YPBgWLwYrJyR0NUMHjyYZs2a5XtfU1JSWLx4Menp6eYEACjXDoPO+pmwlD1mmazIXCYBOHHiBCEhITRs2JAnn3ySpFJWNsvOzsZgMOTbROVUXAJQUntxefTULktiYq87xIIxFfceWZqkxMXFceTIEb788ku1f8D996ObOxe348fhzBl1WOHAgeDrC5cuwZdfqkMQ69WDJk3g2WfVfX/9VeLzOFOtiSUCAwMZM2ZMvvd1yZIlnD17lm+//db8ey3vEQOu3ofCNEIkMDDQ0aFoouA1LDs7W9PyXSIB6NKlC0uWLGHDhg0sWLCAxMREevTowZUrV4o9JzY2Fl9fX/MWFhZWjhELV1fwi9ceFxhbEhOT8koAimPpnWJUVBSNGjVi0KBBhX8YHg5PPw1ffQUXL6qzEk6erC5frNfDyZPwn/+otQJ16kCbNjB+PLrvvkNf4LNvSkhceZTSyJEjad26NS+88IJ0lBUAhIWF5buOxcbGalq+1cMAHaFfv37m/7dp04YuXboQHh7Ol19+WewX0KRJk4iJiTE/NhgMkgQIm5kuMEuWLCEmJsbylfkq6MgHS+8ULZ7cxt1dXZfgjjvUdQgyM9U1CxIS4Jdf1HUJ/vwT/vwTtw8+oKpOxzp/f8Iee4x6AwYQ9fjjXLx4kcjISO1eZDlz9btvZ5Gens6KFSt48cUXXW4yqIKSk5Px8fExP847AZ8WXPJbqWbNmjRt2pSTJ08We4ynpyc+Pj75NiFsZbrjHTlypHVL81qQAFSEiWE0fw2+vuo8ArNnw/79cOGC2hTwr3+hNGnCJ4rCsb/+In7BArj3XsLbtSN2yxbCFi7Ea9MmKKF20BL2bFLIu+iP0F5pHSxdScFrmCQAQFZWFqdOnSI4ONjRoYhKQsu7s4IXS0kALBAQAIMGwcKF5B4+TJ+vvqJujx6MePRRqFsX3c2b+Ozfj8e77xIUGYmbvz907AgxMbBmjdrEYAV7dsQsj7UMbOXqNVamGTFduSNjeXKJJoAJEybw0EMPER4eTkpKClOmTMHNzY0hQ4Y4OjThRIxGIwaDAR8fn3L/ArNmXLqiKGRmZlKjRg27xOLqX+Kl0ev1NOzenceDgwno0AE8Pck+fJhzy5ZR7+xZ+O03qiQlwe7d6jZ7tnpiy5Zw553q1qMHhIYW+xzR0dEYDAa7XERKWsvA0fKOPnBFeWfElKaU0rlEAnDu3DmGDBnC33//TUBAAN27d2fbtm0OXTRDOB9HLgRiaQJgqv4dPHgwOp2OefPmMXr0aIviLS3BMX15u7u7u/SXeGn0en3+iZR0OpRGjfjroYcIbNWKixcvEgq4b90Kmzer2+HD/2wLF6rn1a+vzkPQtava96B1a7UvAvZtjzeV7YhE1REcmZiLkrnEb2P58uWkpKSQnZ3NuXPnWL58OY0aNXJ0WKICs1cbcN7qX9PCNZb2XC+qmj1vnBX9zt8qoaEwZAgsWACHDqnDCFevhvHjISJCHWVw5gwsXQovvAAdOkDNmtCrF7z2Grr169Ffvmy38Cz9Xbna8MaiVIQmropKvimEKIK1bcApKSnExcUV+UWd9wsw70QlUVFRNG/evEyr9xWM09m+bFNSUvjwww955513NL2IFfd+F3th9feH/v3h/fdh1y64fFmdsvjNN9WZC3184OpV+PVXePtt3B55hHodOuDWqhVERcGiRXDwIJTz++pMk0KJikcSAFHh2XIXZe2MaPHx8Zw8ebLIu/m8F+W8E5VoUc1cME5bE4DSzivq55Y8V3x8PFu3bmXz5s2aXsSKe78trgXx8VEv/FOmqInA5cvqMMOPP4bISJRmzQDQHTsGS5bAM8+oTQR+fnDvvfDqq2qNwvnz+YotSwJWVOyuPjOfI0hNmOVcog+AEGWR9y7K0jXirb04R0VFkZaWVqa7eVuUNYkwtc9Wq1aNjIwMvLy8ij0ub/8KvV6Pj48PBoOhxP4GUVFRnDx5kqpVq2p6ESv4fptqBCZNmkTDhg2tL1Cvh1at1O3pp8nNySHlzz8JOnOGawkJ+Bw6hH7nTnV+gp9/VjeT4GDo1Ak6dUKJiMAQHEyNVq2svgAV1QFP674IlaE93tU7MpYnSQBEhWfPHt0mISEhREdHO0XPY9PF2s/Pr9QvedOxpgt/SkoK8+bNY/z48SVOomL6ks3KyjKfN2fOnELnhYSE8MILLxAQEKDpe1Pw/c5bIzB9+vQyl6/X6/Fp0ABuu42M9u2pERqqVpcePAg7dsDOnep28KC6/PHatbB2LW5APUBp1EhNCjp3Vv9t316dCtlKWq966MiOssL5yF+AqPCKu4tytapCS5syFEXh8uXLXLp0yeoLh6ljorXV9Y5uq46KiqJx48aa1cAU+bfh7g7t2qnNAYsWwb59YDDA77+rQw2HDEFp3BgA3alTsHy5Og9Bjx5qk0OrVjB0KLrZs/Hats2i5ZDLY9lje6sIHRkrKtf45hPCAtZ+0WiZAJTUCbC0uCyNu7iLbFHtzqa5Bqy9cJg6JppqSyx9jxzdVh0aGsrYsWMJLWFsv11Uq6YOJRw3DpYtI/fIEZL27iV3/Xp46y145BG1icBoVEcjfP45bhMnEjRkCO4BAdCwITz2GMyYAevXq7UJFYyjk0NRPEkARIVhzy+a0i6EJXUCLBhXamoq7777LikpKdSsWZMlS5bk+3lxz1XcRTbvxb6s08wWrC2xNAEoqZalRo0a6HQ6m+KxlOl5nKE2x1izJsq998Jrr6mzEKakwLlz8O23MG0aPPqouhASQGIifP01vP46PPAAhISoCcP998Nrr6FfvRrPc+fKffSBlhydHIriOf7TIoRG7PlFU9qFsKQq6IJxffnllxw5coT4+Hhq1qzJyJEjLbrrtqRDWFxcHPv372fChAmazjVfVC2FJTUXznRhLg/Fvid166prG7zxBqxapc5B8Pff6mJH//43PPEEtGgBOh2kpcH338Pbb+PxxBO0HzQIzzp14Pbb1dUT582DTZvU5ZP/x9mGf+Ylixw5r8rxqRSVgpZfNNY2D5TUCTBvXEajkQceeICGDRuakwUt446OjubGjRt4eXlZNNd8WZofpGq3sE8++YSDBw+yaNGi0i/Ifn7qxEMTJqgTEh0+rC5i9McfMH8+REdzrlUr3tTrSb16FbZvh08+gTFjoGdPqF1bnfCoXz94+WVuxcdj3LcPbt4s9imlPV7kJaMAhCiCvYYSGY1GqlWrRmRkpEUX/LzDtiwREhLC+++/zxdffMGDDz5Y6h2h6SIeHx9fYs1JUSMpymN0RXHKq2ahuNENxYmMjCQ1NZVhw4bZ1tu+enV1auKuXQH4ePJkttapg7FpU6b17GleEpkDB9RahPPn4fx59Bs2YJ4Y3d0dmjWDNm3UuQtatlS3Bg1sGhIrKi5JAIRwYnmHbRUn7/rnISEhBAYG8vzzz1s0CsB0Eb///vuZM2dOsesSFFVLUbBmw5Le6raMQy/qnPJKAKy9YJpqgoKDg7lo5QqERRk5ciRXrlxh5PjxUK8e5F1AyGBQhyH++SfG/fu5uWsXnidOoMvIUDscHjoEX3zxz/EeHjzToAE5isLIli1h5Uo1MWjSBDw8yhyrcD2SAAjhIoprlsi7/vkbb7wBwIULF/jss88IDg4ucWIc00Xc0pqA4liTAFh7Z2zJOfYa0unIWg4opXnIx0ddxOiOOziXlMTs2bMZv3Il9fT6f2oJDh6EI0fU7fp16h47xgxQp0Q2cXNTk4AWLf6pLWjZUq1FqFq1vF6qcABJAESF50zj/fV6PdeuXWPFihWEhoZaNWtdcc0SgwcPZsWKFfkuUitXruTUqVMWT4wTFRVFVlZWsePobR2PrtfrqVWrlvmxvTqq2avJpuAF2NIajNTUVD766KNimw60npHPXFOxeLFaUxEWpo4k+OcJ4exZcg8eJHPrVmqmpKA/evSffgdHj6rb6tX/nKPTqaMVmjWDpk3/2Zo1U8u3MG5n+vyJ/CQBEBVeeUwNWlKVdMEvwHXr1nH69GlNZq3T6/U0bdq00F3i3XffzR9//MH9eS8CJTBd6IKCgkhLSyv089ISgOJev16vx8/PD4CcnByOHz9ubq6wpE0drG+HtydLazDyDu0squlA6xn5Sq2p0OuhQQOUsDAMrVvjExqK3t0dFEUdpph3uWTTdumS2s/gzBl1vYS8PD3VWoOCiUHTpmrnxDzDPmVqXuclCYCoMBx5p2FJAgDqF/+gQYO4evWqJrPWFZyS13THvXXrVnx8fFi/fj1d/9ehzJ7S0tKIi4sr9SKdt7mipDb1vHfIRbXDO+tdpWlCqOHDh3P16lUGDRqE0Wi0e5w2jyTR6dQhinXrqoscmSgKXLgAx4//sx07pv578iRkZ6vNCwcPFi6zVq38iUGjRv9sfn75kgPhWJIAiArDkXca1qzvvnLlSh5++GGLvqytvfs13XE///zzxMbGajo1bsEEJ+9rtrSzXFHNFUXJe4dc1N2ts95VLlmyhDNnzrBhwwZGjx4NUGQC4Gy1GoWaI3Q6CAxUtx498p+QkwNJSfmTAtOWlKSurLh9u7oV5OsLjRqhNGhAdmgoHi1aoG/SRE0OQkPV/gii3EgCIIQGLL0gxcXFcezYMa5evUrfvn0tOr60C2tRyYelixMVPLe4RKakBAAs6yyXkpLCihUrGDx4sFV3qlrMk6BVm3tpF25Tr/0RI0agKEqx5ViaMFmaWJZ2XN5VH4v6WWnNEfneP3d3dQrjhg3hvvvyH3jtGpw69U9icOKE+vjUKbWpITMT9uxBt2cPhdadrFIFGjRQy81ba9CoEdSvb9NiSqJkkgAIUY6io6O5fPkyDz30kEUXIksurGW5Gzadm5OTU6ayLJ2l8NixY6xYsYLOnTvbFK+ttGpzL+3Cnfd9MBgMGAyGImMZNGgQmZmZpdaEWPr7MB1nep0FE51z584xe/Zsxo4dWyhRKC2pMRqNXLp0CYPBUPr7V62aOvdA69aFf3btmjr18alT5B4/ztU//6TGhQvoT59W99+69U9NQlH8/dVEIDw8/7+m/1s4V4b4hyQAQmispLtN0wUCKPWuO+/xjp5GVYsq6+joaDIyMhicdyy7Bc/n6Neel6XDAgv2zcjLaDRStWpVRo8erflrM13oC/6e8g7zLJi4lJbUGI1GMjMz8z22qTalWjW47Ta47TbOJyUx+/x5xk+frsaZm6uul2CqLTh9+p//nzql1hxcvKhuu3YVXX6tWoWTgrz/1qwp/Q8KkARACI1Ze7fprO3ZeWkxg5wpmbHkwpH3+aZMmVLmDn/WJDAFL3B5n1urhEyn02k2BDCv4n5PJSUuliQ1BSebKmttSqE43dzUC3V4uDo9ckEZGXD2rLqdOfPPv6b///232vfg8mV1meai+PiowxdDQ4v/19vbptfjqiQBEKKc2evLPy9rZ8orrQ25tIuEVm3VRT2fFgmSNQlM3gTO3d3d5ucuKeko74mLQkNDmTJlSpFTSpeW1KSkpBATE4OXl1e+yabsEWexatZUt7Zti/55Vlbh5CBvknDhgjpzommGxOL4+pacIISFQQmzcroaSQCEKGflMYTN3d2d0NBQ3N0t+4iXdpEt7SJhyUXamtetddOHpRccLVfVc8S8+8Vd6Ev6/ZT2e4mLi8PLy4sbN24QHR1tTmxM/QlsSWY1b9qqUcPcvFCka9fUEQrnzkFy8j//5v1/ZuY/W0lJQs2aajIQEvLPFhyc/9+gIHWuBCcnCYAQFZAzNis4MiZLLzgF27utlfdiWlzSYc8hgLa8x3q9Hh8fn2Lb9fP23QgJCWH69OnmxCY6Oto1lnuuVg2aN1e34ly5UjhBKPivwaA2R2RkFD0HQl61axdOEAomC8HBDl2HQRIAITSm5Re8s054U5Hlbe+2tuo+7wW4uKTDGVfkK6nfSkhICGPGjDE/LriAlC1/53mTDqf5+/b2VtdDaNGi+GMMhn+SgdRUdWij6d+8/791S+2X8Pff6roMJfH3/ycZCApSt3JKlCUBEEJjWn7BO+OdfHnQOvGxZqKmvO3dBX9/Wvw+HL3AUFFKS1rz9lspuICULX/npgTg3Llzdu8Poykfn38WSyqOoqgX/pIShLyJgml0Q2mJgh1IAiCExkr6gnfVO/ryjlvrxMeaiZrytnfbg7MM7cyrpIu5aXbJvBfqmjVrEh0dTVZWllMlMk5Bp1Pv6v39i54PwaRgopCW9s+WnAxff233UCUBEEJjJX3Bu+odvavGba2C7d2VhSVJa8HHNWrUKFMi40zTITtESYmCwVAuCYBr3YYIIYQdmdq7AwMDHR1KubKlVqKstUJ5ax2EY0gNgBBC2IGrNvdYqqy1Qs7YF6KykQRACCH+x7Scsun/ZS2rqAukMyYGjojJGftCVDaSAAihMWf8gheWMXV4s/dzOFt/CmeMSdifJABCaEy+TIUonSTKjicJgBBCiHInibLjSeolhBBCVEKSAAghhBCVkCQAQgghRCUkCYAQQghRCUkCIIQQQlRCLpUAzJ8/n/r16+Pl5UWXLl3YsWOHo0MSQgghXJLLJAArVqwgJiaGKVOmsGfPHtq2bUvfvn25cOGCo0MTQgghXI5OURTF0UFYokuXLnTq1IkPP/wQAKPRSFhYGKNHj+aVV14pdHx2djbZ2dnmx5mZmdSrV4/k5GR8fHzKLW4hhBDCGgaDgbCwMJKSkvD19TXv9/T0xNPTU7snUlxAdna24ubmpqxevTrf/mHDhikPP/xwkedMmTJFAWSTTTbZZJOtQmxTpkzR9NrqEjMBXrx4kdzc3EJLdAYGBnL06NEiz5k0aRIxMTHmx0ajkYiICPbs2YNOp7NrvFrp1KkTO3furLDPr3X5WpRXljJsOdfSc0x3BFKDZTlHf36s5eh45fNu/3MtPcdUY3327Nl8syVqevdPBZ4KuKiqEk9Pz3zVKc7Ozc3NoV/29n5+rcvXoryylGHLudae4+PjIwmAhRz9+bGWo+OVz7v9z7X2nJo1a9r1d+ISnQD9/f1xc3MjPT093/709HSCgoIsLmfUqFFah2ZXjo7X3s+vdflalFeWMmw519G/44rM1d5bR8crn3f7n+vo33FBLtUJsHPnzsybNw9Qq/Tr1avHCy+8UGQnQCEqEoPBgK+vL5mZmS51VyuEsF55fd5dpgkgJiaG4cOH07FjRzp37swHH3zA1atXiYqKcnRoQtidp6cnU6ZM0bwNUAjhfMrr8+4yNQAAH374If/+979JS0ujXbt2zJ07ly5dujg6LCGEEMLluFQCIIQQQghtuEQnQCGEEEJoSxIAIYQQohKSBEAIIYSohCQBEEIIISohSQCEcDHJycn07NmTli1b0qZNG1auXFnifiGE6yruc52RkUHHjh1p164drVq1YtGiRVaXLaMAhHAxqamppKen065dO9LS0oiIiOD48eMYDIYi91evXt3RIQshbFTc593Ly4vs7GyqVavG1atXadWqFbt27aJ27doWl+0yEwEJIVTBwcEEBwcDEBQUhL+/P5cuXSIsLKzI/ZIACOG6Svq8V6tWDYDs7GwURcHa+3lpAhCinG3evJmHHnqIkJAQdDoda9asKXTM/PnzqV+/Pl5eXnTp0oUdO3YUWdbu3bvJzc0lLCzMov1CiPJlz897RkYGbdu2JTQ0lJdeegl/f3+rYpMEQIhydvXqVdq2bcv8+fOL/PmKFSuIiYlhypQp7Nmzh7Zt29K3b18uXLiQ77hLly4xbNgwPv74Y4v2CyHKnz0/7zVr1mT//v0kJiaybNmyQgvmlUoRQjgMoKxevTrfvs6dOyujRo0yP87NzVVCQkKU2NhY874bN24oPXr0UD777LN85xa3XwjheFp/3vN67rnnlJUrV1oVj9QACOFEbt68ye7du+ndu7d5n16vp3fv3mzduhUARVGIjIykV69eDB061HxccfuFEM6pLJ/39PR0rly5AkBmZiabN2+mWbNmVj2/dAIUwolcvHiR3NxcAgMD8+0PDAzk6NGjAGzZsoUVK1bQpk0bc3vif//7XzIzM4vc37p16/J8CUIIC5Xl8379+nWeeeYZc+e/0aNHW/1ZlwRACBfTvXt3jEZjkT8rbr8QwjWV9Hnft29fmcqWJgAhnIi/vz9ubm6FOvOkp6cTFBTkoKiEEPbg6M+7JABCOBEPDw8iIiJISEgw7zMajSQkJNC1a1cHRiaE0JqjP+/SBCBEOcvKyuLkyZPmx4mJiezbtw8/Pz/q1atHTEwMw4cPp2PHjnTu3JkPPviAq1evEhUV5cCohRC2cOrPu1VjBoQQZfbrr78qQKFt+PDh5mPmzZun1KtXT/Hw8FA6d+6sbNu2zXEBCyFs5syfd1kLQAghhKiEpA+AEEIIUQlJAiCEEEJUQpIACCGEEJWQJABCCCFEJSQJgBBCCFEJSQIghBBCVEKSAAghhBCVkCQAQgghRCUkCYAQQghRCUkCIIQQQlRCkgAIIYQQlZAkAEII4QSSk5Pp2bMnLVu2pE2bNqxcudLRIYkKThYDEkIIJ5Camkp6ejrt2rUjLS2NiIgIjh8/TvXq1R0dmqigpAZAiCL07NmTcePGOToMp1dR3qe///6bOnXqcObMGYfFEBwcTLt27QAICgrC39+fS5cumX/++OOP89577zkoOlERSQIgKpy0tDRGjx5Nw4YN8fT0JCwsjIceeoiEhARHh6aJyMhIdDodOp2OKlWqEBgYyL333svixYsxGo3lGsuqVauYPn26+bGWCcFdd92FTqfjiy++yLd/3rx5hISEaPIcJjNmzOCRRx6hfv36mpZrq927d5Obm0tYWJh53+uvv86MGTPIzMx0YGSiIpEEQFQoZ86cISIigl9++YV///vf/Pnnn2zYsIG7776bUaNGOTo8zdx3332kpqZy5swZvv/+e+6++27Gjh3Lgw8+SE5OTrnF4efnh7e3t+blKorC3r17CQ4O5uuvv873s927d9OhQwfNnuvatWvExcURHR2tWZnFadeuHa1atSq0paSkmI+5dOkSw4YN4+OPP853bqtWrWjUqBGff/653eMUlYQiRAXSr18/pW7dukpWVlahn12+fFlRFEW5ceOGMnr0aCUgIEDx9PRUunXrpuzYsSPfsXfddZcyduxY8+Pw8HBl9uzZ+Y5p27atMmXKFPPxL7zwgjJ27FilZs2aSp06dZSPP/5YycrKUiIjI5UaNWoojRo1UtavX1/oeUaPHq289NJLSq1atZTAwEBzmcUZPny48sgjjxTan5CQoADKokWLzPtyc3OVt99+W6lfv77i5eWltGnTRlm5cqVVz79y5UqlVatWipeXl+Ln56fcc8895vc37/s0fPhwBci3JSYmKp9++qni5+en3LhxI1+5jzzyiPLUU08V+RqPHTumAMr8+fOVatWqKVevXjX/rHXr1srkyZNLfI+ssXLlSiUgIKDQ/tzcXGXmzJlKo0aNFA8PDyUsLEx56623FEWx/fddmhs3big9evRQPvvssyJ/PnXqVKV79+7Wv0ghiiA1AKLCuHTpEhs2bGDUqFFFdpyqWbMmABMnTuTrr7/m008/Zc+ePTRu3Ji+ffvma2+1xaeffoq/vz87duxg9OjRPPfccwwaNIg77riDPXv20KdPH4YOHcq1a9cKnVe9enW2b9/OrFmzmDZtGj/99JPVz9+rVy/atm3LqlWrzPtiY2P57LPPWLhwIYcOHWL8+PE89dRTbNq0yaLnT01NZciQIYwYMYIjR46wceNGBgwYgFJE3+E5c+bQtWtXnn76aVJTU0lNTSUsLIxBgwaRm5vL2rVrzcdeuHCBdevWMWLEiCJfy+7du/Hy8mLkyJH4+Pjw/fffA3Djxg2OHDmiaQ3Ab7/9RkRERKH9kyZN4p133uGNN97g8OHDLFu2jMDAQPPPbf19F0dRFCIjI+nVqxdDhw4t8pjOnTuzY8cOsrOzbXuxQuTl6AxECK1s375dAZRVq1YVe0xWVpZSpUoVZenSpeZ9N2/eVEJCQpRZs2aZ99lSA5D3ziwnJ0epXr26MnToUPO+1NRUBVC2bt2a73kK3tF16tRJefnll4t9DcXVACiKogwePFhp0aKFoijq3WS1atWUP/74I98x0dHRypAhQyx6/t27dyuAcubMmSKfr+D7VPCxyXPPPaf069fP/Pi9995TGjZsqBiNxiLLnTBhgtK5c2fzuY8//riiKIqybds2BVCSkpKKPM8WjzzyiDJixIh8+wwGg+Lp6ZmvNiUvW3/fJfntt98UnU6ntG3b1rwdOHAg3zH79+8v8fchhDXcHZl8CKElxYIRradOneLWrVt069bNvK9KlSp07tyZI0eOlOn527RpY/6/m5sbtWvXpnXr1uZ9prvHCxcuFHseqL3BCx5jKUVR0Ol0AJw8eZJr165x77335jvm5s2btG/f3qLnb9u2Lffccw+tW7emb9++9OnTh8cee4xatWpZFdfTTz9Np06dOH/+PHXr1mXJkiXmzoxF2bNnj/kuf8CAAQwYMIDs7Gz27NlDQEBAvs5x77//PpcuXeKtt94CICcnh9q1a/P333/j7l76V9z169fx8vLKt+/IkSNkZ2dzzz33FHuerb/v4nTv3r3UTpxVq1YFsLhWQYiSSAIgKowmTZqg0+k4evSo5mXr9fpCCcatW7fyPa5SpUq+x6Ze+nkfA4W+5Is6z9be/EeOHKFBgwYAZGVlAbBu3Trq1q2b7zhPT0+Lnt/NzY2ffvqJP/74gx9//JF58+bx2muvsX37dvPzWKJ9+/a0bduWzz77jD59+nDo0CHWrVtX7PF79uxhyJAhgDqyoEqVKvzwww9FdgA8ePAg999/f773oH79+hZd/AH8/f25fPlyvn2mC21JbP19l4WpmSogIECzMkXlJX0ARIXh5+dH3759mT9/PlevXi3084yMDBo1aoSHhwdbtmwx77916xY7d+6kZcuWxZYdEBBAamqq+bHBYCAxMVHbF1BGv/zyC3/++ScDBw4EoGXLlnh6epKUlETjxo3zbXnvoEuj0+no1q0bU6dOZe/evXh4eLB69eoij/Xw8CA3N7fIn40cOZIlS5YQHx9P7969i43h9OnTZGRkmC/07u7uPPzww3z99dfs2bOnUHv9wYMHadu2rfnxvn378t2Jl6Z9+/YcPnw4374mTZpQtWpVpxs6evDgQUJDQ/H393d0KKICkBoAUaHMnz+fbt260blzZ6ZNm0abNm3Iycnhp59+YsGCBRw5coTnnnuOl156CT8/P+rVq8esWbO4du1aicPAevXqxZIlS3jooYeoWbMmkydPxs3NrRxfWX7Z2dmkpaWRm5tLeno6GzZsIDY2lgcffJBhw4YB4O3tzYQJExg/fjxGo5Hu3buTmZnJli1b8PHxYfjw4aU+z/bt20lISKBPnz7UqVOH7du389dff9GiRYsij69fvz7bt2/nzJkz1KhRAz8/P/R69T7jiSeeYMKECSxatIjPPvus2OfcvXs3Hh4etGrVyrxv4MCB5g51r732mnm/oigcPXo0XzPH5cuXmTRpUqmvzaRv375MmjSJy5cvm5s2vLy8ePnll5k4cSIeHh5069aNv/76i0OHDpXLcMHi/Pbbb/Tp08dhzy8qFkkARIXSsGFD9uzZw4wZM3jxxRdJTU0lICCAiIgIFixYAMA777yD0Whk6NChXLlyhY4dO/LDDz+U2K49adIkEhMTefDBB/H19WX69OkOrQHYsGEDwcHBuLu7U6tWLdq2bcvcuXMZPny4+YILMH36dAICAoiNjeX06dPUrFmTDh068Oqrr1r0PD4+PmzevJkPPvgAg8FAeHg47733Hv369Svy+AkTJjB8+HBatmzJ9evXSUxMNE+u4+vry8CBA1m3bh39+/cv9jn37NlDq1at8PDwMO+79957yc3N5ebNm/maABITE2natCm7du0y7+vXrx+tW7dm5cqV/PTTT1y+fJkpU6bkSyjyat26NR06dODLL7/kX//6l3n/G2+8gbu7O5MnTyYlJYXg4GCeffZZS942u7hx4wZr1qxhw4YNDotBVCyyFoAQotzcc8893HbbbcydO1eT8tauXcvKlSv573//a94XHh7O5s2bCQ8PB9QmgR9//JGJEycWW866det46aWXOHjwYL4EypksWLCA1atX8+OPPzo6FFFBSA2AEMLuLl++zMaNG9m4cSMfffSRZuUePHgwX9+Nq1evkpmZab74G41GPvroIyZPnlxiOQ888AAnTpzg/PnzVvWPKE9VqlRh3rx5jg5DVCCSAAgh7K59+/ZcvnyZmTNn0qxZM83KPXjwIP/3f/9nfnz06FFzQmA0Ghk/fjzPPfccoaGhpZbl7IsajRw50tEhiApGmgCEEBXSrFmzWLNmDW3atOHee+81j44QQqgkARBCCCEqIefs7SKEEEIIu5IEQAghhKiEJAEQQgghKiFJAIQQQohKSBIAIYQQohKSBEAIIYSohCQBEEIIISohSQCEEEKISkgSACGEEKISkgRACCGEqIQkARBCCCEqIUkAhBBCiEro/wFqSvvWO5tzSAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the column density\n",
    "PvN_plot_Step2 = HtoH.PvN(CatP, CatdP, CatN, Nscale=[2e+21,10e+22], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "c49996f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfMAAAF2CAYAAAB6cnEgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABPdUlEQVR4nO3deXgT1f4/8HdCN5YmUAq0pQsFkcVCKVgRUUBAFmXHe5HLUhBcy9Yq0iqCgBr21V746QWKyqJ4AXdQKhQXKNCGpYishRbZZGlDWUrbzO8Pv8lt2rRJJpNl0vfrefJAJpOZz0nS+cw5c+YchSAIAoiIiEi2lK4OgIiIiOzDZE5ERCRzTOZEREQyx2ROREQkc0zmREREMsdkTkREJHNM5kRERDLHZE5ERCRzTOZEREQyV22SeVFREd555x0UFRW5OhSH8fQyenr5AJbRU3h6GT29fID8yqhwh+FcV65ciZUrV+LcuXMAgIceeggzZsxA3759AQD37t3Da6+9hk2bNqGoqAi9e/fGv//9bzRq1Mjqfeh0OqjVahQUFEClUjmiGC7n6WX09PIBLKOn8PQyenr5APmV0S1q5qGhoZg7dy4yMzNx8OBBdO/eHQMHDsSxY8cAAAkJCfj666+xefNmpKen4+LFixgyZIiLoyYiInIPXq4OAAD69+9v8vy9997DypUrsW/fPoSGhmL16tXYsGEDunfvDgBYu3YtWrVqhX379uHRRx91RchERERuwy2SeVmlpaXYvHkzbt++jU6dOiEzMxPFxcXo2bOncZ2WLVsiPDwce/furTSZFxUVmVzryM/PBwAUFBQ4NH5X0ul0Jv96Gk8vH8AyegpPL6Onlw9wfhkFQcD169cREBAApfJ/jea+vr7w9fW1agNu4ciRI0Lt2rWFGjVqCGq1Wvj2228FQRCE9evXCz4+PhXWj42NFd54441Ktzdz5kwBAB988MEHH3zI9jFz5kyrcqjb1MxbtGiBQ4cOoaCgAF988QXi4uKQnp4uenvJyclITEw0Pi8oKEB4eDjy8vJk0ZmBiMwrKSnBn3/+icaNG8PLy6vS5SUlJcjNzQUAhIeHm6xL5G50Oh3CwsKQm5sLtVptXG5VrRxu1Mzu4+ODBx54AADQoUMHHDhwAMuWLcOwYcNw//595Ofno27dusb1r1y5gqCgoEq3V1nThEqlYjInkrGSkhLodDqoVKoKybzs8pKSEvj7+wNAhXWJ3JVarRaVo9yiN7s5er0eRUVF6NChA7y9vZGWlmZ87cSJE8jNzUWnTp1cGCEREZF7cItT1eTkZPTt2xfh4eG4desWNmzYgN27d2PHjh1Qq9UYN24cEhMTERAQAJVKhYkTJ6JTp07syU5ERAQ3SeZXr17F6NGjcenSJajVarRt2xY7duzAU089BQBYsmQJlEolhg4dajJoDBEREbnJCHDOILfRfIjIvJKSEly4cAGhoaEVrpmXXV5SUoLz588DACIiInjNnNyavTnKba+ZExERkXWYzImIiGSOyZyIiEjmmMyJiIhkjsmciIhI5pjMiYiIZI7JnIiISOaYzImIiGSOyZyIiEjmmMyJiIhkjsmciIhI5pjMiYiIZI7JnIg8VnZ2NkaMGAGtVuvqUIgcismciDyWRqNBTk4Opk6d6upQiByKyZyIPFZycjIiIyOxYMECV4dC5FCc4JeIPFZUVBTWr1+PiIgIV4dC5FCsmRMREckckzkREVUqNTUVCoUCCoXCuKxbt25QKBQYM2aM6wIjE0zmRORQer0e+fn50Ov1brEd+h9DUjb32LZtGwCgQYMG6NixIzp27OjaYCU0ZswYKBQKdOvWzer3bNu2DV26dIG/vz9q1qyJ5s2bY+7cuY4L0kZM5kTkUK5K5lqtFiNGjEB2drZd+60OfHx8jAnb8AgICAAAPPPMM9i3bx/27dvn4ihdZ9GiRRg8eDB+/vln1KlTB61atcLdu3eRlpbm6tCMmMyJSFa0Wi2GDx9u8d7xqVOnIicnBxqNxkmRyVdwcLAxYRseXbp0AWC+md2cuLg4NG/eHP7+/vDx8UFERAQmTZoEnU5nXMdQI27SpAnWrVuHsLAwqFQqTJkyBXfv3sWUKVOgVqsRHh6OVatWmWz/4sWLeP755xESEgIfHx80bdoUc+bMQUlJiXEdQ0vD6NGjMXPmTAQHB6NevXoYOXIkbt26BQDGfQNAenq6sWy7d+82W668vDwkJSUBAJYvX46LFy8iKysLFy5cwJYtW2z7oB2IyZyIZCUpKcmqe8cXLFiAyMhIJCcnOymy6u3LL7/EzZs30axZM4SFhSE3NxcrVqzAuHHjKqx76dIlvPrqq/D19cWtW7ewbNkydOjQAWvXroVarUZeXh7i4+Nx/PhxAMD169fx6KOPYu3atSgsLESrVq2Ql5eHGTNm4MUXX6yw/U2bNmHJkiWoWbMm8vPzsX79emOTeExMDAIDAwEA/v7+xpYIlUpltlxbtmxBSUkJateujX379iEwMBDBwcEYNWoUbt++LdXHZzcmcyKSlblz51p173hMTAzWr1+PqKgoJ0UmX+fPn69wzdxW6enpuHbtGg4dOoQzZ87grbfeAvD3teZ79+6ZrHv//n388MMPOHnyJMLCwgAAp0+fxuHDh/H777/Dz88Per0e6enpAIAPPvgAeXl5aNSoEc6cOYPDhw/jiy++APB3y8Hp06dNtu/n54fjx4/j9OnT6NChAwAYm8S3bt2KZ555BgDQvn17Y0tE+/btzZbrxIkTAIDbt29j8+bNCA4OxvXr1/Hpp5/i6aefRnFxsc2flSMwmRORrMTExGDjxo2IiYlxdSgew9w1c1vt3LkTUVFRqFmzJhQKBd577z0AQElJCf766y+TdevVq4fOnTtDqVQiPDwcwN9jAjRp0gR16tRBw4YNAQBXrlwBAOzfv9/4vGHDhlAoFBg0aBAAQBAEZGRkmGy/e/fuaNy4MZRKJVq2bGmyLVuVbcZfs2YNsrOzsXr1agB/X/L59ddfRW1Xahw0hoiomjNcMxdr/fr1eP31143bCgsLw7Vr13D27FkAQGlpqcn6ZZu0vby8KiwztAwIgmDyPn9/f7Ru3brC/mvVqmXyvG7duhW2X35b1mrcuLHx/7GxsQCARx55xLjs3LlzorYrNdbMiYjILoYTAX9/f+Tk5CAjIwO9evWSbPuGJOrl5YVNmzYZm8Z//PFHvPrqqxg8eLBN2zMkf2uueffs2dP4/4MHD5r8CwDNmze3ad+OwmRORER2adu2LQDg1q1baNq0KZo2bYrPP/9csu3Hx8ejcePGuHnzJlq0aIF27dqhWbNmqF+/PuLi4mzenqHp/eDBg2jTpg0effRR3L171+y6nTt3xsCBAwEAY8eORZs2bTB27FgAQI8ePdC5c2eRpZIWkzkREdll3LhxSExMRGBgIG7duoVu3bph9uzZkm2/QYMG2LdvH8aOHYv69evj2LFjuHv3Lp544gksWbLE5u09//zzGDp0KNRqNbKzs5GRkVHhUkBZmzZtwrRp0xAUFIRTp04hMjISb7/9Nr7++mt7iiUphSD2QoLM6HQ6qNVqFBQUVHoLAhFJr6SkBBcuXEBoaKjx+qU92wkKCsLly5crbK/8fkpKSnD+/HkAQEREhF37JnI0e3MUa+ZE5BGsHUyGyBMxmRORRzCM+GZpMBkiT8RkTkQewTDim6XBZIg8ES8iEZFHMAwmExoa6upQiJyONXMiIiKZYzInIiKSOSZzIiIimWMyJyIikjkmcyIiIpljMiciIpI5t0jmGo0GsbGx8Pf3R8OGDTFo0CDjhPAG3bp1g0KhMHm8/PLLLoqYiIjIfbhFMk9PT0d8fLxxSrvi4mL06tWrwvR0L7zwAi5dumR8zJ8/30URExERuQ+3GDRm+/btJs9TU1PRsGFDZGZmokuXLsbltWrVQlBQkLPDIyIicmtuUTMvr6CgAAAQEBBgsnz9+vUIDAxEVFQUkpOTcefOnUq3UVRUBJ1OZ/IgIiJyZ+XzVlFRkVXvc7tkrtfrMWXKFHTu3BlRUVHG5f/617/w6aefYteuXUhOTsYnn3yCkSNHVrodjUYDtVptfISFhTkjfCIiItHCwsJMcpdGo7HqfW7RzF5WfHw8srOz8csvv5gsf/HFF43/b9OmDYKDg9GjRw+cOXMGzZo1q7Cd5ORkJCYmGp/rdDomdCIicmt5eXkm85n7+vpa9T63SuYTJkzAN998gz179licLKFjx44AgNOnT5tN5r6+vlZ/CERERO5ApVKZJHNruUUyFwQBEydOxNatW7F7925ERkZafM+hQ4cAAMHBwQ6OjoiIyL25RTKPj4/Hhg0b8OWXX8Lf3x+XL18GAKjVatSsWRNnzpzBhg0b8PTTT6N+/fo4cuQIEhIS0KVLF7Rt29bF0RMREbmWWyTzlStXAvh7YJiy1q5dizFjxsDHxwc7d+7E0qVLcfv2bYSFhWHo0KGYPn26C6IlIiJyL26RzAVBqPL1sLAwpKenOykaIiIieXG7W9OIiIjINkzmREREMsdkTkREJHNM5kRERDJnVwe44uJiXL58GXfu3EGDBg0qjKVOREREjmdzzfzWrVtYuXIlunbtCpVKhSZNmqBVq1Zo0KABIiIi8MILL+DAgQOOiJWIiIjMsCmZL168GE2aNMHatWvRs2dPbNu2DYcOHcLJkyexd+9ezJw5EyUlJejVqxf69OmDU6dOOSpuIiIi+j82NbMfOHAAe/bswUMPPWT29UceeQTPP/88Vq1ahbVr1+Lnn39G8+bNJQmUiIiIzLMpmW/cuNGq9Xx9ffHyyy+LCoiIiIhsI8kIcNeuXUNGRgZKS0sRGxvLyU+IiIicyO5k/t///hfjxo3Dgw8+iOLiYpw4cQIpKSkYO3asFPERERGRBTb3Zi8sLDR5PmvWLOzfvx/79++HVqvF5s2b8dZbb0kWIBEREVXN5mTeoUMHfPnll8bnXl5euHr1qvH5lStX4OPjI010REQOpNfrkZ+fD71e7+pQiOxiczP7jh07EB8fj9TUVKSkpGDZsmUYNmwYSktLUVJSAqVSidTUVAeESkQkLUMyr1OnDpRKDohJ8mVzMm/SpAm+/fZbbNy4EV27dsWkSZNw+vRpnD59GqWlpWjZsiX8/PwcESsRERGZIfpUdPjw4Thw4AAOHz6Mbt26Qa/Xo127dkzkRERETiYqmX/33XdYtGgRDh48iP/85z+YP38+RowYgalTp+Lu3btSx0hEMqbVajF8+HBotVpXh0LksWxO5q+99hrGjh2LAwcO4KWXXsKcOXPQtWtXZGVlwc/PDzExMfj+++8dESsRydDUqVORk5ODqVOnujoUIo9lczJPTU3Fd999h02bNuHAgQP45JNPAAA+Pj6YM2cOtmzZgvfff1/yQIlInhYsWIDIyEgsWLDA1aEQeSybk3nt2rWRk5MDAMjLy6twjbx169b4+eefpYmOiGQvJiYGGzduRExMjKtDIfJYNidzjUaD0aNHIyQkBF27dsWcOXMcERcRERFZyeZb00aMGIE+ffrg7NmzaN68OerWreuAsIjIVfR6PXQ6HVQqFe+9JpIJm/5Sc3NzAQD169dHbGysxUT+559/ig6MiFyDo6IRyY9NyTw2NhYvvfQSDhw4UOk6BQUF+OijjxAVFYX//ve/dgdIREREVbOpmf3333/He++9h6eeegp+fn7o0KEDQkJC4Ofnh5s3b+L333/HsWPH0L59e8yfPx9PP/20o+ImIqqSUqmEWq2GTqdzdShEDmdTzbx+/fpYvHgxLl26hA8++ADNmzfHtWvXcOrUKQB/X0/PzMzE3r17mciJyKWUSiXq1q3L6/5ULYiaz7xmzZp49tln8eyzz0odDxEREdmIp6xEREQyx2ROREQkc0zmREREMsdkTkREJHOiOsCdP38eR44cQaNGjfDII49IHRMRkSTYo52qC5uT+caNGzFmzBgUFxdDoVAYpzxt0KCBI+IjIhLNkMyJPJ3Np6uzZs3Cv/71L/zxxx/44YcfAABJSUmSB0ZERETWsTmZnz17FjNnzsSDDz6IHj164NNPP8WmTZscERsRkZFWq8Xw4cOh1WpdHQqR27E5mZeUlKBWrVrG5y1btoRer8fly5clDYyIqKypU6ciJyeHLYFEZojqFbJu3Tr89ttvKCwsBAB4eXnhzp07kgZGRFTWggULEBkZifnz57NTG1E5NneAe+KJJ/Duu+/i1q1bUCqViIyMxL1797B69Wr07NkTDz/8MPz9/R0RKxFVYzExMdi4cSNCQ0Ph5SXqRhwij2XzX0R6ejoA4NSpU8jMzERWVhaysrKwcuVKaDQaKJVKNG/eHMePH5c8WCIiIqpI9Olt8+bN0bx5czz33HPGZTk5OTh48CA7qBCRS+n1euh0OqhUKjbHU7Ug6a88MjIS//jHP/D+++/b9D6NRoPY2Fj4+/ujYcOGGDRoEE6cOGGyzr179xAfH4/69eujTp06GDp0KK5cuSJl+ETkIfR6PfLz86HX610dCpFTuMUpa3p6OuLj47Fv3z78+OOPKC4uRq9evXD79m3jOgkJCfj666+xefNmpKen4+LFixgyZIgLoyYiInIPbtGLZPv27SbPU1NT0bBhQ2RmZqJLly4oKCjA6tWrsWHDBnTv3h0AsHbtWrRq1Qr79u3Do48+WmGbRUVFKCoqMj7X6XSOLQQREZGdyucqX19f+Pr6WnyfW9TMyysoKAAABAQEAAAyMzNRXFyMnj17Gtdp2bIlwsPDsXfvXrPb0Gg0UKvVxkdYWJjjAyciIrJDWFiYSe7SaDRWvc8tauZl6fV6TJkyBZ07d0ZUVBQA4PLly/Dx8akwxnKjRo0qHawmOTkZiYmJxuc6nY4JnYiI3FpeXh5UKpXxuTW1csDOZJ6Wloa0tDRcvXq1QkeTNWvWiNpmfHw8srOz8csvv9gTmtVNE0RERO5CpVKZJHNriW5mnzVrFnr16oW0tDRcu3YNN2/eNHmIMWHCBHzzzTfYtWsXQkNDjcuDgoJw//595Ofnm6x/5coVBAUFiS0CEf0f9v4mkjfRNfNVq1YhNTUVo0aNsjsIQRAwceJEbN26Fbt370ZkZKTJ6x06dIC3tzfS0tIwdOhQAMCJEyeQm5uLTp062b1/ourOkMzr1Knj6lCISATRyfz+/ft47LHHJAkiPj4eGzZswJdffgl/f3/jdXC1Wo2aNWtCrVZj3LhxSExMREBAAFQqFSZOnIhOnTqZ7clORERUnYhuZh8/fjw2bNggSRArV65EQUEBunXrhuDgYOPjs88+M66zZMkS9OvXD0OHDkWXLl0QFBSELVu2SLJ/IiIiORNdM7937x4+/PBD7Ny5E23btoW3t7fJ64sXL7Z6W4IgWFzHz88PKSkpSElJsTlWIvJ8SqXSOJsar/1TdSM6mR85cgTt2rUDAGRnZ5u8plAo7AqKiMhWhmQOgMmcqh3RyXzXrl1SxkFEREQiueUIcERERGQ9uwaNyc/Px+rVq41zl7du3Rrjxo2DWq2WJDgiIiKyTHTN/ODBg2jWrBmWLFmCGzdu4MaNG1iyZAmaNWuGrKwsKWMkIiKiKoiumSckJGDAgAH46KOP4OX192ZKSkowfvx4TJkyBXv27JEsSCIiIqqc6GR+8OBBk0QOAF5eXnjjjTfw8MMPSxIcERERWSa6mV2lUiE3N7fC8ry8PPj7+9sVFBFVbxwrnsg2opP5sGHDMG7cOHz22WfIy8tDXl4eNm3ahPHjx2P48OFSxkhE1QyTOZFtRDezL1y4EAqFAqNHj0ZJSQkAwNvbG6+88grmzp0rWYBERERUNdHJ3MfHB8uWLYNGo8GZM2cAAM2aNUOtWrUkC46IiIgss+s+cwCoVasW2rRpI0UsRFRN6fV66HQ6qFQqKJUcy4rIVjYl88TERMyZMwe1a9dGYmJilevaMtEKEVVvZedTZzInsp1NyVyr1aK4uNj4/8pwohUiqi7YqkDuwKZkXnZylXXr1iE0NLTCj1cQBOTl5UkTHRGRm2OrArkD0b+8yMhIXLt2rcLyGzduIDIy0q6giMhzlJ1nnIgcQ3QHOEEQzC4vLCyEn5+f6ICIyLOUnWeciBzD5mRu6PimUCgwY8YMk1vRSktLkZGRgXbt2kkWIBGRo2i1WkyaNAnLly9HbGysq8MhEs3mdi+tVgutVgtBEHD06FHjc61Wiz/++APR0dFITU11QKhErseRyeRBq9Vi+PDhVXbUBYCpU6ciJycHU6dOdVJkRI5hc83c0Alu7NixWLZsGVQqleRBEbkrdnaSh7JJevfu3ZWut2DBAkyaNAkLFixwXnBEDiD6aLR27VomciJySwsWLEBkZKTFJB0TE4ONGzciJibGSZEROYboZK7RaLBmzZoKy9esWYN58+bZFRQRkT2YpKm6EZ3M/9//+39o2bJlheUPPfQQVq1aZVdQRNUJr8O7H34nJDeik/nly5cRHBxcYXmDBg1w6dIlu4Iiqk6YONyPO34n7hgTuQ/RyTwsLAy//vprheW//vorQkJC7AqKiIhMMZlTVUQn8xdeeAFTpkzB2rVrcf78eZw/fx5r1qxBQkICXnjhBSljJCInsva2LiJyH6JHgJs6dSquX7+OV199Fffv34cgCKhZsyamTZuGpKQkKWMkIiey9rYuT5hgRK/Xo7CwkLVdkj3Rf4EKhQLz5s3DX3/9hX379uHIkSO4ceMGZsyYwVnTiGTM2tu6PKHZl8mcPIXomrlBbm4url+/jvv37+Ps2bPG5QMGDLB300TkAobbukJDQ10dChFZSXQyP3v2LAYPHoyjR49CoVAYJ14x1MpLS0uliZCIiIiqJLqZffLkyYiMjMTVq1dRq1YtZGdnY8+ePXj44YervM5GRERE0hJdM9+7dy9++uknBAYGQqlUokaNGnj88ceh0WgwadIk9oQlIiJyEtE189LSUvj7+wMAAgMDcfHiRQBAREQETpw4IU10REREZJHomnlUVBQOHz6MyMhIdOzYEfPnz4ePjw8+/PBDNG3aVMoYiYiIqAqik/n06dNx+/ZtAMDs2bPRr18/PPHEE6hfvz4+++wzyQIkIiKiqolO5r179zb+/4EHHsAff/yBGzduoF69erzPnIiIyIlEXTMvLi5Gjx49cOrUKZPlAQEBTOTk0TjICBG5I1HJ3NvbG0eOHJE6FiK3x2RORO5IdG/2kSNHYvXq1VLGQkQuwslViORNdDIvKSnBypUr8fDDD+Oll15CYmKiycMWe/bsQf/+/RESEgKFQoFt27aZvD5mzBgoFAqTR58+fcSGTkTllJ1cxRV4MkFkH9Ed4LKzs9G+fXsAwMmTJ01es/W6+e3btxEdHY3nn38eQ4YMMbtOnz59sHbtWuNzX19fGyMmososWLAAkyZNsji5iqNYO1MbEZknOpnv2rVLsiD69u2Lvn37VrmOr68vgoKCJNsnEf2PqydXsXQy4QkztBE5ks3N7GfPnjVOquJMu3fvRsOGDdGiRQu88soruH79epXrFxUVQafTmTyIyD0ZTiZiYmLMvq7X61FQUFDpsYfJnjxF+bxVVFRk1ftsTubNmzfHX3/9ZXw+bNgwXLlyxdbN2KRPnz74+OOPkZaWhnnz5iE9PR19+/atcmY2jUYDtVptfISFhTk0RiJyDGsSNZM5eYqwsDCT3KXRaKx6n83N7OXPjL/77jurdybWc889Z/x/mzZt0LZtWzRr1gy7d+9Gjx49zL4nOTnZpCOeTqdjQieSIUu1ciJPkpeXB5VKZXxubf8w0dfMXalp06YIDAzE6dOnK03mvr6+7CRHRESyolKpTJK5tWxuZjfcGlZ+mTNduHAB169fR3BwsFP3S0TOwVvViGwjqpl9zJgxxlrvvXv38PLLL6N27dom623ZssXqbRYWFuL06dPG5zk5OTh06BACAgIQEBCAWbNmYejQoQgKCsKZM2fwxhtv4IEHHjAZH56IXE+v10On00GlUkGpFD2Mhcmtajt37pQwQiLPZHMyj4uLM3k+cuRIu4M4ePAgnnzySeNzw7XuuLg4rFy5EkeOHMG6deuQn5+PkJAQ9OrVC3PmzGEzOpGbMXREq1Onjl3J3NX3vRPJjc3JvOzALVLp1q1blZ1bduzYIfk+icTQarWYMGECPvjgA3Tq1EmybU6aNAnLly9HbGysJNuUO1ff904kN+JPnYmqoaSkJOTm5iIpKUmybbp6KFW5UCgUUKvVVtX4lUol6tata1frAJGc8JdOZIO5c+ciPDwcc+fOlWybCxYsQGRkJJuULTAkaAAW7ylnMqfqRpa3phG5SkxMDFJSUtCuXTtJt8kmZetJdV2eyJPwL4GIiEjmmMyJSBRPvhfck8tGnsmuZva0tDSkpaXh6tWrFa5frVmzxq7AiMi9efK0pbaUjXcjkDsQXTOfNWsWevXqhbS0NFy7dg03b940eRCRZ3Onjnt6vR6FhYWSTbRiS9l4NwK5A9E181WrViE1NRWjRo2SMh4ikgl36rgndTK3pWwc4Ibcgehkfv/+fTz22GNSxkJEDibVcKv0P+50UkPVl+i/5vHjx2PDhg1SxkJEDsZ5v4k8k+ia+b179/Dhhx9i586daNu2Lby9vU1eX7x4sd3BERHR39jRjqoiOpkfOXLEOHBGdna2yWvOnhKViEgMR4y17yiefPcA2U90Mt+1a5eUcRAROV3ZsfbT09NdHU6V2NGOqsIeMERUbTlirH1HMXS0i4mJcXUo5IbsGjQmPz8fq1evxvHjxwEArVu3xrhx46BWqyUJjoiqD0PnvICAAKft0xFj7RO5guia+cGDB9GsWTMsWbIEN27cwI0bN7BkyRI0a9YMWVlZUsZIRNWAIAgoKCgw29M+OzsbI0aM4PCqRJUQncwTEhIwYMAAnDt3Dlu2bMGWLVuQk5ODfv36YcqUKRKGSOQ+tFot4uPjmVScTKPRcJQ1oirYVTOfNm0avLz+11Lv5eWFN954AwcPHpQkOCJ3U7bDlKeSejQ1KSQnJ1c5vCpPsqi6E53MVSoVcnNzKyzPy8uDv7+/XUERuauqOkx5yoAsjkrm9nw+UVFRWL9+faWdvyo7yfKU74TIEtHJfNiwYRg3bhw+++wz5OXlIS8vD5s2bcL48eMxfPhwKWMkchuGDlPmkgoTR9VTh1b1+SiVSqjVatFjVFR2ksXvhKoL0b3ZFy5cCIVCgdGjR6OkpAQA4O3tjVdeeUUWt3kQkfTEDmyiVCpRt25dFBYWitove6VTdSe6Zu7j44Nly5bh5s2bOHToEA4dOmTs0e7r6ytljEQkE+40LSpRdWL3oDG1atVCmzZt0KZNG9SqVUuKmIhIppwxsElVTflE1ZVNzeyJiYmYM2cOateujcTExCrX5UQrROQIHKOcqCKbkrlWq0VxcbHx/5XhRCtE5Cgco5yoIpuSednJVdatW4fQ0FAolaYt9YIgIC8vT5roiEhScppG09ATXaVSmSw3NOWHhobavQ+lUok6depUOI4RyY3oX3BkZCSuXbtWYfmNGzcQGRlpV1BEnqrsrVKuuG2qbBO1uzMM75qZmYkRI0ZUmGpZCkzm5ClE/4IFQTC7vLCwEH5+fqIDIvJkrk7mcuxtnpSUhJycHGg0GleHQuS2bL7P3NDxTaFQYMaMGSY92EtLS5GRkcF7PYlsYLjH2hm1QymbqAHTmc4cFf/cuXORkJCA5ORkh2yfyBPYnMwNHd8EQcDRo0fh4+NjfM3HxwfR0dF4/fXXpYuQyMMZkrkcGZrCbTkZKXst3Jr3xMTEYP369faGSuTRbE7mhk5wY8eOxbJlyyp0TiEiqoohmfNaNZF0RA/nunbtWgDA77//jtzcXNy/f9/k9QEDBtgXGZEbqg4dpqpDGYk8jehknpOTg0GDBuHo0aNQKBTGDnGGe8xLS0uliZDIAfR6PXQ6ndVNvQZyT3TWzIgmxzJqtVpMmDABH3zwATp16mR2HbHfOZEciP5FT5o0CZGRkbh69Spq1aqFY8eOYc+ePXj44Yc5KhO5veo6m5Y1yVyOw6VaM898df3OqXoQncz37t2L2bNnIzAwEEqlEkqlEo8//jg0Gg0mTZokZYxE5GBle9Rbey+6vdOWSqmqeeaJqgPRyby0tBT+/v4AgMDAQFy8eBEAEBERgRMnTkgTHZGMyLFGa1A2mVt7L7ozb6mzpKp55omqA9F/hVFRUTh8+DAAoGPHjpg/fz5+/fVXzJ49G02bNpUsQCK5kNPoalVxxsxnZB4vBZBYopP59OnTjT+42bNnIycnB0888QS+++47LF++XLIAiZzF3gOpHEdXc1dSN+HLJUnKJU5yP6J7s/fu3dv4/wceeAB//PEHbty4gXr16rnFNTQiW9l7/7PUo6tVZ4Ym/MLCQkm2p9frUVBQIMm2iNyRpBe7AgICRCXyPXv2oH///ggJCYFCocC2bdtMXhcEATNmzEBwcDBq1qyJnj174tSpUxJFTUTOIKc+Be7UH4DIGqJ/qd27d8esWbMqLL958ya6d+9u07Zu376N6OhopKSkmH19/vz5WL58OVatWoWMjAzUrl0bvXv3xr1790TFTlQduFuTraFPQVW3j7kLJnOSG9G/1N27d+ODDz7AoEGDcPv2bePy+/fvIz093aZt9e3bF++++y4GDx5c4TVBELB06VJMnz4dAwcORNu2bfHxxx/j4sWLFWrwJC/ulmw8jbt9voY+BZXdPianmjuRu7HrtHPnzp24fPkyHn30UZw7d06ikEzl5OTg8uXL6Nmzp3GZWq1Gx44dsXfv3krfV1RUBJ1OZ/Ig9+JuyYYcy1Iv+fJ3A5RN7gqFAmq1mjVl8njl81ZRUZFV77PrLyM4OBjp6elo06YNYmNjHTLy2+XLlwEAjRo1MlneqFEj42vmaDQaqNVq4yMsLEzy2Ig8gZQnVfY0T5e/G6Bsszybvam6CAsLM8ldGo3GqveJ/sswdHTz9fXFhg0bMHnyZPTp0wf//ve/xW5SUsnJySgoKDA+8vLyXB0SkcuZG3fdXZJ5+Zq7pWZ5a7DpnuQmLy/PJHclJydb9T7Rt6YZJlYxmD59Olq1aoW4uDixmzQrKCgIAHDlyhUEBwcbl1+5cgXt2rWr9H2+vr7w9fWVNBYiuZPTJCqG5B4UFGTSCmfuhKGycpWt3a9Zs8ZpsROJpVKpRE0tLvovOicnB4GBgSbLhg4din379kn6RxMZGYmgoCCkpaUZl+l0OmRkZFQ6OxIR/Y+n9U2wJZlLUbsnkgPRNfOIiAizy6OiohAVFWXTtgoLC3H69Gnj85ycHBw6dAgBAQEIDw/HlClT8O6776J58+aIjIzE22+/jZCQEAwaNEhs+EROU3bqTVft3zAYTnVTtnb/559/Vnid1+LJU9iUzBMTEzFnzhzUrl0biYmJVa67ePFiq7d78OBBPPnkkyb7AYC4uDikpqbijTfewO3bt/Hiiy8iPz8fjz/+OLZv3w4/Pz9bwicSzZCQ69SpY/PB352SqVarxaRJk7B8+XLExsZaXN/T5wA3JHMiubMpmWu1WhQXFxv/XxlbR4Hr1q1bhWvw5bc3e/ZszJ4926btEkmlbEKW6uBva2KVQtnbv6y5+8TeIW7N0Wq1SExMdGq5PQFbEagqNiXzXbt2AQCKi4uhVCqxatUqNG/e3CGBETmbo5OroZZruHZta2I1ty1ba8wLFizApEmTXDoZTFJSkuhyuyNnJVm2IlBVRP36vL29ceTIEaljIXIpR09hapjsw9AKZWmWtao6ront1GbN9Ka2JCcxt37NnTvXo2aXY42Z3IHoX9/IkSOxevVqKWMhchpzScjZU5haSqyu6oVuS3IScwLkivnStVotRowYgezsbKftk8iZRPdmLykpwZo1a7Bz50506NABtWvXNnndlg5wRM5mrolb7BSmhibvWrVqWVw3OzsbGo0GK1askDSZuaqjmjs021vD0LSv0WjwzDPPuDocIsmJ/qvPzs5G+/bt4e/vj5MnT0Kr1Rofhw4dkjBEIulJWQs31KABWKzRajQahzTlm6vFO2P0M1fUssUwNO1bO5oWkdyIrpkbOsMRyZHYWnhVrOmglJycDI1GU+VJhC01/aqUHf1s3bp1dm1L7mJiYoyffUhICHvRk8dhjw0iJ4qKisL69eurrMlKda3cE0Y/k7LfgKNaRYjcgeiaucHvv/+O3Nxc3L9/32T5gAED7N00eThX3GddnVQ2trmcCIIAnU6HGjVq2L0ta1pFiORKdDI/e/YsBg8ejKNHj0KhUBhvtzEMGFNaWipNhOSx7LnPmsjA2tq7oVWksqGoieRMdDP75MmTERkZiatXr6JWrVo4duwY9uzZg4cffpgHZrKKs28FI8/k6UPOEllD9C9/7969mD17NgIDA6FUKqFUKvH4449Do9Fg0qRJUsZIHkouPaHt5e5zart7fNbgwC1U3Yn+5ZeWlsLf3x8AEBgYiIsXLwL4eza1EydOSBMdkQco26vcHTl65Dtr6fV6FBYWGpvLlUol1Gq1yb+ezhNOrMg1RP91REVF4fDhwwCAjh07Yv78+fj1118xe/ZsNG3aVLIAieRObK9yZx3Y3eVyh7lkXrduXXh5eVWbWre7nFiR/Ij+65g+fbrxj2727NnIycnBE088ge+++w7Lly+XLEAiuTNcTgBg05CiZSckcaSqLndIeULBpnDL3OXEiuRH9F9V7969MWTIEADAAw88gD/++APXrl3D1atX0b17d8kCJHIWscnG2oRXdkhRa7jDhCRS1hSZzC2rLv1ISHqS/lUFBATYPJc5kbsQm2ysTXi2DinqDgd2OdYUXTVBDZEr2XSfeWJiotXrcqIVqi6snWwkJiYG69evd1JU0nDEsLeOZkjmderUcXUoRE5jUzK39roZa+dUnUid8Awj4y1evBjBwcE2vcfSaHrVqanbcP85a+hUHdiUzDm5CpHjlZ8gxVDTLDsoSvnkXdVoelqtFomJicZ1LU0G4yn0ej0KCgpcHQaRU9h1ep6fn49FixZh/PjxGD9+PJYsWcI/HvJIzrz/t/ytbOauAZe/Tm/u2rahFu6sXvHuzHCvOlsNyVOJTuYHDx5Es2bNsGTJEty4cQM3btzA4sWL0axZM2RlZUkZI5HLOfP+X2s6vpVP3ubeY0jmCxculF0nNqlVp8sLVD2J/mUnJCRgwIABOHfuHLZs2YItW7YgJycH/fr1w5QpUyQMkcj1pOjVXX5QFHvY0tPdHXrFW0OpVKJOnTpMuEQi2FUznzZtGry8/nfZ3cvLC2+88QYOHjwoSXBE7sLWhGiuadyQzA0zDMqF1Ld6VVZLNpfM7b28oVAoqs1QsFS9if6Fq1Qq5ObmVliel5dnHLOdqLrypHudBUFAQUGBw5O5OfZe3mDzOlUXon/hw4YNw7hx4/DZZ58hLy8PeXl52LRpE8aPH4/hw4dLGSNRBfYmS2ce5JlQxJPjoDVEriD66LJw4UIMGTIEo0ePRpMmTdCkSROMGTMGzz77LObNmydljOSh7ElyVSVzaxK9UqmESqVyyn3IVTUrM8lXzZrLG+Wb4jnzGFVHoo8iPj4+WLZsGW7evIlDhw7h0KFDuHHjBpYsWQJfX18pYyQP5ahkZm2t3dVN4ZbKb+/n44zyWROjo+Mo3xQv9ZSzrv6dEFnD5qOEXq/HvHnz0LlzZ8TGxmL27Nl44IEH0KZNG9SqVcsRMRJVS9UtmYutUZdvihc75WxlmMxJDmw+Srz33nt48803UadOHTRu3BjLli1DfHy8I2IjkgVrk66ht7ZCoWCCMENsZ7fyTfFyuRWPSEo2J/OPP/4Y//73v7Fjxw5s27YNX3/9NdavX8+DElVbtiZzpVJpUzKvLteA2dmNSDybk3lubi6efvpp4/OePXtCoVDg4sWLkgZGJDeOqm07evS5qk5GnDkMqqNr1OxwSJ7M5l91SUkJ/Pz8TJZ5e3ujuLhYsqCI7GFtTVbqGq+jkrmja6yWkrlcE2D571fOZSGyxKZZ04C/B5AYM2aMSY/1e/fu4eWXX0bt2rWNy7Zs2SJNhEQ2qmoGMTHrlWeYWrPsLGZSKp905DinuDsoO8GMLd8vkRzZfCSKi4tDw4YNoVarjY+RI0ciJCTEZBmRI1VVq7ZUkzXUoOfNmyeqxmtPDdya4UXF1CA9pdYpZevG3Llz3fYaPDtAktRsrpmvXbvWEXGQG7O3JuqImmxVtWpLNVnDgTQ6Otq4nqNr24DpHORSXxc2JHNHKD93uiMZho6V4sTEnVs0DL9BTixDUuGviCyytxbhiFqI1NeRnVFTcuY0qlKSa9xE1QmTOcmSuZ7P7tR0ae4ygFxvvaosbnf6vMvylEsORLbgr508hjslF3O1WbkOZlJZ3O70eZcl52Qu59jJtfiLIXIAV95OZg93TdDVBZM5icVfDFlkaDLOzMzkgd5KhtpsdHS0Qz4zJnMiKksWyfydd96BQqEwebRs2dLVYVUbhibjN954Q9SB3tnDkVpKdIZ4Dh8+LCoh6vV6FBYWWvU5uDI5GsqZlZVldbyOjEPM9+8ONVVn/X55IkX2kEUyB4CHHnoIly5dMj5++eUXV4dUbdg7C5Wze0NbSgCGeKZNm+bwZF4ZRyYIQ1J4/fXXjVOBujKZ2/P9u0Myd9bvl8mc7CGbZO7l5YWgoCDjIzAwsMr1i4qKoNPpTB4kjr0dt1zRi7uqA6O5eCwlV2uSii0J2pYEYWtCM5TdMGjK+++/b9X7HMWW7z87OxsjRoxwq0ll3PEuBCZ+z1U+bxUVFVn1Ptkk81OnTiEkJARNmzbFiBEjkJubW+X6Go3GZES6sLAwJ0VK5TmiF7elBFfVwc5cPJaSqzUJ1ZYEbUuCEFs7NZSzffv2Ng1OInWisOX712g0bndPu72/X0ckXkvbZLKXr7CwMJPcpdForHqfLJJ5x44dkZqaiu3bt2PlypXIycnBE088gVu3blX6nuTkZBQUFBgfeXl5ToyYHE3q5lcpal+2bMOZt6mVnXrVGoZR2FyRCJKTk23+HtyhKb4q5hKro6/DM5nLV15enknuSk5Otup97vnrL6dv3774xz/+gbZt26J379747rvvkJ+fj88//7zS9/j6+kKlUpk8yD7uPq+2PfHZkly1Wi3i4+Mr7MeRCbo6HJyVSiU6d+6MDRs22PQZunsyNyj7HXJUPapM+bxVdlKzqrj3r78SdevWxYMPPojTp0+7OpRqwXCwLDsLlTty1gEyPj4eR48eRXx8vEP3U5Yzk7kUHfzEkEtSFqvsd+iO1+FJ3mT5V1NYWIgzZ84gODjY1aFUC4aD7MKFC0UdgJx1kOYBUhqCILi097tc2HOCFRMTg/Xr1yMyMpKfM0lCFsn89ddfR3p6Os6dO4fffvsNgwcPRo0aNTB8+HBXh1atiG1GdlYyd9Z16JSUFLRp0wYpKSk2vc/Ta57uxtGft72tJa7sm2BOdbiU48lkcVS5cOEChg8fjhYtWuCf//wn6tevj3379qFBgwauDo0scMcDhL0H+ZiYGKSkpLjtSY3UKovb3ftQuMvnrdVq8dxzz2HPnj1u9XdQnjv+rZL1ZHFU2bRpEy5evIiioiJcuHABmzZtQrNmzVwdFlnBVQeIqhKNuxzk3ZFWq8WoUaNw8uRJ47LKPi924rKO4XN6++23IQiCU/bp7idaJD0ezcgjuTrRuPKEofyB3NbBbM6dO4cVK1ZYXFeOfRRccXJp+JySkpKq3L+UCdjVv39yPiZzspqcarS2Jhqpy2bY3uHDh41NrCUlJU5JJIYDeVJSEgDYdBfCggUL0KRJE0ycONHiurZ04rImiToj0boimRs+p4ceeqjK6+S2JGBLiV+OJ1pkH/c/KpPbkEMyNxyso6OjbeoM56iyGQ7Qb775pl3J3JZaW/mx9A3Dulo7mM0nn3yCBx980Kq4rO3EZc16er0eOp0OKpXKru9Brtd+bUnAlhK/MwclIvfgvkdlIhGk7iFsLjHYMqKa4QBt7ShOlbGl1lb+QG7rgf3rr79Gnz59sHHjRrtiFkOKkypHJ3NHbd+W74k1byqPyZwcSu4dcexN5oYm1qioKLvicObBe9asWbh9+zYSEhIcvi+pOSLRlj/BELOP7OxsjBo1CseOHZMkJta8qTwmc3IoZ3bEUSqVUKvVUCgUDt9XZft31GUIMQdvsfHMnDkTtWvXxpIlS2wN0+X0ej0KCgok7TUuxfeq0Whw6tQpjBw5EoMGDYJWq3W7GeLkfuJd3TGZk0M5s0bp6mv6tuzfGbGK3Uf//v2xfft2yQZlcvVJVllSJCzDcLdZWVlWbUupVOL9999HSUkJ7t+/jxMnTiApKQkLFiyw+UTXkZcQ2ANe3pjMyaHYHGieq088nEmpVCInJ8diLdQZNUMpEpYhmVt7l4BSqUSXLl3wySefoEWLFmjRogUWLlyIJUuWoGnTpjad6DoymfM6vLx5/pGEyMVsrcl5ImuSqFQ1w7LN1+WTnxQJyzBr3pgxY2zaVlRUFLZs2YJt27YhJibGoSe6Yk4WeeItb0zm5FEMtbvDhw9LUvM1V1u05UCp1+tx4MABjB07FvHx8U5vxnSX27SsSaJS1Qw1Go3xcy5ffikSVlJSEnJzc5GamlphW/Z83lLOVledWn7ob/ymyaMYanfTpk2T5GBmrrZoazJfsGAB8vLyAMDpzZiO6BBWVlW14LKsGWDG1kRb2f6Sk5Md+jnPnTsX4eHheP/9962OyRrWJHN2UqPKMJmTR5H6up8U2+vXrx/u3buHiRMnelwzZlW14PKcMQYA8Hdz9vr16x32ORsm2mnXrp3V75GqE2D50f2IDJjMyaGc3dwn9XU/Kbb3zTffwM/PD6mpqaK3YfgcAVis+ZX9zLVaLUaMGIHz589L/h0Yemk7ohZsSNS2jppnKG92drak8ZTfR3x8PA4dOmT1e6T6OzCcXM6fP5/N6GSCvwRyKF67A+Lj4xEaGmq2WdZaYpO5YeKU2bNn23TLnL+/v8X1Db20N23aJHktWGwyN9RcNRqNw/oLGK6Zv/nmmza/V6FQQK1W2zX97saNG9GhQ4dq/3dFpvhLIHKwBx98EMuXL7epWVYqYiacCQgIEDU+uqXruVIPkmJuf2WHzzXXX0CKk8uqrplb+gyq2r8tIws6Ak+85Y3fGpGdqqoBHj58GImJiTh16pQLIrNtZjN7TZ06FWfPnkVCQkKFfSmVSouDpNiaTMx1Tiw7fK65JncpElaHDh2wbt06tG/f3qqYKlM+lsOHD+OFF17A4cOHRcdmDyZzeeO3Rh7F3tt7xDTNVvWeadOmIS8vDytXroSXl5dLDpb2djyz9jMx1IoNneHKUiqVFgdJKX+t31Kv7cpaHQxN2UlJSTh58iRGjRolae/vqmrQtrSElE+eHIGN7MFkTtWKpcRU/nV7ayuGg/v06dMlSeaZmZn45z//iczMTKvWVyqVOH/+vE3N27ZMLFL25MnSpDK2dCa0ptd2ZdszxL9w4UKUlJTA29vbaQnSng6T1pwIsPZMleEvgjyKpeuOtta87T14mktw9nTMeuWVV3Do0CG88sorVq2vVCoxe/ZsnDt3zuqEZut99FINdFKWIbGNGTNG9H3VhrnZmzdv7pB7zqVOrNacCDCZU2X4iyCPInZQDUeOlFa+B7OzR2Wz9175qj5Tw21aWq3WYiuALYnIkNhSU1NFNz0rlUp07twZGzZscIshU6ta311G6iP5YjInjyL2uqPhYJqZmSn5CFtS1qZSUlLQpk0bpKSkWP2e6OhofPTRR4iOjha1z6o+U8NtWklJSRZbAcR8DnPnzhV9IuKoWqzY7VpK5o4cqY88H5M5eRRLtdDKapmG5WLGT7e2NUCK2pehJ3WHDh1Eb8NWVX2m8+fPNw5iYmldMapqejaXHF09tSybwclVFEI1ORXU6XRQq9UoKCiASqVydTjkICUlJbhw4QJCQ0Ph5eVV4fVu3brhjz/+QMuWLbF79+4Ky0NCQuDr64vly5cjNjbW+Lper4dOpzN7/3Vl2ywfDwBcuHABQUFBuHz5cqUx2lM+S+9RKpWVlkPMPst/LmLiM8ew3Vq1auHixYuitydVPI5WUlKC8+fPAwAiIiLcOlZyDHtzFE8fySPYevtU+ZqjYXlKSorZmmBV26+qNmqupmbPZBn21vykvjbrbk3ZRNUV/1LIIxiSFIAqk0Blzbb23FJkbVOwPU355rZnLWfOtMUkLI4jx9Cn6oG/GvIo9iaTyt4vRUI0dCQDnDsVatkObGLKYctn6m7J3N3iqYytY+gTlcdfDXkEa5NUZQd3w/LKBnaRYnQuS035jlL2MoCYcrgyITrq5MzdSN1xkKofdoAjj1BVJzQpHDhwAJMmTarQMQ6wvpOVYT17OsCJUTY+rVZbaTnIdeTSUY8chx3giOD4mo0U85q7Qy1R6vneicg98BSQPIIhSRluAXMma5O0Yb2SkhInRWa6X6VSyRHGiDwUkzmRFawZKESKbTlC+fhc3TpARNJjMieygq0J21nbktO+qXLucAmG5I3JnDwCD4YkZzzJInsxmZNH4MGQiKozVmOIiIhkjsmciIhI5pjMiYiIZI7JnIiISOaYzImIiGROVsk8JSUFTZo0gZ+fHzp27Ij9+/e7OiQiIiKXk00y/+yzz5CYmIiZM2ciKysL0dHR6N27N65everq0IiIiFxKNsl88eLFeOGFFzB27Fi0bt0aq1atQq1atbBmzRpXh0ZERORSshg05v79+8jMzERycrJxmVKpRM+ePbF3716z7ykqKkJRUZHxeUFBAYC/p5kjIiJyJ4bcZMhVBr6+vvD19bX4flkk82vXrqG0tBSNGjUyWd6oUSP88ccfZt+j0Wgwa9asCsvDwsIcEiMREZG9wsPDTZ7PnDkT77zzjsX3ySKZi5GcnIzExETj8/z8fERERCA3NxdqtdqFkTmOTqdDWFgY8vLyRE1u7+48vXwAy+gpPL2Mnl4+wPllFAQB169fR0BAgMkcE9bUygGZJPPAwEDUqFEDV65cMVl+5coVBAUFmX1PZU0TarXaY398BiqVyqPL6OnlA1hGT+HpZfT08gHOLaM9FU1ZdIDz8fFBhw4dkJaWZlym1+uRlpaGTp06uTAyIiIi15NFzRwAEhMTERcXh4cffhiPPPIIli5ditu3b2Ps2LGuDo2IiMilZJPMhw0bhr/++gszZszA5cuX0a5dO2zfvr1Cp7jK+Pr6YubMmVZff5AjTy+jp5cPYBk9haeX0dPLB8ivjApBEARXB0FERETiyeKaOREREVWOyZyIiEjmmMyJiIhkjsmciIhI5jwume/Zswf9+/dHSEgIFAoFtm3bZvK6IAiYMWMGgoODUbNmTfTs2ROnTp1yTbAiVFW+4uJiTJs2DW3atEHt2rUREhKC0aNH4+LFi64LWARL32FZL7/8MhQKBZYuXeq0+KRgTRmPHz+OAQMGQK1Wo3bt2oiNjUVubq7zgxXJUhkLCwsxYcIEhIaGombNmsYJlORCo9EgNjYW/v7+aNiwIQYNGoQTJ06YrHPv3j3Ex8ejfv36qFOnDoYOHVph8Ct3ZqmMN27cwMSJE9GiRQvUrFkT4eHhmDRpUoXxxd2VNd+hgSAI6Nu3r8Vjkqt4XDK/ffs2oqOjkZKSYvb1+fPnY/ny5Vi1ahUyMjJQu3Zt9O7dG/fu3XNypOJUVb47d+4gKysLb7/9NrKysrBlyxacOHECAwYMcEGk4ln6Dg22bt2Kffv2ISQkxEmRScdSGc+cOYPHH38cLVu2xO7du3HkyBG8/fbb8PPzc3Kk4lkqY2JiIrZv345PP/0Ux48fx5QpUzBhwgR89dVXTo5UnPT0dMTHx2Pfvn348ccfUVxcjF69euH27dvGdRISEvD1119j8+bNSE9Px8WLFzFkyBAXRm0bS2W8ePEiLl68iIULFyI7OxupqanYvn07xo0b5+LIrWPNd2iwdOlSKBQKF0RpJcGDARC2bt1qfK7X64WgoCBhwYIFxmX5+fmCr6+vsHHjRhdEaJ/y5TNn//79AgDh/PnzzglKYpWV8cKFC0Ljxo2F7OxsISIiQliyZInTY5OKuTIOGzZMGDlypGsCcgBzZXzooYeE2bNnmyxr37698NZbbzkxMulcvXpVACCkp6cLgvD3scXb21vYvHmzcZ3jx48LAIS9e/e6Kky7lC+jOZ9//rng4+MjFBcXOzEyaVRWPq1WKzRu3Fi4dOmSVcddV/C4mnlVcnJycPnyZfTs2dO4TK1Wo2PHjpVOpSp3BQUFUCgUqFu3rqtDkYxer8eoUaMwdepUPPTQQ64OR3J6vR7ffvstHnzwQfTu3RsNGzZEx44d3bJpzx6PPfYYvvrqK/z5558QBAG7du3CyZMn0atXL1eHJoqhaTkgIAAAkJmZieLiYpPjTcuWLREeHi7b4035Mla2jkqlgpeXbMYkMzJXvjt37uBf//oXUlJSKp0LxB1Uq2R++fJlADA7larhNU9y7949TJs2DcOHD/eoyRDmzZsHLy8vTJo0ydWhOMTVq1dRWFiIuXPnok+fPvjhhx8wePBgDBkyBOnp6a4OTzIrVqxA69atERoaCh8fH/Tp0wcpKSno0qWLq0OzmV6vx5QpU9C5c2dERUUB+Pt44+PjU+FEWq7HG3NlLO/atWuYM2cOXnzxRSdHZ7/KypeQkIDHHnsMAwcOdGF0lsnv1ImsUlxcjH/+858QBAErV650dTiSyczMxLJly5CVleXe16/soNfrAQADBw5EQkICAKBdu3b47bffsGrVKnTt2tWV4UlmxYoV2LdvH7766itERERgz549iI+PR0hIiEltVg7i4+ORnZ2NX375xdWhOIylMup0OjzzzDNo3bq1VfNvuxtz5fvqq6/w008/QavVujAy61SrmrmhicSWqVTlyJDIz58/jx9//NGjauU///wzrl69ivDwcHh5ecHLywvnz5/Ha6+9hiZNmrg6PEkEBgbCy8sLrVu3NlneqlUrWfVmr8rdu3fx5ptvYvHixejfvz/atm2LCRMmYNiwYVi4cKGrw7PJhAkT8M0332DXrl0IDQ01Lg8KCsL9+/eRn59vsr4cjzeVldHg1q1b6NOnD/z9/bF161Z4e3u7IErxKivfTz/9hDNnzqBu3brG4w0ADB06FN26dXNRtOZVq2QeGRmJoKAgk6lUdTodMjIyPGYqVUMiP3XqFHbu3In69eu7OiRJjRo1CkeOHMGhQ4eMj5CQEEydOhU7duxwdXiS8PHxQWxsbIVbZE6ePImIiAgXRSWt4uJiFBcXQ6k0PQTVqFHD2DLh7gRBwIQJE7B161b89NNPiIyMNHm9Q4cO8Pb2NjnenDhxArm5ubI53lgqI/D3MbRXr17w8fHBV199Jas7LiyVLykpqcLxBgCWLFmCtWvXuiDiynlcM3thYSFOnz5tfJ6Tk4NDhw4hICAA4eHhmDJlCt599100b94ckZGRePvttxESEoJBgwa5LmgbVFW+4OBgPPvss8jKysI333yD0tJS47W5gIAA+Pj4uCpsm1j6DsufoHh7eyMoKAgtWrRwdqiiWSrj1KlTMWzYMHTp0gVPPvkktm/fjq+//hq7d+92XdA2slTGrl27YurUqahZsyYiIiKQnp6Ojz/+GIsXL3Zh1NaLj4/Hhg0b8OWXX8Lf39/4t6ZWq1GzZk2o1WqMGzcOiYmJCAgIgEqlwsSJE9GpUyc8+uijLo7eOpbKaEjkd+7cwaeffgqdTgedTgcAaNCgAWrUqOHK8C2yVL6goCCzrSjh4eFmT2xcyqV96R1g165dAoAKj7i4OEEQ/r497e233xYaNWok+Pr6Cj169BBOnDjh2qBtUFX5cnJyzL4GQNi1a5erQ7eape+wPDnemmZNGVevXi088MADgp+fnxAdHS1s27bNdQGLYKmMly5dEsaMGSOEhIQIfn5+QosWLYRFixYJer3etYFbqbK/tbVr1xrXuXv3rvDqq68K9erVE2rVqiUMHjxYuHTpkuuCtpGlMlb2HQMQcnJyXBq7Naz5Ds29xx1vTeMUqERERDJXra6ZExEReSImcyIiIpljMiciIpI5JnMiIiKZYzInIiKSOSZzIiIimWMyJyIikjkmcyIiIpljMiciIpI5JnMiIiKZYzInIpLA9evX0bBhQ5w7d86m9z333HNYtGiRY4KiaoPJnKqVMWPGQKFQQKFQwNvbG40aNcJTTz2FNWvWSDr1Zrdu3TBlypQq1zHEUdnjnXfekSwed2PN5+NIO3bssPj5//DDDzZt87333sPAgQPRpEkT47IxY8ZUmJHxiy++gJ+fnzGBT58+He+99x4KCgrsLRZVY0zmVO306dMHly5dwrlz5/D999/jySefxOTJk9GvXz+UlJQ4LY5Lly4ZH0uXLoVKpTJZ9vrrrzstFqncv39fFvvr0qWLyWddv359vP322ybLevToYfX27ty5g9WrV2PcuHFVrvef//wHI0aMwMqVK/Haa68BAKKiotCsWTN8+umnospCBMDzpkAlqkpcXJwwcODACsvT0tIEAMJHH30kCIL5aVWjo6OFmTNnGp9v3rxZiIqKEvz8/ISAgAChR48eQmFhoRAXF2fzdJBr164V1Gq12ddKS0uF999/X2jSpIng5+cntG3bVti8ebPx9a5duwoTJkwQJk+eLNStW1do2LCh8OGHHwqFhYXCmDFjhDp16gjNmjUTvvvuO5Ptdu3aVYiPjxfi4+MFlUol1K9fX5g+fbpxClJL+y27jcmTJwv169cXunXrJgiCIHz//fdC586dBbVaLQQEBAjPPPOMcPr0aeN3UNnnY+lzr2x/1sRamQsXLggAhB9++MGq9c3ZvHmz0KBBgwrLy/7e5s2bJ/j5+QlbtmypsN6sWbOExx9/XPT+iVgzJwLQvXt3REdHY8uWLVatf+nSJQwfPhzPP/88jh8/jt27d2PIkCEQBAHLli1Dp06d8MILLxhreWFhYaJj02g0+Pjjj7Fq1SocO3YMCQkJGDlyJNLT043rrFu3DoGBgdi/fz8mTpyIV155Bf/4xz/w2GOPISsrC7169cKoUaNw584dk22vW7cOXl5e2L9/P5YtW4bFixfjP//5j9X7NWzDx8cHv/76K1atWgUAuH37NhITE3Hw4EGkpaVBqVRi8ODB0Ov1dn8+5vZnbazmaLVaAED79u2tjqG8n3/+GR06dKj09WnTpmHOnDn45ptvMHjw4AqvP/LII9i/fz+KiopEx0DVnKvPJoicqbKauSAIwrBhw4RWrVoJgmC5hpiZmSkAEM6dO2d2W127dhUmT55sdVyV1czv3bsn1KpVS/jtt99Mlo8bN04YPny4cV9la3UlJSVC7dq1hVGjRhmXXbp0SQAg7N271yTGVq1aGWvigiAI06ZNE1q1amXVfg3biImJsVi+v/76SwAgHD161Pg+c5+PNTXz8vuzNtbKzJo1SwgLC7O4Xo8ePYSTJ0+afW3gwIHC888/X2F5XFyc4OPjIwAQ0tLSKt324cOHq/w9EVni5dIzCSI3IggCFAqFVetGR0ejR48eaNOmDXr37o1evXrh2WefRb169SSN6fTp07hz5w6eeuopk+X3799HTEyM8Xnbtm2N/69Rowbq16+PNm3aGJc1atQIAHD16lWT7Tz66KMmZe7UqRMWLVpk9X4BmK2Rnjp1CjNmzEBGRgauXbtm7FyYm5uLqKgoq8pemfL7syVWc7KysqyqlZ86dQrNmjUz+9rdu3fh5+dn9rW2bdvi2rVrmDlzJh555BHUqVOnwjo1a9YEgAotJ0TWYjIn+j/Hjx9HZGQkAECpVEIQBJPXi4uLjf+vUaMGfvzxR/z222/44YcfsGLFCrz11lvIyMgwbkMKhYWFAIBvv/0WjRs3NnnN19fX+H9vb2+T1wy99cs+B2B1j31r9wsAtWvXrvD+/v37IyIiAh999BFCQkKg1+sRFRVlscOapc/d3P5sidWcrKwsjB8/vsLyY8eO4fnnn8fdu3cxbNgwBAUFQak0f2UyMDAQN2/eNPta48aN8cUXX+DJJ59Enz598P3338Pf399knRs3bgAAGjRoYDFeInOYzIkA/PTTTzh69CgSEhIA/H1QvXTpkvF1nU6HnJwck/coFAp07twZnTt3xowZMxAREYGtW7ciMTERPj4+KC0ttTuu1q1bw9fXF7m5uejatavd2ysvIyPD5Pm+ffvQvHlzu/Z7/fp1nDhxAh999BGeeOIJAMAvv/xisk5ln481n3t59sR67do15OXlVaiZFxUVYdiwYdi4cSPatGmDQYMGmbR+lBcTE1Nlb/SIiAikp6cbE/r27dtNEnp2djZCQ0MRGBhoU/xEBkzmVO0UFRXh8uXLKC0txZUrV7B9+3ZoNBr069cPo0ePBvB3h7jU1FT0798fdevWxYwZM1CjRg3jNjIyMpCWloZevXqhYcOGyMjIwF9//YVWrVoBAJo0aYKMjAycO3cOderUQUBAQKW1uqr4+/vj9ddfR0JCAvR6PR5//HEUFBTg119/hUqlQlxcnF2fRW5uLhITE/HSSy8hKysLK1aswKJFi+zab7169VC/fn18+OGHCA4ORm5uLpKSkkzWqezzsfS5S/0ZZWVlAajY+W3btm3o2rWr8VJFq1atEBISUul2evfujeTkZNy8ebPSSy1hYWHYvXs3nnzySfTu3Rvbt2+HSqUC8HcHul69elVZTqKqMJlTtbN9+3YEBwfDy8sL9erVQ3R0NJYvX464uDhjwk1OTkZOTg769esHtVqNOXPmmNQQVSoV9uzZg6VLl0Kn0yEiIgKLFi1C3759AQCvv/464uLi0Lp1a9y9exc5OTkmg4nYYs6cOWjQoAE0Gg3Onj2LunXron379njzzTft/ixGjx6Nu3fv4pFHHkGNGjUwefJkvPjii3btV6lUYtOmTZg0aRKioqLQokULLF++HN26dTOuU9nnY+lzr4zYWLVaLRo1alQhUR89ehTt2rUzPs/MzESfPn0q3U6bNm3Qvn17fP7553jppZcqXS80NNQkoe/YsQM+Pj7Ytm0btm/fbrGcRJVRCOUvUBFRtdCtWze0a9cOS5cudXUobmfx4sU4f/48li1bhh9++AF9+vTB9evXq+zg+O2332Lq1KnIzs62qRVm5cqV2Lp1q80jzhGVxZo5EVE5I0eORN++fRETE4OoqChERERYvFPhmWeewalTp/Dnn3/adN+8t7c3VqxYYW/IVM2xZk5UTbFmTuQ5mMyJiIhkjsO5EhERyRyTORERkcwxmRMREckckzkREZHMMZkTERHJHJM5ERGRzDGZExERyRyTORERkcwxmRMREckckzkREZHM/X9oPpINTTTX1gAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the dust column density\n",
    "PvT_plot_Step2 = HtoH.PvT(CatP, CatT, Tscale=[10,25], Pscale=[0.0,30.0], CatdP=CatdP, errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "491d0620",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Saving figures\n",
    "# PNG\n",
    "PvI_plot_Step2.savefig('Plots/'+SourceName+'_PvI_plot_Step2.png',dpi=300)\n",
    "PvN_plot_Step2.savefig('Plots/'+SourceName+'_PvN_plot_Step2.png',dpi=300)\n",
    "PvT_plot_Step2.savefig('Plots/'+SourceName+'_PvT_plot_Step2.png',dpi=300)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "PoLiteWIP",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
